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Abstract 

Different theoretical methods used for the description of diffractive processes in 
small-a; deep inelastic scattering are reviewed. The semiclassical approach, where a par- 
tonic fluctuation of the incoming virtual photon scatters off a superposition of target 
colour fields, is used to explain the basic physical effects. In this approach, diffraction 
occurs if the emerging partonic state is in a colour singlet, thus fragmenting indepen- 
dently of the target. Other approaches, such as the idea of the pomeron structure function 
and two gluon exchange calculations, are also discussed in some detail. Particular atten- 
tion is paid to the close relation between the semiclassical approach and the method of 
diffractive parton distributions, which is linked to the relation between the target rest 
frame and the Breit frame point of view. While the main focus is on diffractive structure 
functions, basic issues in the diffractive production of mesons and of other less inclusive 
final states are also discussed. Models of the proton colour field, which can be converted 
into predictions for diffractive cross sections using the semiclassical approach, are pre- 
sented. The concluding overview of recent experimental results is very brief and mainly 
serves to illustrate implications of the theoretical methods presented. 
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1 Introduction 



1.1 Preface 

The term diffraction is derived from optics, where it describes the deflection of a beam 
of hght and its decomposition into components with different frequencies. In high energy 
physics it was originally used for small-angle elastic scattering of hadrons. If one of the 
hadrons, say the projectile, is transformed into a set of two or more final state particles, 
the process is called diffractive dissociation or inelastic diffraction. Good and Walker 
have pointed out that a particularly intuitive physical picture of such processes emerges 
if the projectile is described as a superposition of different components which scatter 
elastically off the target [Q]. Since the corresponding elastic amplitude is different for 
each component, the outgoing beam will contain a new superposition of these components 
and therefore, in general, new physical states. These are the dissociation products of the 
projectile. 

Even at very high energy, the above processes are generically soft, i.e., the momentum 
transfer is small and the dissociation products have small p±. Therefore, no immediate 
relation to perturbative QCD is apparent. 

By contrast, diffractive jet production, observed at the CERN SppS collider in proton- 
antiproton collisions p[, involves a hard scale. Although one of the hadrons escapes 
essentially unscathed, a high-p^ jet pair, which is necessarily associated with a high 
virtuality in the intermediate states, is produced in the central rapidity range. The 
cross section of the process is parametrically unsuppressed relative to non-diffractive jet 
production. This seems to contradict a naive partonic picture since the colour neutrality 
of the projectile is destroyed if one parton is removed to participate in the hard scattering. 
The interplay of soft and hard physics necessary to explain the effect provides one of the 
main motivations for the study of these 'hard diffractive' processes. 

The present review is focussed on diffraction in deep inelastic scattering (DIS), which 
is another example of a hard diffractive process. This process became experimentally 
viable with the advent of the electron-proton collider HERA, where DIS at very small 
values of the Bjorken variable x can be studied. In the small-x or high-energy region, 
a significant fraction of the observed DIS events have a large rapidity gap between the 
photon and the proton fragmentation region 0,^. In contrast to the standard DIS process 
7*p — ^ X, the relevant reaction reads 7*p —>■ XY, where X is a high-mass hadronic state 
and Y is the elastically scattered proton or a low- mass excitation of it. Again, these events 
are incompatible with the naive picture of a partonic target and corresponding simple 
ideas about the colour flow. Naively, the parton struck by the virtual photon destroys the 
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colour neutrality of the proton, a colour string forms between struck quark and proton 
remnant, and hadronic activity is expected throughout the detector. Nevertheless, the 
observed diffractive cross section is not power suppressed at high virtualities with 
respect to standard DIS. 

The main theoretical interest is centered around the interplay of soft and hard physics 
represented by the elastic or almost elastic scattering of the proton and the scattering 
of the highly virtual photon respectively. Diffractive DIS is much simpler than hard 
diffraction in hadronic reactions since only one non-perturbative object is involved. In a 
large fraction of the events, the momentum transfer to the proton is very small. Therefore, 
one can hope to gain a better understanding of the bound state dynamics of the proton by 
studying diffractive DIS. Inclusive reactions of the virtual photon with hadronic targets 
are well-studied theoretically and constrained by a large amount of DIS data. It is a 
challenge to utilize this knowledge for the investigation of the non-perturbative dynamics 
of diffraction and, thereby, of the proton structure. 

The paper is organized as follows. In the remainder of the Introduction, different 
approaches to diffractive DIS are put into historical and physical perspective. The semi- 
classical model, which is particularly close to the interests of the author, is emphasized. 
Chapter ^ is concerned with the fundamental observations, the basic concepts required 
for their understanding, and the necessary kinematic considerations. The semiclassical 
approach, which is used as the starting point for the discussion of other models, is intro- 
duced in Chapter ^. In Chapter ^, the concepts of soft pomeron, triple pomeron vertex 
and pomeron structure function are discussed. Diffractive parton distributions, which are 
the more fundamental objects from the point of view of perturbative QCD, are intro- 
duced. The semiclassical approach is put in relation to the method of diffractive parton 
distributions, and explicit formulae for these quantities are derived. In Chapter |^, two 
gluon exchange calculations are discussed. The emphasis is on their rigorous validity in 
specific kinematic situations and on their partial correspondence to the semiclassical ap- 
proach. Chapter ^ introduces three models for the colour field of the proton, which can be 
converted into predictions for diffractive cross sections using the methods of Chapters § 
and 1^. A discussion of recent experimental results and their description by theoretical 
models is given in Chapter ^, followed by the Conclusions. 

1.2 Models for diffraction 

The claim that diffraction in DIS should be a leading twist effect, and the understanding 
of the fundamental mechanism underlying such processes can be traced back to the 
famous paper of Bjorken and Kogut Their argument is based on a qualitative picture 
of DIS in the target rest frame, where the incoming virtual photon can be considered as 
a superposition of partonic states. The large virtuality sets the scale, so that states 
with low-px partons, i.e., aligned configurations, are suppressed in the photon wave 
function. However, in contrast to high-p^ configurations, these aligned states have a large 
hadronic interaction cross section with the proton. Therefore, their contribution to the 
DIS cross section is expected to be of leading twist. Naturally, part of this leading twist 
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contribution is difFractive since the above \ow-p± configurations represent transversely 
extended, hadron-like objects, which have a large elastic cross section with the proton. 

Note that the basic technical methods date back even further, namely, to the calcu- 
lation of yU^/i~ pair electroproduction off an external electromagnetic field by Bjorken, 
Kogut and Soper They derive the transition amplitude of the incoming virtual pho- 
ton into a /i^/U^ pair, which corresponds to the qq wave function of the virtual photon 
employed by Nikolaev and Zakharov in their seminal work on diffraction at HERA |^. 

The above physical picture, commonly known as the aligned jet model of diffractive 
and non-diffractive DIS, can naturally be extended to account for a soft energy growth 
of both processes. (Note that we are dealing with the underlying energy dependence of 
the soft scattering process, not with the additional x-dependence induced by Altarelli- 
Parisi evolution, that is added on top of it.) Since a large sample of hadronic cross 
sections can be consistently parametrized using the concept of the Donnachie-Landshoff 
or soft pomeron 0, it is only natural to assume that this concept can also be used to 
describe the energy dependence of the interaction of the aligned jet component of the 
small-x virtual photon and the target. This gives rise to the desired non-trivial energy 
dependence although, as will be discussed in more detail later on, this energy dependence 
is not sufficiently steep in the diffractive case. Note, however, that the above soft energy 
dependence may also be logarithmic, which has the advantage of explicit consistency 
with unitarity at arbitrarily high energies. 

A more direct way of applying the concept of the soft pomeron to the phenomenon of 
hard diffraction was suggested by Ingelman and Schlein in the context of diffractive jet 
production in hadronic collisions Their idea of a partonic structure of the pomeron, 
which can be tested in hard processes, applies to the case of diffractive DIS as well [llO| . 
Essentially, one assumes that the pomeron can, like a real hadron, be characterized by 
a parton distribution. This distribution factorizes from the pomeron trajectory and the 
pomeron-proton vertex, which are both obtained from the analysis of purely soft hadronic 
reactions. The above non-trivial assumptions are often referred to as 'Regge hypothesis' 
or 'Regge factorization'. 

The Ingelman- Schlein approach described above is based on the intuitive picture of a 
pomeron flux associated with the proton beam and on the conventional partonic descrip- 
tion of the pomeron photon collision. In the limit where not only the total proton-photon 
center-of-mass energy but also the energy of the pomeron-photon collision becomes large, 
the concept of the triple pomeron vertex can be applied |TI|]. At HERA, this limit cor- 
responds to rapidity gap events with very large diffractive masses. 

The concept of fracture functions of Veneziano and Trentadue or, more specifi- 
cally, the diffractive parton distributions of Berera and Soper |[T^ provide a framework 
for the study of diffractive DIS that is firmly rooted in perturbative QCD. Loosely speak- 
ing, diffractive parton distributions describe the probability of finding, in a fast moving 
proton, a parton with a certain momentum fraction x, under the additional requirement 
that the proton remains intact losing only a certain fraction of its momentum. This idea 
is closely related to the concept of a partonic pomeron described above, but it gives up 
the Regge hypothesis, thus being less predictive. This weaker factorization assumption 
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is sometimes referred to as 'diffractive factorization', as opposed to the stronger Regge 
factorization assumption. 

An essential feature of diffractive DIS is the colour singlet exchange in the t channel, 
necessary to preserve the colour neutrality of the target. None of the above approaches 
address this requirement explicitly within the framework of QCD. Instead, the colour 
singlet exchange is postulated by assuming elastic scattering in the aligned jet model, by 
using the concept of the pomeron, or by defining diffractive parton distributions. If one 
insists on using the fundamental degrees of freedom of QCD, the simplest possibility of 
realizing colour singlet exchange at high energy is the exchange of two gluons |T^. One 
could say that two gluons form the simplest model for the pomeron. 

The fundamental problem with this approach is the applicability of perturbative 
QCD to the description of diffractive DIS. As will be discussed in more detail below, the 
hard scale of the photon does not necessarily justify a perturbative description of the 
t channel exchange with the target. 

Diffractive processes where the t channel colour singlet exchange is governed by a 
hard scale include the electroproduction of heavy vector mesons , electroproduction 
of light vector mesons in the case of longitudinal polarization or at large t ||17|| , 
and virtual Compton scattering |]TB|-pU[. In the leading logarithmic approximation, the 



relevant two-gluon form factor of the proton can be related to the inclusive gluon dis- 
tribution |15|. Accordingly, a very steep energy dependence of the cross section, which 
is now proportional to the square of the gluon distribution, is expected from the known 
steep behaviour of small-x structure functions. The origin of this steep rise itself may 
be attributed to a combination of the x-dependence of the input distributions and their 
Altarelli-Parisi evolution, or to the BFKL resummation of large logarithms of x. 

To go beyond leading logarithmic accuracy, the non-zero momentum transferred 
to the proton has to be taken into account. This requires the use of 'non- forward' or 



'off-diagonal' parton distributions (see [0,|2TI] and refs. therein), which were discussed 
in 1 19, 20] within the present context. Although their scale dependence is predicted by 
well-known evolution equations, only limited information about the relevant input dis- 
tributions is available. In particular, the simple proportionality to the square of the 
conventional gluon distribution is lost. 

The perturbative calculations of meson electroproduction discussed above were put 



on a firmer theoretical basis by the factorization proof of where the conditions 
required for the applicability of perturbation theory are discussed in detail. As shown 
explicitly in the simple model calculation of ||23|, QCD gauge invariance ensures that 



the non-perturbative meson formation process takes place after the scattering off the 
hadronic target. 

The situation becomes even more complicated if two gluon exchange calculations are 
applied to more general diffractive final states. Examples are the exclusive electroproduc- 
tion of heavy quark pairs p^- p7|] or high-p_|_ jets p8|-pT|. A straightforward perturbative 



analysis shows that the virtual photon side of the process is dominated by small trans- 
verse distances, establishing two gluon exchange as the dominant mechanism. However, 
it is not obvious whether the very definition 'exclusive', meaning exclusive on a partonic 
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level, can be extended to all orders in perturbation theory. As will become clear later 
on, additional soft partons in the final state may destroy the hardness provided by the 
large p± or the heavy quark mass. 

From the point of view of the full diffractive cross section, perturbative two gluon 
exchange corresponds to a higher twist effect. However, it can be argued that it dominates 
the kinematic region where the mass M of the diffractively produced final state X is 
relatively small, M <^ Q [B^]. More generally, higher twist effects in diffraction, their 



calculability and possibilities for their experimental observation are discussed in [^^-|34|. 

Eventually it is possible to attempt the description of the full cross section of diffrac- 
tive DIS within the framework of two gluon exchange. In such a general setting, pertur- 
bation theory can not be rigorously justified and two gluons are basically used as a model 
for the soft colour singlet exchange in the t channel. Nevertheless, it is interesting to see 
to what extent concepts and formulae emerging from perturbation theory can describe a 
wider range of phenomena. In particular, much attention has been devoted to the energy 
dependence of diffractive cross sections both in the framework of conventional BFKL 
summation of gluon ladders P5|,|55| and within the colour dipole approach to small-x 



resummation |P8|,|37||. 

An apparently quite different approach emerged with the idea that soft colour in- 
teractions might be responsible for the large diffractive cross section at HERA. The 
starting point is Buchmiiller's observation of the striking similarity between x and 
dependence of diffractive and inclusive DIS at small x |3^. It is then natural to assume 
that the same hard partonic processes underlie both cross sections and that differences 
in the final states are the result of non-perturbative soft colour exchange 



In Ref. [^], boson-gluon fusion was proposed as the dominant partonic process, and 



diffraction was claimed to be the result of soft colour neutralization of the produced qq 
pair. A simple statistical assumption about the occurrence of this colour neutralization 
lead to a surprisingly good fit to the data. Note that similar ideas concerning the rotation 



of quarks in colour space had previously been discussed by Nachtmann and Reiter |5T 
in connection with QCD vacuum effects on hadron-hadron scattering. 

A closely related approach was introduced in [^, where the assumption of soft colour 
neutralization was implemented in a Monte Carlo event generator based on perturbation 
theory. Normally, the partonic cascades underlying the Monte Carlo determine the colour 
of all partons produced and, with a certain model dependence, the hadronic final state. 
The introduction of an ad-hoc probability for partons within the cascade to exchange 
colour in a non-perturbative way leads to a significant increase of final states with rapidity 
gaps, resulting in a good description of the observed diffractive events. 

Both of the above models share a fundamental theoretical problem: the soft colour 
exchange is introduced in an ad-hoc manner, independently of kinematic configuration 
or space-time distances of the partons involved. If all relevant distances are large, our 
ignorance of the true mechanism of non-perturbative interactions justifies the naive as- 
sumption of random colour exchange. However, in many other situations, e.g., if a small 
colour dipole is involved, colour exchange is suppressed in a perturbatively calculable 
way. This is known as colour transparency. The above simple soft colour models do not 
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properly account for this effect. One possibility of incorporating the knowledge of the 
perturbative aspects of QCD while keeping the essential idea of soft colour exchange as 
the source of leading twist diffraction is the semi classical treatment. 



1.3 The semiclassical method and its relation to 
other approaches 

The semiclassical approach to diffractive electroproduction evolved from the attempt to 



justify the phenomenologically successful idea of soft colour exchange ^ within 
the framework of QCD. The starting point is the proton rest frame picture of DIS 
advocated long ago in and developed since by many authors (see in particular 0]). 
To describe diffraction within this framework, the scattering of energetic partons in 
the photon wave function off the hadronic target needs to be understood. While the 
scattering of small transverse size configurations is calculable perturbatively, i.e., via two 
gluon exchange, different models may be employed for larger size configurations. A rather 
general framework is provided by the concept of the dipole cross section a{p) utilized in 
the analysis of [0. 

It is helpful to begin by recalling the discussion of high-energy electroproduction 
of fi~^fi~ pairs off atomic targets given by Bjorken, Kogut and Soper in [^]. There, the 
closely related QED problem was solved by treating the high-energy scattering of fj,~^fj,~ 
pairs off a given electromagnetic field in the eikonal approximation. An analogous QCD 



process was considered by Collins, Soper and Sterman in where the production of 



heavy quark pairs in an external colour field was calculated. In Ref. |^^, Nachtmann 
developed the idea of eikonalized interactions with soft colour field configurations as a 
method for the treatment of elastic high-energy scattering of hadrons. 



The semiclassical approach to diffraction, introduced in combines the concepts 



and methods outlined above. The proton is modelled by a soft colour field, and the 
interactions of the fast partons in the photon wave function with this field are treated in 
the eikonal approximation. As a result, the non-perturbative proton structure is encoded 
in a combination of non-Abelian phase factors associated with the partons. In the simplest 
case of a qq fluctuation, two such phase factors are combined in a Wegner- Wilson loop, 
which carries all the information about the target. Since the outgoing qq pair can be 
either in a colour singlet or in a colour octet state, both diffractive and non-diffractive 
events are naturally expected in this approach. Thus, both types of events are described 
within the same framework. The colour state of the produced pair is the result of the 
fundamentally non-perturbative interaction with the proton colour field, very much in 



the spirit of the soft colour neutralization of |]3^ . 

The essential difference between the semiclassical calculation and the soft colour 
proposal is the recognition that the possibility of soft colour exchange is intimately 
related to parton level kinematics. The semiclassical approach allows for a consistent 
treatment of the scattering of both small and large transverse size configurations. The 
requirement of colour neutrality in the final state of diffractive events is explicitly shown 
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to suppress the former ones. This confirms the famihar quahtative arguments on which 
the ahgned jet model is based 0. 

A further essential step is the inclusion of higher Fock states in the photon wave 
function. In the framework of two gluon exchange, corresponding calculations for the qqg 
state were performed in I^Sj. The semiclassical analysis of the qqg state demonstrates 



that, to obtain leading twist diffraction, at least one of the three partons has to have 
small p± and has to carry a small fraction of the photon's longitudinal momentum. This 
is a natural generalization of the aligned jet model, which was previously discussed on a 
qualitative level in P^]. 

A fundamental prediction derived from the semiclassical treatment of the qqg state 
is the leading twist nature of diffractive electroproduction of heavy quark pairs [0 and 
of high-px jets These large cross sections arise from qqg configurations where the 



gluon is relatively soft and has small p±, so that the t channel colour singlet exchange 
remains soft in spite of the additional hard scale provided by the qq pair. 

As mentioned previously, considerable work has been done attempting to employ the 
additional hardness provided by final states with jets or heavy quarks in order to probe 
mechanisms of perturbative colour singlet exchange [p^ -|3T|. In fact, the semiclassical 
approach is also well suited to follow this line of thinking. As demonstrated in |^ , the 



features of diffractive jet production from qq and qqg fluctuations of the photon, both 
treated consistently in the semiclassical framework, are very different. In the qq case, a 
small-size colour dipole tests the small distance structure of the proton colour field. It can 
be shown within the semiclassical framework that the relevant operator is the same as 
in inclusive DIS, thus reproducing the well-known relation with the square of the gluon 



distribution [15|. In the qqg case, the additional soft gluon allows for a non-perturbative 



colour neutralization mechanism even though the final state contains two high-p^ jets. 

The semiclassical treatment of leading twist diffraction is equivalent to a treatment 
based on diffractive parton distributions [^. To be more precise, the results of lead to 



the conclusion that a leading twist contribution to diffraction can arise only from virtual 
photon fluctuations that include at least one soft parton. Starting from this premise, 
the diffractive production of final states with one soft parton and a number of high-p^ 
partons was considered in |^9[. It was shown that the cross section can be written as a 



convolution of a hard partonic cross section and a diffractive parton distribution. The 
latter is given explicitly in terms of an average over the target colour field configurations 
underlying the semiclassical calculation. 

From this point of view, the semiclassical prediction of the dominance of qqg final 
states and soft colour neutralization in high-p^ jet and heavy quark production appears 
natural. It corresponds simply to the dominant partonic process, boson-gluon fusion, 
with the gluon taken from the diffractive gluon distribution. The additional gluon in the 
final state is, from the point of view of parton distributions, merely a consequence of 
the preserved colour neutrality of the target. From the point of view of the semiclassical 
calculation, it is the scattering product of the original gluon from the photon wave 
function that is necessary for the softness of the colour singlet exchange. 

Even without an explicit model for the proton colour field, the semiclassical approach 
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provides an intuitive overall picture of inclusive and diffractive small-x DIS and predicts a 
number of qualitative features of the final state. However, the ultimate goal is to develop 
the understanding of non-perturbative colour dynamics that is needed to calculate the 
required field averages. 

A very successful description of a number of hadronic reactions was obtained in the 
framework of the stochastic vacuum. Different effects of the non-trivial structure of the 
QCD vacuum in high-energy hadron-hadron scattering were considered in [^. Field 



strength correlators in a theory with vacuum condensate were employed in |5y,^ to de- 
scribe the pomeron. The model of the stochastic vacuum was introduced by Dosch and 



Simonov |^2[ in the context of gluon field correlators in the Euclidean theory. Based 



on ideas of [50,^ and, in particular, on Nachtmann's suggestion to describe high- 



energy hadron-hadron scattering in an eikonal approach the stochastic vacuum 



model was adapted to Minkowski space and applied to high-energy processes by Dosch 
and Kramer [Q. In the present context of diffractive and non-diffractive electroproduc- 
tion, the work on inclusive DIS [0, C- and P-odd contributions to diffraction |5^, and 
exclusive vector meson production pB|-p^ is particularly relevant. Results on diffractive 
structure functions have recently been reported ||59||. 



In the special case of a very large hadronic target, the colour field averages required 
in the semiclassical approach can be calculated without specifying the details of the non- 
perturbative colour dynamics involved. McLerran and Venugopalan observed that the 
large target size, realized, e.g., in an extremely heavy nucleus, introduces a new hard 



scale into the process of DIS |Q . From the target rest frame point of view, this means 
that the typical transverse size of the partonic fluctuations of the virtual photon remains 
perturbative Thus, the perturbative treatment of the photon wave function in the 



semiclassical calculation is justified. Note that the small size of the partonic fluctuations 
of the photon does not imply a complete reduction to perturbation theory. The long 
distance which the partonic fluctuation travels in the target compensates for its small 
transverse size, thus requiring the eikonalization of gluon exchange. 

For a large target it is natural to introduce the additional assumption that the gluonic 
fields encountered by the partonic probe in distant regions of the target are not corre- 
lated. In this situation, inclusive and diffractive DIS cross sections become completely 
calculable. A corresponding analysis of inclusive and diffractive parton distributions was 
performed in [^|. Starting from the above large target model, expressions for the in- 



clusive quark distribution, which had previously been discussed in a similar framework 
53| , and for the inclusive gluon distribution were obtained. Diffractive quark and gluon 



m 



distributions were calculated, within the same model for the colour field averaging, on 
the basis of formulae from |^9|. The resulting structure functions, obtained from the 



above input distributions by Altarelli-Parisi evolution, provide a satisfactory description 
of the experimental data . 

A number of further approaches to diffractive DIS were proposed by several authors. 
For example, it was argued in that a super-critical pomeron with colour-charge parity 



Cc = —1 plays an essential role. As a result, the main features of the naive boson-gluon 
fusion model of are reproduced while its problems, which were outlined above, are 



avoided. In a different approach, a geometric picture of diffraction, based on the idea 
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of colourless gluon clusters, was advertised (see, e.g., There can be no doubt that 

far more material relevant to the present subject exists which, because of the author's 
bias or unawareness, is not appropriately reflected in this introduction or the following 
chapters. 
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2 Basic Concepts and Phenomena 



In this chapter, the kinematics of diffractive electroproduction at small x is explained in 
some detail, and the notation conventionally used for the description of this phenomenon 
is introduced. The main experimental observations are discussed, and the concept of the 
diffractive structure function, which is widely used in analyses of inclusive diffraction, 
is explained. This standard material may be skipped by readers familiar with the main 
HERA results on diffractive DIS. 



2.1 Kinematics 

To begin, recall the conventional variables for the description of DIS. An electron with 
momentum k collides with a proton with momentum P. In neutral current processes, 
a photon with momentum q and virtuality = —Q^ is exchanged, and the outgoing 
electron has momentum k' = k — q. In inclusive DIS, no questions are asked about 
the hadronic final state Xw, which is only known to have an invariant mass square 
W'^ = (P + g)^. The Bjorken variable x = Q'^/{Q'^ + W'^) characterizes, in the naive 
parton model, the momentum fraction of the incoming proton carried by the quark 
that is struck by the virtual photon. If a; <C 1, which is the relevant region in the present 
context, Q is much smaller than the photon energy in the target rest frame. In this sense, 
the photon is almost real even though ^ A^ (where A is some soft hadronic scale). 
It is then convenient to think in terms of a high-energy '~f*p collision with centre-of-mass 
energy W. 

The small-x limit of DIS became experimentally viable only a few years ago, with 
the advent of the HERA accelerator (Hadron-Elektron-Ring-Anlage) at the DESY accel- 
erator centre in Hamburg. In this machine, 27 GeV electrons (or positrons) are collided 
with 920 GeV protons at interaction points inside the HI and ZEUS detectors. The large 
centre-of-mass energy of the ep collision, 300 GeV, allows for a very high hadronic 

energy W and thus for the observation of events with both very small x and high Q"^. 
For values of well in the perturbative domain, a statistically usable number of events 
with X as small as 10~^ is observed. 

Loosely speaking, diffraction is the subset of DIS characterized by a quasi-elastic in- 
teraction between virtual photon and proton. A particularly simple definition of diffrac- 
tion is obtained by demanding that, in the 7*p collision, the proton is scattered elastically. 
Thus, in diffractive events, the final state contains the scattered proton with momentum 



P' and a diffractive hadronic state Xjv/ with mass M (see Fig. p.ll) . Since diffractive 



events form a subset of DIS events, the total invariant mass of the outgoing proton and 
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Figure |2.1| : Diffractive electroproduction. The full hadronic final state with invariant 
mass W contains the elastically scattered proton and the diffractive state with invariant 
mass M. 



the diffractive state Xm is given by the standard DIS variable W . 

The following parallel description of inclusive and diffractive DIS suggests itself. 

In the former, virtual photon and proton collide to form a hadronic state Xy/ with 
mass W . The process can be characterized by the virtuality and the scaling variable 
X = Q^/{Q^ + W^), the momentum fraction of the struck quark in the naive parton 
model. 

In the latter, a colour neutral cluster of partons is stripped off the proton. The virtual 
photon forms, together with this cluster, a hadronic state Xm with mass M. The process 
can be characterized by Q^, as above, and by a new scaling variable P = Q'^ /{Q"^ + M^), 
the momentum fraction of this cluster carried by the struck quark. 

Since diffraction is a subprocess of inclusive DIS, the struck quark from the colour 
neutral cluster also carries a fraction x of the proton momentum. Therefore, the ratio 
^ = x/ (3 characterizes the momentum fraction that the proton loses to the colour neutral 
exchange typical of an elastic reaction. This exchanged colour neutral cluster loses a mo- 
mentum fraction (3 to the struck quark that absorbs the virtual photon. As expected, the 
product X = (3^is the fraction of the original proton's momentum carried by this struck 
quark. Since the name pomeron is frequently applied to whichever exchange with vac- 
uum quantum numbers dominates the high-energy limit, many authors use the notation 
X]p = ^, thus implying that the proton loses a momentum fraction ^ to the exchanged 
pomeron. 

Therefore, x, and (3 or, alternatively, x, and ^ are the main kinematic variables 
characterizing diffractive DIS. A further variable, t = {P — P')^, is necessary if the 
transverse momenta of the outgoing proton and the state Xm relative to the 7*P axis 
are measured. Since the proton is a soft hadronic state, the value of \t\ is small in most 
events. The small momentum transferred by the proton also implies that M <^ W. 

To see this in more detail, introduce light-cone co-ordinates g± = go ± Is and q± = 
(gi, ^2)- It is convenient to work in a frame where the transverse momenta of the incoming 
particles vanish, q± = P± = 0. Let A be the momentum transferred by the proton. 
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A = P — P', and 171"^ = = P'^ the proton mass squared. For forward scattering, 
P'j_ = 0, the relation 

t = A' = A^A^ = -eml (2.1) 

holds. Since ^ = (Q^ + M^)/((5^ + W^), this means that small M implies small \t\ and 
vice versa. Note, however, that the value of |t| is larger for non- forward processes, where 
t = A+A_ - Ai. 

So far, diffractive events have been characterized as those DIS events which contain 
an elastically scattered proton in their hadronic final state. An even more striking feature 
is the large gap of hadronic activity seen in the detector between the scattered proton 
and the diffractive state Xm- It will now be demonstrated that this feature, responsible 
for the alternative name 'rapidity gap events', is a direct consequence of the relevant 
kinematics. 

Recall the definition of the rapidity y of a particle with momentum k, 

y=-ln^ = -ln^^. (2.2) 

This is a convenient quantity for the description of high-energy collisions along the z- 
axis. Massless particles moving along this axis have rapidity —00 or +00, while all other 
particles are characterized by some finite intermediate value of y (for a detailed discussion 
of the role of rapidity in the description of diffractive kinematics see, e.g., ||66|| . ) 

In the centre-of-mass frame of the 7*p collision, with the z-axis pointing in the proton 
beam direction, the rapidity of the incoming proton is given by t/p = ln(P+/mp). At small 
^, the rapidity of the scattered proton is approximately the same. This is to be compared 
with the highest rapidity y^ax of any of the particles in the diffractive state Xm- Since 
the total plus component of the 4-momentum of Xm is given by (^ — x)P+, and the pion, 
with mass m^, is the lightest hadron, none of the particles in Xm can have a rapidity 
above T/max = ln((^ — a;)P+/m^). Thus, a rapidity gap of size Ay = ln(m^/(,^ — x)mp) 
exists between the outgoing proton and the state Xm- For typical values of ,^ ~ 10^^ (cf. 
the more detailed discussion of the experimental setting below) the size of this gap can 
be considerable. 

Note, however, that the term 'rapidity gap events' was coined to describe the appear- 
ance of diffractive events in the HERA frame, i.e., a frame defined by the electron-proton 
collision axis. The rapidity in this frame is, in general, different from the photon-proton 
frame rapidity discussed above. Nevertheless, the existence of a gap surrounding the 
outgoing proton in the 7*p frame clearly implies the existence of a similar gap in the 
ep frame. The exact size of the ep-frame rapidity gap follows from the specific event 
kinematics. The main conclusion so far is the kinematic separation of outgoing proton 
and diffractive state Xm in diffractive events with small ^. 



The appearance of a typical event in the ZEUS detector is shown in Fig. |2.2| , where 
FCAL, BCAL, and RCAL are the forward, central (barrel) and rear calorimeters. The 
absence of a significant energy deposit in the forward region is the most striking feature 
of this DIS event. In the naive parton model of DIS, a large forward energy deposit is 
expected due to the fragmentation of the proton remnant, which is left after a quark has 
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been knocked out by the virtual photon. In addition, the whole rapidity range between 
proton remnant and current jet is expected to fill with the hadronization products of 
the colour string that develops because of the colour charge carried by the struck quark. 
Thus, the event in Fig. p.2| shows a clear deviation from typical DIS, even though the 
outgoing proton left through the beam pipe and remains undetected. Hence, the name 
rapidity gap events is also common for what was called diffractive electroproduction in 
the previous paragraphs. 
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Figure |2.2| : Diffractive event in the ZEUS detector (figure from 



The definition of diffraction used above is narrower than necessary. Without losing 
any of the qualitative results, the requirement of a final state proton P' can be replaced 
by the requirement of a low-mass hadronic state Y, well separated from the diffractive 
state Xm- In this case, the argument connecting elastically scattered proton and rapidity 
gap has to be reversed: the existence of a gap between Xm and Y becomes the distinctive 
feature of diffraction and, under certain kinematic conditions, the interpretation of Y as 
an excitation of the incoming proton, which is now almost elastically scattered, follows. 

However, this wider definition of diffraction has the disadvantage of introducing a 
further degree of freedom, namely, the mass of the proton excitation. If one insists on 
using only the previous three parameters, x, and ^, the definition of a diffractive 
event becomes ambiguous. More details are found in the discussion of the experimental 
results below and in the relevant experimental papers. 
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2.2 Fundamental observations 



Rapidity gaps are expected even if in all DIS events a quark is knocked out of the 
proton leaving a coloured remnant. The reason for this is the statistical distribution of 
the produced hadrons, which results in a small yet finite probability for final states with 
little activity in any specified detector region. However, the observations described below 
are clearly inconsistent with this explanation of rapidity gap events. 

The first analysis of rapidity gap events at HERA was performed by the ZEUS 
collaboration 0. More than 5% of DIS events were found to possess a rapidity gap. 
The large excess of the event numbers compared to naive parton model expectations was 
soon confirmed by an HI measurement . The analyses are based on the pseudo-rapidity 
rj = — lntan(^/2), where 6 is the angle of an outgoing particle relative to the beam axis. 
Pseudo-rapidity and rapidity are identical for massless particles; the difference between 
these two quantities is immaterial for the qualitative discussion below. 

In the ZEUS analysis, a rapidity r/max was defined as the maximum rapidity of a 
calorimeter cluster in an event. A cluster was defined as an isolated set of adjacent cells 
with summed energy higher than 400 MeV. The measured r/max distribution is shown in 
Fig. p73| . (Note that the smallest detector angle corresponds to ?7max = 4.3; larger values 
are an artifact of the clustering algorithm.) 
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Figure |2.3| : Distribution of ?7max, the maximum rapidity of a calorimeter cluster in an 
event, measured at HERA (figure from 0]). 

To appreciate the striking qualitative signal of diffraction at HERA, the measured 
^max distribution has to be compared with naive expectations based on a purely partonic 
picture of the proton. This is best done using a parton-model-based Monte Carlo event 
generator. The corresponding r/max distribution, which is also shown in Fig. p.3| , is strongly 
suppressed at small //max- This qualitative behaviour is expected since the Monte Carlo 
(for more details see and refs. therein) starts from a partonic proton, calculates the 
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hard process and the perturbative evolution of the QCD cascade, and finally models the 
hadronization using the Lund string model (see, e.g., ||6^). According to the Lund model, 
the colour string, which connects all final state partons and the coloured proton remnant, 
breaks up via qq pair creation, thus producing the observed mesons. The rapidities of 
these particles follow a Poisson distribution, resulting in an exponential suppression of 
large gaps. 

It should be clear from the above discussion that this result is rather general and does 
not depend on the details of the Monte Carlo. QCD radiation tends to fill the rapidity 
range between the initially struck quark and the coloured proton remnant with partons. 
A colour string connecting these partons is formed, and it is highly unlikely that a large 
gap emerges in the final state after the break-up of this string. 

However, the data shows a very different behaviour. The expected exponential de- 
crease of the event number with ?7jnax is observed only above r/max — 1-5; below this value 
a large plateau is seen. Thus, the naive partonic description of DIS misses an essential 
qualitative feature of the data, namely, the existence of non-suppressed large rapidity 
gap events. 

To give a more specific discussion of the diffractive event sample, it is necessary to 
define which events are to be called diffractive or rapidity gap events. It is clear from 



Fig. |2.3| that, on a qualitative level, this can be achieved by an r^max cut separating the 
events of the plateau. The resulting qualitative features, observed both by the ZEUS |^ 
and HI collaborations [Q, are the following. 

There exists a large rapidity interval where the r/max distribution is flat. For high 7*p 
energies W, the ratio of diffractive events to all DIS events is approximately independent 
of W. The dependence of this ratio is also weak, suggesting a leading-twist contribu- 
tion of diffraction to DIS. Furthermore, the diffractive mass spectrum is consistent with 
a distribution. 

A number of additional remarks are in order. Note flrst that the observation of a 
flat r^max distribution and of a spectrum are interdependent as long as masses 

and transverse momenta of final state particles are much smaller than M^. To see this, 
observe that the plus component of the most forward particle momentum and the minus 
component of the most backward particle momentum are largely responsible for the total 
invariant mass of the diffractive final state. This gives rise to the relation dM^/M^ = 
(ilnM^ ~ c??7max, which is equivalent to the desired result. 

Furthermore, it has already been noted in that a significant contribution from 
exclusive vector meson production, e.g., the process YP ^ pp, is present in the rapidity 
gap event sample. A more detailed discussion of corresponding cross sections, which have 
by now been measured, and of relevant theoretical considerations is given in Chapters ^ 
and 1^. The discussion of other, more specific features of the diffractive final state, such 
as the presence of charmed mesons or high-p^ jets, is also postponed. 



19 



2.3 Diffract ive structure function 



The diffractive structure function, introduced in ||68| and first measured by the HI cof 



laboration is a powerful concept for the analysis of data on diffractive DIS, which 
is now widely used by experimentalists and theoreticians. 



Recall the relevant formulae for inclusive DIS (see, e.g., |7^). The cross section for 



the process ep eX can be calculated if the hadronic tensor, 

W,,{P, q) = ^J2< P\jlmX X X\j,mP > (27r) V(g + P - Px) , (2.3) 
47r ^ 

is known. Here j is the electromagnetic current, and the sum is over all hadronic final 
states X. Because of current conservation, q-W = W-q = 0, the tensor can be decomposed 
according to 

W,.{P, q) = (^g,. - W,{x, Q') + [P, + ^q,^ [p, + ^g,) W,{x, Q') . (2.4) 

The data is conveniently analysed in terms of the two structure functions 

F2(x,Q2) = {P-q)W2{x,Q^) (2.5) 
Fl{x,Q^) = {P-q)W2{x,Q^)-2xWi{x,Q^). (2.6) 



Introducing the ratio R = Fl/{F2 — Fl), the electron-proton cross section can be written 

as 

where y = Q'^/sx, and s is the electron-proton centre-of-mass energy squared. In the 
naive parton model or at leading order in as in QCD, the longitudinal structure function 
Fl{x,Q'^) vanishes, and R = 0. Since R corresponds to the ratio of longitudinal and 
transverse virtual photon cross sections, ai/cTT, it is always positive, and the corrections 
associated with a non-zero R are small at low values of y. 

In the simplest definition of diffraction, the inclusive final state X is replaced by the 
state Xm P', which consists of a diffractively produced hadronic state with mass M and 
the scattered proton. This introduces the two additional kinematic variables ^ and t. 
However, no additional independent 4-vector is introduced as long as the measurement 
is inclusive with respect to the azimuthal angle of the scattered proton. Therefore, the 
decomposition in Eq. (|2.4|) remains valid, and the two diffractive structure functions 
F^^'^\x,Q'^,C,,t) can be defined. The diffractive cross section reads 



dxdQ^didt xQ^ \ ^ 2[l + i?^(4)(x,Q2,e,t)]J ^"^'^'^''^^ ^^-^^ 

where R^ = F^ / {F2 —Fj^). In view of the limited precision of the data, the dominance of 
the small-y region, and the theoretical expectation of the smallness of F^ , the corrections 
associated with a non-zero value of R^ are neglected in the following. 
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Figure |2.4|: First measurement of tlie diffractive structure function F^^^^^ [x-^Q^ The 



fit is based on a factorizable ^ dependence of tlie form ^ ^'^^ (figure from |]69|). 

A more inclusive and experimentally more easily accessible quantity can be defined 
by performing the t integration, 

F^'^^Xx^Q^i) = I dtF^^'\x,Q\i,t). (2.9) 

The results of the first measurement of this structure function, performed by the HI 
collaboration, are shown in Fig. |2.4 (The underlying cross section includes events with 
small-mass excitations of the proton in the final state.) Far more precise measurements, 
since performed by both the HI and ZEUS collaborations, are discussed in Chapter ^ 
together with different theoretical predictions. 

The main qualitative features of diffractive electroproduction, already discussed in 
the previous section, become particularly apparent if the functional form of Fj^''^'* is 
considered. The (3 and dependence of Fj^*-^'* is relatively flat. This corresponds to the 
observations discussed earlier that diffraction is a leading twist effect and that the mass 
distribution is consistent with a 1/M|- spectrum. The success of a 1/^" fit, with n a 
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number close to 1, reflects the approximate energy independence of the diffractive cross 
section. More specifically, however, and in view of more precise recent measurements, it 
can be stated that the fitted exponent is above 1, so that a slight energy growth of the 
diffractive cross section is observed. This will be discussed in more detail later on. 



Note finally that the formal definition of ^2^*^^"* as an integral of F^^'^' is not 
easy to implement since the outgoing proton or proton excitation is usually not tagged. 
Therefore, most measurements rely on different kinematic cuts, in particular an rj^ax cut, 
and on models of the non- diffractive DIS background. A somewhat different definition 
of F^^^\ based on the subtraction of 'conventional DIS' in the distribution, was 



introduced in the ZEUS analysis of |71| . From a theoretical perspective, the direct mea- 
surement of F^^^^ by tagging the outgoing proton appears most desirable. Recently, such 
a measurement has been presented by the ZEUS collaboration |72| although the statistics 



are, at present, far worse than in the best available direct analyses of ^2^*^^"* [ [75| , [7^ 
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3 Semiclassical Approach 



In this chapter, the main physical idea and the technical methods of the semiclassical 
approach to diffraction are introduced. Although historically the semiclassical approach 
is not the first relevant model, it is, in the author's opinion, well suited as a starting point 
for the present review. On the one hand, it is sufficiently simple to be explained on a 
technical level within the limited space available. On the other hand, it allows for a clear 
demonstration of the interplay between the relevant kinematics and the fundamental 
QCD degrees of freedom. 

The underlying idea is very simple. From the proton rest frame point of view, the 
very energetic virtual photon develops a partonic fluctuation long before the target. 
The interaction with the target is modelled as the scattering off a superposition of soft 
target colour fields, which, in the high-energy limit, can be calculated in the eikonal 
approximation. Diffraction occurs if this partonic system is quasi-elastically scattered off 
the proton. This means, in particular, that both the target and the partonic fluctuation 
remain in a colour singlet state. 

3.1 Eikonal formulae for high-energy scattering 

The amplitude for an energetic parton to scatter off a given colour field configuration 
is a fundamental building block in the semiclassical approach. This amplitude is the 
subject of the present section. The two other basic ingredients for the diffractive cross 
section, i.e., the amplitudes for the photon to fluctuate into different partonic states and 
the integration procedure over all colour field configurations of the target proton are 
discussed in the remainder of this chapter. 

The essential assumptions are the softness and localization of the colour field and the 
very large energy of the scattered parton. Localization means that, in a suitable gauge, 
the colour field potential A^{x) vanishes outside a region of size ~ 1/A, where A is a 
typical hadronic scale. Later on, it will be assumed that typical colour field configurations 
of the proton fulfil this condition. Softness means that the Fourier decomposition of Afj_{x) 
is dominated by frequencies much smaller than the energy of the scattered parton. The 
assumption that this holds for all fields contributing to the proton state is a non-trivial 
one. It will be discussed in more detail at the end of this chapter. 

The relevant physical situation is depicted in Fig. |3.1| , where the blob symbolizes 
the target colour field configuration. Consider first the case of a scalar quark that is 
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minimally coupled to the gauge field via the Lagrangian 



C,c.i.r = (D^^y {D'^^) (3.1) 

with the covariant derivative 

D^ = d^ + igA^ . (3.2) 

In the high-energy limit, where the plus component of the quark momentum becomes 
large, the amplitude of Fig. ^]T] then reads 

i2'K5{k'^ - ko)T = 2TT6{k'Q - ko)2ko \u{k'^ - k^) - {27ry6\k'^ - k^)] . (3.3) 



It is normalized as is conventional for scattering processes off a fixed target (see, e.g.. 
The expression in square brackets is the Fourier transform of the impact parameter space 
amplitude, U{x±) — 1, where 

U{x±) = P exp (^—-^ J A_{x+,x±)dx^ (3.4) 

is the non-Abelian eikonal factor. The unit matrix 1 G SU{Nc), with Nc the number of 
colours, subtracts the field independent part, and the path ordering operator P sets the 
field at smallest x+ to the rightmost position. 




Figure |3.1| : Scattering of a quark off the target colour field. 



This formula or, more precisely, its analogue in the more realistic case of a spinor 
quark was derived by many authors. In the Abelian case, the high-energy amplitude was 
calculated in ^ in the framework of light-cone quantization. This result was taken over 
to QCD in [Q. A derivation in covariant gauge, based on the solution of the equation 



of motion for a particle in the colour background field, was given in |^3|. In the 



amplitude for the scattering of a fast gluon off a soft colour field was derived by similar 
methods. 



For completeness, a derivation of the amplitude in Eq. (^.31), based on the summation 
of diagrams of the type shown in Fig. |3.2| , is given in Appendix ^ of the present review. 
For the purpose of this section, it is sufficient to explain the main elements of Eq. ( ^.3] ) 
in a physical way, without giving the technical details of the derivation. 

To begin with, it is intuitively clear that the eikonal factor U{x±) appears since 
the fast parton, travelling in x+ direction and passing through the target at transverse 
position x±, is rotated in colour space by the field Af^{x) that it encounters on its way. 

Furthermore, the above amplitude is given for the situation in which the parton is 
localized at x_ ^ 0. This co-ordinate of the fast parton does not change during the 
scattering process. Thus, is always evaluated at a;_ ~ 0, and the a;_ dependence is 
not shown explicitly. 
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Figure |3.2| : Typical diagrammatic contribution to the eikonal amplitude, Eq. ( p.3| ). 
Attachments of gluon lines with crosses correspond to vertices at which the classical 
external field appears. 



The energy 5-function in Eq. (|3.3|) is an approximate one. It appears because the 
energy of the parton can not be significantly changed by the soft colour field. 

Finally, due to the explicit factor ko, the amplitude grows linearly with the par- 
ton energy, as is expected for a high-energy process with t-channel exchange of vector 
particles, in this case, of gluons. 



The amplitude of Eq. ( |3.3D is easily generalized to the case of a spinor quark, where 
the new spin degrees of freedom are characterized by the indices s and s' (conventions 
of Appendix ^). In the eikonal approximation, which is valid in the high-energy limit, 
helicity fiip contributions are suppressed by a power of the quark energy (see Appendix ^ 
for more details). Thus, the generalization of Eq. (|3.3|) reads 



i27r6{ko - ko)Tss' = 27r5(fc[, - ko)2ko \u{k'^_ - k^) - (27r)2(5^(fcl - k 



5. 



(3.5) 



Similarly, the amplitude for the scattering of a very energetic gluon off a soft colour 
field is readily obtained from the basic formula, Eq. ( |3.3D . Note that, although the fast 
gluon and the gluons of the target colour field are the same fundamental degrees of 
freedom of QCD, the semiclassical approximation is still meaningful since an energy cut 
can be used to define the two different types of fields. The polarization of the fast gluon is 
conserved in the scattering process. The main difference to the quark case arises from the 
adjoint representation of the gluon, which determines the representation of the eikonal 



factor. Thus, the amplitude corresponding to Fig. |3.3| reads 



i27c6{k'o - ko)Txy = 2T^5{k'^ - ko)2ko \lf-^{k'^ - k^) - {27Ty6\k'^ - k^) 



6 



XX' 



(3.6) 



where A, A' are the polarization indices and U-^ is the Fourier transform of U-^{x±), the 
adjoint representation of the matrix U{x±) defined in Eq. (|3.4|). 




Figure 3.3: Scattering of a gluon off the target colour field. 



Clearly, the above eikonal amplitudes have no physical meaning on their own since 
free incoming quarks or gluons cannot be realized. However, they serve as the basic 
building blocks for the high-energy scattering of colour neutral objects discussed in the 
next section. 
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3.2 Production of qq pairs 

In this section, the eikonal approximation is used for the calculation of the amplitude for 



qq pair production off a given target colour field . Both diffractive and inclusive cross 
sections are obtained from the same calculation, diffraction being defined by the require- 
ment of colour neutrality of the produced pair. The qualitative results of this section 
are unaffected by the procedure of integrating over all proton colour field configurations, 
which is discussed in Sect. 




Figure |3.4| : Electroproduction of a qq pair off the target colour field. 



The process is illustrated in Fig. |3.4 The corresponding T matrix element has three 
contributions, 

T = Tgg+Tg + Tg, (3.7) 

where Tgg corresponds to both the quark and antiquark interacting with the field, while 
Tq and Tg correspond to only one of the partons interacting with the field. 

Let Vq{p',p) and Vq{k',k) be the effective vertices for an energetic quark and anti- 
quark interacting with a soft gluonic field. For quarks with charge e, the first contribution 
to the amplitude of Eq. (p^) reads 



i2TT5{kQ + ko- qo)Tqq = ie [ j!r^Us'ip')Vq{p ,p)—^ — /(?)^r Vq{k, k')vr'ik') , 

J [ZTTj* p — m —f — m 

(3.8) 

where q = p + k hj momentum conservation, e(g) is the polarization vector of the 
incoming photon, and r', s' label the spins of the outgoing quarks. 

At small X, the quark and antiquark have large momenta in the proton rest frame. 
Hence, the propagators in Eq. ( p.8|) can be treated in a high-energy approximation. In a 
co-ordinate system where the photon momentum is directed along the z-axis, the large 
components are p+ and k^. It is convenient to introduce, for each vector k, a vector k 
whose minus component satisfies the mass shell condition, 

k. = ^4^, (3.9) 



while the other components are identical to those of k. The propagators in Eq. ( |3.8|) can 
be rewritten according to the identities 

1 _ ErVr{k)Vr{k) 



It + m fc^ — 2k 

1 EsUs{p)us{p) , 7+ 



pl — m — 2p_( 



(3.11) 
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To obtain the first term in a high-energy expansion of the scattering amphtude Tgg, the 
terms proportional to 7+ in Eqs. (|3.1CI| ) and ( |3.11|) can be dropped. 

After inserting Eqs. ( |3.1CI| ) and into Eq. (p.8|), the relations 



Us'ip')Vq{p',p)usip) 

Vr{k)Vq{k,k')Vr'{k') 



2n5ip', - po)2po - p^) - i27rY5'ip'^ - p^) 

2TT6{k'o - ko)2ko \u\k±-k'^) - {2nf6\k'^ - k^ 



5r 



(3.12) 
(3.13) 



which correspond to Eq. ( |3.5| ) and its antiquark analogue, can be applied. Writing the 
loop integration as d'^k = {l/2)dk+dk-(fk± and using the approximation 6{lo) ~ 26{l+) 
for the energy 5-functions, the fc+ integration becomes trivial. The k_ integral is done 
by closing the integration contour in the upper or lower half of the complex /c_ plane. 
The result reads 



le 



47r^ 



94 



d'^k^ 



N^ + k^, 



(3.14) 



Uip'± - Pi) - {2tx)H\p'^ - p^)] p{k^ - k'^) - {2ixY5\k'^ - k^) 



where 



P+ = (1 - a) q+ , fc+ 



a q+ 



(3.15) 



Thus, a and 1—a characterize the fractions of the photon momentum carried by the two 
quarks, while (iV^ + k\) measures the off-shellness of the partonic fluctuation before it 
hits the target. In the following, the quark mass is set to m = 0. 

The above expression for Tqq contains terms proportional to UU\ U, and t/^, as 
well as a constant term. The amplitudes Tg and Tq, which contain terms proportional to 
U and [/^ and a constant term, are derived by the same methods. Calculating the full 
amplitude according to Eq. ( p.7|) , the terms proportional to U and W cancel. Thus, the 
colour field dependence of T is given by the expression 



U{p'i. - P±)UHk± - fcl) - {27r)'6\p'^ - PA_)5\k'^ - k^) 
Introducing the fundamental function 



(3.16) 



(3.17) 



which encodes all the information about the external field, the complete amplitude can 
eventually be given in the form 



T 



le 



47r2 



q+ i d^kA 



ail 



a] 



Us'{j))4{q)Vr'{k) 



-iA I X 



W^Ak'^-ki 



(3.18) 

k'^+p'^ 



where Wx^ is the Fourier transform of Wx^iyi.) with respect to j/^, and Aj 
is the total transverse momentum of the final qq state. 

From the above amplitude, the transverse and longitudinal virtual photon cross 
sections are calculated in a straightforward manner using the explicit formulae for 
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Us'{p)4{q)Vr-i{k) given in Appendix Summing over all qq colour combinations, as ap- 
propriate for the inclusive DIS cross section, the following result is obtained. 



da 



L 



da dk'l 

dax 
da dk'? 



(a^ + n-af) 



2u W^Ak'^-k^] 



d'k 



N^ + k^, 



2(27r) 



d'^k^ 



k±W^Akl-k, 



m + k^. 



(3.19) 



(3.20) 



Note that only a single integration over transverse coordinates appears. This is a con- 
sequence of the 5-function induced by the phase space integration over A^, applied to 
exp[— -iAxXx] from Eq. (p.l8|) and to the corresponding exponential from the complex 
conjugate amplitude. The contraction of the colour indices of the two W matrices is 
implicit. 

Consider the longitudinal cross section in more detail. The integrand can be expanded 
around k± = k'^. Shifting the integration variable k± to l± = k^ — k'^, the Taylor 
expansion of the denominator in powers of l± yields 



1 



1 



21 1 k\ 



m + {k'^ + /x)2 + k'l (iV2 + k'lY 



+ 



(3.21) 



From the definition of the colour matrix W^AV-^^ -^l- ( P.lTp , it is clear that 



dH^w^A-^A.) = (277)2 ly, JO) ^ Q 



dH^ W.J-/±) 



(3.22) 
(3.23) 



Using rotational invariance, i.e., k'M Id^k''^, , the result 



dar 



Ae^a{l-a) N^kf 



da dk'l {2TiY{m + k'lY 
is obtained. It evaluates to the total longitudinal cross section 



(3.24) 



I \di_W^M\^ . 

J X I 



(3.25) 



The transverse contribution can be evaluated in a similar way. In the perturbative 
region, where a(l — a) 3> A^/Q^ g^^d k'l ^ A^, the integrand can again be expanded 
around k± = k'j_. The analogue of Eq. ( |3.24| ) reads 



dai 



da dk'l 



e^[a^ + {l-aY] [N^ + k 
2{2T^f{N^ + k'lY 



Ml 



(3.26) 



While, in the longitudinal case, the a and k'^ integrations were readily performed, a 
divergence is encountered in the transverse case, Eq. ( [3.261 ). The k'^ integration gives 
rise to a factor l/a(l — a), so that the a integral diverges logarithmically at a — and 



28 



at a — >■ 1. Thus, the region of small a{l — a), where the expansion around fc^ = fc^ 
does not work, is important for the total cross section 0"^, which therefore can not be 
obtained by integrating Eq. ( |3.26| ). Instead, the endpoint region, where the large distance 
structure of Wxj_{yj_) is important, can be separated by a cutoff /i^ (A^ ^ ^ Q'^). 
The complete leading twist result for the transverse cross section, where contributions 
suppressed by powers of A^//i^ or jj? /Q"^ have been dropped, reads 



67r2g2 



In 



- 1 




(3.27) 



, da I dk'f 
(27r)6 7o J ^Jx^ 



The resulting physical picture can be summarized as follows. For longitudinal pho- 
ton polarization, the produced qq pair has small transverse size and shares the photon 
momentum approximately equally. Only the small distance structure of the target colour 
field, characterized by the quantity |9_lVF2^^(0)P, is tested. For transverse photon polar- 
ization, an additional leading twist contribution comes from the region where a or 1—a 
is small and k'1 ~ A^. In this region, the qq pair penetrating the target has large trans- 
verse size, and the large distance structure of the target colour field, characterized by 
the function Wx^ijji) at large y^, is tested. This physical picture, known as the aligned 
jet model, was introduced in on a qualitative level and was used more recently for a 
quantitative discussion of small-x DIS in 0]. 

It is now straightforward to derive the cross sections for the production of colour 
singlet qq pairs corresponding, within the present approach, to diffractive processes. Note 
that the cross sections in Eqs. (|3.19|) and ( p.20|) can be interpreted as linear functionals of 
tT{Wx^{ili.)Wl^{il'^), where the trace appears because of the summation over all colours 
of the produced qq pair in the final state. Introducing a colour singlet projector into the 
underlying amplitude corresponds to the substitution 



_(yx)trW^t (yl 



(3.28) 



in Eqs. ( |3.19| ) and ( |3.2CI| ). This change of the colour structure has crucial consequences 
for the subsequent calculations. 



Firstly, the longitudinal cross section, given by Eqs. ( p.24| ) and ( p.25| ), vanishes at 
leading twist since the derivative d±Wx^ (0) is in the Lie-algebra of SU (N^), and therefore 
trS^W^.jO) =0. 

Secondly, for the same reason the InQ^ term in the transverse cross section, given by 
Eq. ( |3.27| ), disappears. The whole cross section is dominated by the endpoints of the a 
integration, i.e., the aligned jet region, and therefore determined by the large distance 
structure of the target colour field. At leading order in l/Q"^, the diffractive cross sections 
read 







i2TiYN^ Jo 



da / dk 



/2 



2, k^tiW^^ik'^ - ki_) 



d'k 



m + k^. 



(3.29) 
(3.30) 
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The cutoff of the a integration, yU^/Q^, has been dropped since, due to the colour singlet 
projection, the integration is automatically dominated by the soft endpoint. 

In summary, the leading-twist cross section for small-x DIS receives contributions 
from both small- and large-size qq pairs, the latter corresponding to aligned jet configu- 
rations. The requirement of colour neutrality in the final state suppresses the small-size 
contributions. Thus, leading twist diffraction is dominated by the production of pairs 
with large transverse size testing the non-perturbative large-distance structure of the 
target colour field. 



3.3 Higher Fock states 



Given the leading order results for qq pair production of the last section, it is natural to 
ask about the importance of radiative corrections. A systematic procedure for calculating 
to all orders in perturbation theory does not yet exist in the semiclassical framework. 
However, as will be seen in the following chapters, the summation of leading logarithms 
in Q"^ is understood. Here, the particularly important case of the diffractive production 

PI. The 



of a quark-antiquark-gluon system (see Fig. p.5|) is discussed in some detail 
purpose of this discussion is to establish, as one of the essential features of diffractive 
DIS, the necessary presence of a wee parton in the wave function of the incoming virtual 
photon. 



a) 



b) 




Figure |3.5| : Diagrams for the process 7* — qqg. Two similar diagrams with the gluon 
radiated from the produced antiquark have to be added. 



Consider the sum of diagram Fig. |3.5| a and its analogue with the gluon radiated from 
the antiquark. If the quark propagators with momenta p and I are rewritten according 
to Eqs. ( 3.101 ) and (|3.11 ) and the 7+ terms are dropped, the corresponding amplitude 
can be given in the form 

i27r%'o + + /[, - go)T(a) = /^/ E^i'i.rV.A'A^iS • (3.31) 

Here T^^^ is the amplitude for the scattering of the qqg system off the target colour field. 



30 



and T^^) is the remainder of the diagram, describing the fluctuation of the virtual photon 
into the partonic state. In the following calculation, the k_ integration will be performed 
in such a way that the gluon propagator goes on shell. Anticipating this procedure, 
the sum over the intermediate gluon polarizations A is restricted to the two physical 
polarizations, defined with respect to the on-shell vector k. According to Eqs. ( |3.5| ) and 
( p.6|) , the contribution to T^^^ where quark, antiquark and gluon interact with the field 
reads 



V 119 J s's,r'r,X'> 



i27i5{p'^ - po)2po [Uip'^ - p^) - {2'nf5\p'^ - p^ 
xi27rS{l'o - lo)2lo [f/^(/± - /I) - {27if5\l'^ - l±) 
xi27c6{k'Q - ko)2ko \u^{k'^ - ki_) - {2Tif5^{k'^ 



6r'r 



ki 



X'X 



(3.32) 
(3.33) 
(3.34) 



As explained at the beginning of Sect. p.2| , contributions where not all of the partons 
interact with the field have to be added. 

The photon-qqg transition amplitude T^'^\ with all colour indices suppressed, reads 



T. 



(2) 



srX 



legs 
p'^k'^P 



Usip) 



m 



-m+m- 



m 



Vr{l) . (3.35) 



The colour structure of T(a) is given by the combination of the U matrices in T^^^ , the 
indices of which are partially contracted according to the perturbative amplitude T^'^\ 
This colour structure is characterized by the colour tensor 



(iy(3))^ 



Kw^-'\x^,y^,z^mf,, (3.36) 



where 



rpa 



and T^p are the conventional SU{Nc) generators with adjoint 



(3.37) 
A^^-i) and 



fundamental (a, /3 = 1 ■ ■ ■ N^) indices. The notation W^^^ is chosen to stress the similarity 
with the qq case, where the external field is tested by the function W = W^'^\ defined 
in Eq. ( |3.17[ ). In analogy to this equation, the last term in Eq. ( |3.37D subtracts the 
unphysical contribution where none of the partons is scattered by the external field. One 
can think of x±, y± and z± as the transverse positions at which quark, gluon and antiquark 
penetrate the proton field, picking up corresponding non-Abelian eikonal factors. The 
indices a, (3 and a correspond to the colours of the produced quark, antiquark and gluon. 

To obtain the complete result, the amplitude T(b), which is the sum of the diagram in 



Fig. p75| b and its analogue with the gluon radiated from the antiquark, has to be added. 
The calculation is similar to the case of the amphtude T(a). Since the gluon is radiated 
after the quark pair passes the target field, the colour structure is determined by the 
same function W that appeared in the qq production amplitude of the last section. 

The diagrams in Fig. |3.5| , with summation over all possible colours of the qqg final 
state, represent an correction to inclusive qq pair production according to Fig. The 



factor as is accompanied by a factor \nQ . Higher Fock states lead to additional factors 
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tts InQ^, and an all-orders summation can be performed using standard renormalization 
group techniques, i.e., Altarelli-Parisi evolution. The implementation of this summation 
in the semi classical framework is explained in more detail at the end of Sect. [4.4| and in 
Appendix 0. Hence, the inclusive case is not further discussed in this section. 

For the diffractive case, the following situation can be anticipated on the basis of the 
experience of the last section. There are three qualitatively different kinematic configu- 
rations. First, if all three transverse momenta, p'^, k'^ and l'^ are soft, the large distance 
structure of the target field is tested, and a leading twist contribution to the DIS cross 
section results. Second, if all three momenta are large, p']_^ A;^, ^ A^, all three po- 
sitions x^, and z±_ in Eq. ( |3.37| ) are close together, and an expansion in powers of 
the target colour field can be performed. The leading term corresponds to colour octet 
exchange, so that no leading twist contribution to diffractive DIS results. Third, if two 
of the produced partons have high transverse momenta and the remaining parton is 
soft, a leading twist diffractive contribution arises. The two possible situations, depicted 
in Figs. |3]^ and |3?7| , are discussed quantitatively below. The results of this calculation 
justify the physical picture outlined above. 




Figure |3.6| : Space-time picture in the case of fast, high-p^ quark and antiquark, passing 
the proton at small transverse separation, with a relatively soft gluon further away. 




Figure |3.7| : Space-time picture in the case of fast, high-p^ quark and gluon, passing the 



proton at small transverse separation, with a relatively soft antiquark further away. 

Consider first the case of high-p_L quark and antiquark, i.e., p^, I'l ^ A^, and rel- 
atively soft gluon, i.e., k'1 ~ A^ and a' ^ 1. Here a' = k'^/qo is the analogue of the 
variable a of the last section. 

The assumption of a smooth external field implies small transverse momentum trans- 
fer from the proton, i.e., |/" | ~ A, where l'^ = — l±. The l± integration in Eq. ( |3.31| ) 
can be replaced by an integration, substituting at the same time 

l± = l'±- I'L and = + l'[-k±. (3.38) 

Neglecting in T^'^\ which is justified since <^ \p±\, l^'ll, the only remaining 
dependence is located in the colour factor W^^\ This simplifies the integration to 

(3,39) 
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Defining A = p' + k' + I' — p — k — I to he the total momentum transferred from the 
proton, ly(^) can be given in the form 



where l± is given by Eq. (|3.38|) . The integration gives a 5-function of the variable z±, 



thus resulting in the final formula 



(2vr)^ 



iW^'%= / e-^^^^W e'y^^'^~'''^^W^''\x^,x^ + y^,x^r,p. (3.41) 



This expression shows that in the kinematic situation with two high-p_|_ quark jets and 
a relatively soft gluon the leading twist contribution is not affected by the transverse 
separation of the quarks. It is the transverse separation between quark-pair and gluon 
which tests large distances in the proton field and which can lead to non-perturbative 
effects. 

The colour singlet projection of the colour tensor W^^^ reads 



S{W^'^) = ^jj;^ (W^'^T^,^. (3.42) 

Using the identity 

(^^Ayb ^ 2 tip-^Tm''] , (3.43) 

where U G SU{Nc), the contribution relevant for diffraction, i.e., the production of a 
colour singlet gggf-system, takes the form 

/ SiW^'^) = I e--^- , ^ tAW^, ik± - k'^)] , (3.44) 
J (27r)^ Jx^ ^2{N^-l) 

where 

W^^{y^) = U\x^)U^\x^ + y^)) - 1 . (3.45) 

This is analogous to the quark pair production of the previous section (cf. Eq. (|3T7|) ). 
However, now the two lines probing the field at positions x±_ and x±_ + y\_ correspond to 
matrices in the adjoint representation. An intuitive explanation of this result is that the 
two high-px quarks are close together and are rotated in colour space like a vector in the 
octet representation (cf. Fig. [3l^ ). 

To make this last statement more precise, recall that an upper bound for the loffe- 
time of the fluctuation with two high-p_|_ quarks is given by qQ/p\- This means that 
the distance between the point where the virtual photon splits into the gg-pair and 
the proton can not be larger than go/Pj_- As long as the pair shares the longitudinal 
momentum of the photon approximately equally, i.e., a{l—a) = 0{1), the opening angle 
is ~ Px/go- Hence, when quark and antiquark hit the proton, their transverse distance 
is ~ l/\p±\ < 1/A. 

The diagram in Fig. |3.5|b and its analogue with the gluon radiated from the antiquark 



do not contribute to the diffractive cross section for high-p^ QQ pair production. This 
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is obvious since the colour field is probed only by the qq system, the small transverse 
size of which prevents colour singlet exchange. Thus, the full cross section follows from 
Eq. (|3.31|) , the simplified colour tensor of Eq. ( p.44|) , and the formulae of Appendix 0, 
which are used for the evaluation of the qq'y and qqg vertices in Eq. (|3.35|) . 

The results for longitudinal and transverse photon polarization read 
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(3.46) 



/l(a'g^A;i), (3.47) 



a'g2 + k\ 



Ar^ = a(l-a)Q^ 



(3.48) 



(3.49) 



a(l — a) 

At this point, the fundamental assumption of the softness of the gluon can be quanti- 
tatively justified. To achieve this, fix a and p'^ and perform the integration over oi and 
k'^. As in Eq. ( ^.3(J ), these two integrations are dominated by the region a' <^ 1 and 
k'1^ ~ A^. Thus, the diffractive production of a qqg system with high-p^ quark and anti- 
quark occurs predominantly in a generalized aligned jet configuration, where the gluon 
is the wee parton in the photon wave function. It is this soft gluon that is responsible 
for the resulting leading twist cross section with colour singlet exchange. 

As already pointed out above, another leading twist contribution to diffraction arises 
from the kinematic region where the quark or the antiquark is the soft particle and, 
accordingly, the outgoing antiquark-gluon or quark-gluon system has large transverse 
momentum. The calculation proceeds along the same lines as in the soft gluon case. The 
main qualitative difference is in the colour structure. In analogy to the discussion leading 
to Eq. ( 3.41| ), the large p^ of quark and gluon results in a small transverse separation 
(cf. Fig. ISTTI) , so that effectively the high-p^ quark-gluon system is colour rotated like 
a single quark. Accordingly, the field is tested by the same function W ^ built from V 
matrices in the fundamental representation, that appeared in the qq production cross 
section. Notice also that, in contrast to the soft gluon case, the diagram in Fig. |3.5| b has 
to be taken into account. 



The final results for the transverse and longitudinal cross sections read 
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(3.50) 
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Ua'Q'M, (3.51) 



(3.52) 



34 



These cross sections include both the soft quark and the soft antiquark regions. The 
kinematic variables do not correspond to Fig. |3.5| . They are generic in the sense that p'j^ 
and — are the transverse momenta of the two hard jets, a and 1— a are the correspond- 
ing momentum fractions, and the soft parton, in this case the quark or the antiquark, is 
characterized by the transverse momentum k'^ and the longitudinal momentum fraction 
a'. These conventions are chosen to emphasize the similarity with the soft gluon result 
of Eqs. ( |3.46| )-( |3.49| ). Again, the a' and k'^ integrations are dominated by the region 
a' -Cl and k'} ~ 



3.4 Field averaging 

So far, electroproduction off a fixed 'soft' colour field has been considered. As a conse- 
quence, electroproduction cross sections approach constant values as a; ^ 0. However, 
a proper treatment of the target requires the integration over all relevant colour field 
configurations. 



Following the discussion of |^5|, consider first the elastic scattering of a quark off 



a proton. Although this process is unphysical since quarks are confined, it can serve to 
illustrate the method of calculation. Therefore in the following, confinement is ignored, 
and quarks are treated as asymptotic states. The generalization to the physical case of 
electroproduction is straightforward and will be discussed subsequently. 

A point-like quark with initial momentum q scatters off the proton, which is a rela- 
tivistic bound state with initial momentum p. Let be the proton mass, and s and t the 
usual Mandelstam variables for a 2 ^ 2 process. In the high-energy limit, s ^ t, m-p, 
the contribution from the annihilation of the incoming quark with an antiquark of the 
proton is negligible. The amplitude is dominated by diagrams with a fermion line going 
directly from the initial to the final quark state. Therefore, the proton can be described 
by a Schrodinger wave functional $p[A] (cf. |T^) depending on the gluon field only. 
Quarks are integrated out, yielding a modification of the gluonic action. 



For a scattering process, the amplitude can be written in the proton rest frame as 



<q'P'\qP>= lim / DAtDA^t^*p,[At\^p[A^t\ [^^ L'Ae^'^t^l <q'\q>A ■ (3.53) 

T^oo J J A-T 



Here the fields A_t and At are defined on three-dimensional surfaces at constant times 
—T and T, and A is defined in the four-dimensional region bounded by these surfaces. 
The field A has to coincide with A^t and At at the boundaries, and the action S is 
defined by an integration over the domain of A. The amplitude <q'\q>A describes the 
scattering of a quark by the given external field A. 

The initial state proton, having well defined momentum P, is not well localized 
in space. However, the dominant field configurations in the proton wave functional are 
localized on a scale A ~ Aqcd- The field configurations A{x, t), which interpolate between 
initial and final proton state, are also localized in space at each time t. Assume that the 
incoming quark wave packet is localized such that it passes the origin x = at time 
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t = 0. At this instant the field configuration A{x,t) is centered at 

x[A] = J d^xEA{x) -x J d^xEA{x) , (3.54) 



where Ea{x) is the energy density of the field A{x, t) aXt = 0. The amplitude ( |3.53[ ) can 
now be written as 

<q'P'\qP>= lim / £x [ DAtDA_t^*p,^p [^^ DAe'^6%x[A] - x) <q'\q>A ■ 

(3.55) 

Using the transformation properties under translations, 

<q'\q>L,A= e'^^-'i'^^ <q\q>A , ^*p.[L^A]^ p[LsA] = e'^''-^'^^^*p\A]^p[A] , 

(3.56) 

where 

L,A{y) = A{y - x) , (3.57) 

one obtains, 

<q'P'\qP>=2m^{2Tif5^{P' + q'-P-q)[ <q'\q>A. (3.58) 

J{A} 

Here /|^| denotes the operation of averaging over all field configurations contributing to 
the proton state which are localized at x = at time t = 0. It is defined by 



/ ^!\A] = lim / DAtDA_t^*p,^p f ^ DAe'^6^(x\Am\A] (3.59) 

J{A} 2mpT^ooJ ^ Ja_t 

for any functional The normalization /|^| 1 = 1 follows from 

<p'|p>= 2Po{2'Kf6%P' -P) . (3.60) 

More complicated processes can be treated in complete analogy as long as the proton 
scatters elastically. In particular, the above arguments apply to the creation of colour 
singlet quark- ant iquark pairs, 

<qqP'\-i''P>=2mp{2Tif5\kf-ki) f <qq\-f*>A, (3.61) 

J{A} 

where ki and kf are the sums of the momenta in the initial and final states respectively. 
The generalization of this simplest diffractive process to a process with an additional 
fast final state gluon, 7* qqg (cf. Sect. p73|) , is straightforward. In contrast to the 
quark-proton scattering discussed above, here a colour neutral state is scattered off the 
proton. Therefore no immediate contradiction with colour confinement arises. However, 
it has to be assumed that the hadronization of the produced partonic state takes place 
after the interaction with the proton, which is described in terms of fast moving partons. 

The amplitude for the scattering off a soft external field contains an approximate 
energy (5-function, giving rise to the definition of a functional F, 

< qq\Y >A= + fcO - gO) F[A] , (3.62) 
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where g, fcg, kg are the momenta of the incoming photon and the outgoing quark and 
antiquark respectively. Since the energy transferred to the proton is small, Eq. (|3.61|) can 
now be written as 

<qqP'\-f*P>=2mp{27Ty6\kf-ki) [ F[A] . (3.63) 

J{A} 



When calculating the cross section from Eq. (|3.63| ), the square of the 4-momentum- 



conserving 5-function translates into one 4-momentum-conserving 5-function using 
Fermi's trick. The spatial part of this (5-function disappears after the momentum in- 
tegration for the final state proton. The squared amplitude, integrated over the phase 
space of the outgoing proton and normalized to the total space-time volume, reads 
, 2 
— J d^^^^l < qqP' I 7*P > p = ATim^5{kl + k^ - g°) F[A] , (3.64) 

as it should for the scattering off a superposition of external fields. 

From the above discussion, a simple recipe for the calculation of diffractive processes 
at high energy follows: 

The partonic process is calculated in a given external colour field, localized at x = 
at time t = 0. The weighted average over all colour fields contributing to the proton 
state is taken on the amplitude level. It is assumed that the typical contributing field 
is smooth on a scale A and is localized in space on the same scale A. Finally, the cross 
section is calculated using standard formulae for the scattering off an external field. 

The above discussion was limited to the case of diffraction and did not address 
the important question of how inclusive and diffractive electroproduction processes are 
interrelated. To arrive at a combined treatment of both processes, additional assumptions 
concerning the treatment of the target colour field have to be introduced. A corresponding 
discussion, following the analysis of |^, is outlined below. 

As explained previously, in the semiclassical framework it is natural to expect diffrac- 
tion to occur whenever the produced qq pair emerges in a colour singlet state. Non- 
diffractive events are expected in the colour octet case. A combined treatment clearly 
forbids the use of the proton wave functional for the description of the final state. Thus, 
instead of the amplitude in Eq. ( |3.61 ), one has to consider the corresponding amplitude 



in a 'mixed' representation, < qqA \ 7*P >, where the final state consists of the outgoing 
qq pair and a colour field configuration A. This formalism allows for both the creation 
of colour singlet and non-singlet partonic states. 

For simplicity, time evolution of the field between the actual scattering process and 
the moment at which the final state field configuration A is defined is neglected. The 
squared amplitude, summed over all fields A and normalized to the total space-time 
volume, reads 

^ J DA\< qqA\ 7*P > 1^ = 47rmp5(A;;+A;°-g°) J DA,,, \ %[A,,,] F[A,,,] |' . (3.65) 

Here the integral over A on the l.h. side replaces the phase space integral for the outgoing 
proton in Eq. ( p.64| ). The index 'loc' symbolizes that, on the r.h. side of Eq. ( p.65|) . 



37 



the integration is restricted to fields localized at, say, a; = 0. This can be achieved 
using translation covariance of the proton wave functional and of the matrix element 
< gg I 7* >A, in an argument similar to the one leading to Eqs. ( p.59| ) and (|3.64|) . 

When writing Wx^^{yi), it has so far always been assumed that the functional depen- 
dence on the classical colour field configuration Ac\ is implicit, so that one should really 
read W4x(l/±)[^ci]- As can be seen from Eq. ( |3.65|) , the full inclusive cross sections are 



obtained from the previous formulae by the substitution 



(3.66) 

The same applies to the diffractive cross sections, which are obtained by introducing a 
colour singlet projector on the amplitude level (cf. Eq. ( |3.28| )). 

Decomposing the field A\oc in Eq. ( p.65| ) into its Fourier modes A\oc{k), the path 
integral can be written as 

n IdAUk), (3.67) 



^1°^= |fc|«k1 



where the cutoff \q\ is required to ensure that the basic precondition for the semiclassical 
treatment, the softness of the target colour field with respect to the momenta of the 
fast particles, is respected. This cutoff induces a non-trivial energy dependence of the 
squared amplitude in Eq. ( |3.65| ) and therefore of both the inclusive and diffractive cross 



sections. 

At present, no complete derivation of the explicit form of that energy dependence 
from first principles exists. However, a number of interesting related developments, ap- 
proaching the high-energy limit of QCD from the perspective of colour fields and eikon- 
alized interactions, can be found, e.g., in |]77|,[r5[]. 

On the basis of the above qualitative picture, a soft, non-perturbative energy growth 
was ascribed to the input parton distributions used in the phenomenological analysis 
of (see Sects. and ^]T] for more details). 

The discussion of the present section is closely related to Nachtmann's original pro- 



posal to treat high-energy hadron-hadron scattering in the eikonal approximation [43 
The above viewpoint corresponds to the rest frame of the proton, where the colour field 
encountered by the fast partons of the projectile is naturally considered to be part of 
the proton state. The viewpoint of |^ corresponds to the centre-of-mass frame of high- 



energy processes. In hadron-hadron scattering, both incoming particles are characterized 
by their parton content. These partons then interact via gluon fields which belong to 
neither of the two hadrons. So far, this is completely general. However, once the eikonal 
approximation is used to treat the parton propagators in the colour field background, 
one is forced to rely on modelling the light-cone wave function of the proton. In this case, 
the method becomes more predictive but less general than the wave functional treatment 
of the proton target discussed above. 
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4 From Soft Pomeron to Diffractive 
Parton Distributions 



In the previous chapter, small-x DIS was described as the eikonahzed interaction of a par- 
tonic fluctuation of the virtual photon with a superposition of colour field configurations 
of the proton. Diffraction occurs whenever the partonic fluctuation preserves its overall 
colour neutrality. An essential feature of this approach is its formulation exclusively in 
terms of the fundamental degrees of freedom of QCD. 

It is now appropriate to step back and take a look at diffraction from the historical 
perspective of soft hadronic physics (see [^] for a recent review of hadronic diffraction). 
It will soon become clear how partonic ideas arise in this framework and in which way 
they are related to the results of the last chapter. 



4.1 Soft pomeron 



Before introducing the pomeron, a concept from soft hadronic physics, it is appropriate 
to give a simple argument why, even at very high photon virtualities, diffractive DIS is 
largely a soft process. Such an argument was presented by Bjorken and Kogut in the 
framework of their aligned jet model 0] (see also [BO]). 




Figure |4.1| : Vector meson dominance inspired picture of inclusive (a) and diffractive (b) 
electropro duct ion . 

The underlying physical picture is based on vector meson dominance ideas. At high 
energy or small x, the incoming photon with virtuality = —Q^ fluctuates into a 
hadronic state with mass M, which then collides with the target (see Fig. [4.1| a). The 
corresponding cross section for transverse photon polarization is estimated by 

i^.^^..(M^), (4.1) 
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where the probabihty for the photon to develop a fluctuation with mass M is given by 



and (t(M^) is the cross section for this fluctuation to scatter off the target. The above 
expression for dP{M'^) is most easily motivated in the framework of old-fashioned per- 
turbation theory, where the energy denominator of the amplitude is proportional to the 
off-shellness of the hadronic fluctuation, -\- M^. If this is the only source for a 
dependence, the numerator factor is necessary to obtain a dimensionless expression. 
(See the original derivation of Gribov [5T| for more details.) 

Bjorken and Kogut assume that, for large M^, the intermediate hadronic state typ- 
ically contains two jets and that cr(M^) is suppressed for configurations with high p_L 
(the latter effect being now known under the name of colour transparency). Consider 
hadronic fluctuations with a certain M^, which, in their respective rest frames, are re- 
alized by two back-to-back jets. Under the assumption that the probability distribution 
of the direction of the jet axis is isotropic, simple geometry implies that aligned con- 
figurations, defined by p\ < (where is a soft hadronic scale), are suppressed by 
A^/M^. If only such configurations are absorbed with a large, hadronic cross section, the 
relations a{M'^) ~ l/M^ and 

dM2 ^ (M2 + Q2)2 ^ • ' 

follow. Thus, the above cross section can be interpreted as the total high-energy cross 
section of target proton and aligned jet fluctuation of the photon, i.e., of two soft hadronic 
objects. Therefore, a similar elastic cross section is expected, cr^ ~ cr^ (cf. Fig. ^?l|b). 
The resulting diffractive structure function, as defined in the previous chapter, reads 

F^^^\i,P,Q')r^^. (4.4) 

It is interesting that the very simple arguments outlined above result in an expression 
for Fj^''^'* that captures two important features of the HERA data: the leading-twist 
nature of diffraction and the approximate 1/^ behaviour. The main problems of the 
model are the precise energy dependence, which is measured to be somewhat steeper 
than 1/^, and the limit /5 — * 0, where a constant or even rising behaviour of ^2^*' is 
observed. 

The above discussion shows, in very simple terms, what was also one of the main 
qualitative results of the more technical treatment of the last chapter: diffractive elec- 
troproduction is based on a soft hadronic high-energy cross section. Such cross sections 
are very successfully described within the framework of pomeron exchange P]. Let us 
recall the basic underlying concepts (see |^ for a detailed discussion or for a brief 
introduction). 

Using analyticity and crossing symmetry, the amplitude T'i2^34(s, t), depicted in 



Fig. |4^, can be related to the amplitude Ti2,^2a{s' it')-, where s' = t, t' = s, and bared 
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Figure |4.2|: Scattering process 12 — > 34 via reggeon exchange. 



numbers denote antiparticles. The partial wave expansion for this crossed amphtude 
reads 

oo 

Tr3^-24is',t') = Y.{21 + l)aiis')Pi{cos9) , (4.5) 

1=0 

where 6 is the centre-of-mass frame scattering angle, which is a function of s', t' and the 
particle masses, and Pi are Legendre polynomials. Let the two functions a^(/,t) with 
rj = +1 and ?7 = — 1 be the analytic continuations to complex / of the two sequences 
{ai{t), I = 0, 2, 4, ...} and {ai{t), / = 1, 3, 5, ...}. In the simplest non-trivial case, the only 
singularity of a.r^{l,t) is a single t-dependent pole at / = a{t). It can then be shown that, 
in the limit s — ^ oo, 

/ e \ Ol{t) 

Tu^s4{s,t) = f3is{t)f324{t)Ua{t)) i^-J , (4.6) 
where sq is an arbitrary scale factor, /5i3 and are two unknown functions of t, and 

CMt)) = ' , , (4.7) 
^ sin7ra(t) ^ ^ 

is the signature factor, depending on the signature rj of the relevant Regge trajectory 
a(t). If a^(/,t) has a more complicated analytic structure, the rightmost singularity in 
the / plane dominates the behaviour at large s. 

Within the present context, the essential predictions of the asymptotic expression 
Eq. ( ^.6| ) are the power-like energy dependence s"'-*^ and the factorization of the unknown 
t dependence into the two vertex factors /3i3 and ^24- This last feature, which underlies 
the graphic representation of reggeon exchange in Fig. [4.2|, is relevant if the same Regge 



trajectory governs different scattering processes. Note also that, for positive t = s' and 
integer /, a(t) describes the positions of poles of the physical amplitude Ti^^24{s' ,t')- 
Such poles are expected whenever an on-shell particle with appropriate mass = 
s' and angular momentum / can be created in the collision of 1 and 3. Indeed, most 
Regge trajectories pass through known physical states with mass rn? = t and angular 
momentum a{t). 

The Froissart bound [0] on the high-energy growth of total cross sections. 



c^tot < ^In^s/so), (4.^ 
mi 
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where is the pion mass and Sq is an unknown scale factor, imphes that a(0) < 1 
for all Regge trajectories. However, it was observed early on that a very good fit to pp 
and pp cross sections, which were measured to rise at high energy, could be obtained 



by assuming the dominance of a single pole with a(0) > 1 ||8^. The corresponding 
trajectory, originally introduced with a{0) = 1 |8^, is known as the pomeron trajectory 
(cf. the Pomeranchuk condition of \W^ )- 

Donnachie and Landshoff demonstrated that a large set of different hadronic cross 
sections can be fitted with an intercept ajp(O) = 1.08 0. In spite of the power-like growth 
of Eq. (|4.(j|), the predicted cross sections are so small that the Froissart bound is not 
violated below the Planck scale. It is then argued that unitarity is not a serious problem 
at all realistic energies. However, one should keep in mind that all following analyses 
based on a pomeron trajectory with a(0) > 1 are, strictly speaking, not self-consistent 
in the framework of Regge theory. In fact, it is likely that the rightmost singularity in 
the complex / plane is not a single pole but a cut, in which case many of the following 
results would have to be reconsidered. 

^ X 




Figure 4.3: Diffractive electroproduction via pomeron exchange. 



Given the universal success of the Donnachie-Landshoff pomeron, it is natural to 
apply the same trajectory to the quasi-elastic amplitude of diffractive electroproduction. 
Thus, the diagram in Fig. [4.1|b is interpreted in terms of pomeron exchange (cf. Fig. 



In this situation, the soft energy dependence of the pomeron and the measured t depen- 
dence of the proton-proton-pomeron vertex are naturally combined with the aligned jet 
model prediction of Eq. ( [4.3| ). Ignoring the vertex factor from the upper part of Fig. [4.3| , 
one can write 

j^^mi!^^ (4 9) 

The non-trivial energy dependence represents a clear improvement compared to Eqs. ( [4.3|) 
and ( ^.4| ). However, our ignorance of the upper vertex prevents an unambiguous predic- 
tion of the distribution. In fact, the obtained suppression at large M^, which cor- 
responds to the region of small (3, appears to be the main qualitative problem of the 
presented model. Furthermore, the treatment of the hard scale is rather naive in view 
of the impressive successes of QCD perturbation theory in the description of inclusive 
structure functions. 

As is shown in the next section, a better understanding of the dependence is 
possible in the framework of Regge theory. The obtained results allow for an interpreta- 
tion of the upper vertex in Fig. ^]3| in terms of a 7*-pomeron collision and thus for the 
application of QCD perturbation theory. 
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4.2 Pomeron structure function 



It is the purpose of this section to explain the interpretation of Fig. [4.3| in terms of a 
7*-pomeron colhsion, thereby providing the background for the ensuing partonic treat- 
ment of the pomeron. Before doing so, it is helpful to recall Mueller's generalization of 
the optical theorem and the resulting triple-pomeron interpretation of single diffractive 
dissociation. 

The optical theorem relates the two-particle total cross section to the imaginary part 



of the forward amplitude. This is illustrated by the first equality in Fig. [4.4| . Here y/s 
is the centre-of-mass energy of particles 12, i.e., the mass of the system X, and the 
discontinuity across the real axis is defined by 



disCsA(s) = A{s + it) - A{s - ie) . 



(4.10) 



The last equality in Fig. illustrates the high-energy limit, which is dominated by 
pomeron exchange. 
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Figure |4.4| : Graphic representation of the optical theorem for a two-particle total cross 
section. In the rightmost diagram, which is relevant in the pomeron-dominated high- 
energy limit, the cut is represented by a dotted line. 



Mueller's generalization 



(see also |^), which is illustrated by the first equality 
in Fig. [4.5| , relates a one-particle inclusive cross section to a six-particle amplitude. The 
derivation uses crossing symmetry to reinterpret the outgoing particle 3 as the incoming 
antiparticle 3. In very much the same way as for the usual optical theorem, completeness 
of the sum over X and unitarity are employed to relate the inclusive cross section to 



the discontinuity of the amplitude. Note, however, that in Fig. [4.5| the discontinuity is 
taken in M^, which is the squared mass of the system X and different from the variable 
s characterizing the original process with incoming particles 1 and 2. 



The rightmost diagram in Fig. is obtained in the double limit s/M"^ — » oo and 



— ^ oo (see also [0]). If s ^ and is much larger than any of the other 
variables, the process 12 3X is dominated by pomeron exchange, and the amplitude 
receives a factor Cr,(a(t)) (s/M2)"^'(*). The appearance of the ratio s/M^ is a non-trivial 
result of the relevant kinematics. The remaining dependence of the amplitude is given 
by the cut that is left after the two pomerons coupled to particles 2 and 3 are factorized. 
If the high-energy limit of this amplitude, which corresponds to oo, is again 

dominated by pomeron exchange, a factor (^^(a(0)) (M^)"^^°) results. The appearance of 
these three 'pomeron propagators' is illustrated by the three zigzag lines on the r.h. side 
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Figure |4.5| : Graphic representation of Mueller's generalization of the optical theorem. 
In the rightmost diagram, relevant for s/M"^ — oo and oo, the cut is only through 

the pomeron coupled to particle 1. 



of Fig. W^. The corresponding cut amplitude reads 



discT 



ajp(i) 



2\ aip{0) 



disc(G(t)C,(«(0)) ( — 1 /5ii(0)|, (4.11) 



where G{t) is the unknown triple-pomeron vertex. Since the t dependence of G{t) is 
measured to be weak and since, for the pomeron, rj = +1, the approximations G{t) ^ 
G(0) and ~ 1 can be used at small \t\. The resulting cross section for the process 
12 —>■ 3X is given by 



da 



1 



dtdM^ 167r2M2so"'^^^^^^' (m2 



|G(0)/?n(0)| f — j 



(4.12) 



This expression for the cross section suggests the following interpretation [Q. In- 
coming particle 2 radiates a pomeron carrying a fraction ^ = M^/s of its momentum. 
Then, this pomeron collides with particle 1 producing the diffractive state X with mass 
M. According to this interpretation, Eq. ([4.12|) can be rewritten as 



where 



da 

dt d^ 



IGtt^so 



2 f l-2ar(t) 



is the pomeron flux factor, characterizing the probability for the transition 2 



a,jp{M') = \G{0)(3nm (— 

V So 



2 \ ajp(0)-l 



(4.13) 

(4.14) 
3iP, and 

(4.15) 



is the total cross section for the collision of the pomeron and particle 1. Since the pomeron 
is not a real particle, the decomposition of Eq. ( [4.12|) into pomeron flux and total cross 
section is ambiguous. The definitions given here correspond to interpreting (323{t) and 
(r]{<y{t)) (s/M^)"^ in Eq. ( [4.11| ) as pomeron-particle-particle vertex and pomeron propa- 
gator respectively. The rest of the diagram corresponds to the pomeron-particle scattering 
amplitude, which is used to define the pomeron-particle cross section. 
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In their seminal paper on diffractive jet production in hadron-hadron scattering, In- 
gelman and Schlein exploited this picture to predict the rate of high-p^ jets within the 
diffractive final state X 0. They suggested that, according to the parton model, hard 
processes in the collision of particle 1 and pomeron can be described in terms of a convolu- 
tion of two parton distributions and a hard partonic cross section. The idea of introducing 
a parton distribution for the pomeron proved to be very successful phenomenologically 
and was widely used in subsequent analyses of hard diffractive processes. 

Donnachie and Landshoff applied this idea to the case of diffractive DIS [|10], |Tl|, 
where particle 1 is a highly virtual photon and can therefore be treated in perturbation 
theory. In this approach, HERA diffraction at the level of a few per cent of the total DIS 



cross section was predicted as early as 1987 [IQ]- In contrast to Ingelman and Schlein, who 



focussed on the gluon distribution of the pomeron, Donnachie and Landshoff introduced 
a quark distribution, which can be directly probed by the virtual photon. The resulting 
physical picture of diffraction is shown in Fig. 

X 




P' 



Figure |4.6| : Diffractive DIS via the pomeron structure function. 

In complete analogy to Eq. (|4.13|) , the diffractive cross section can be written as 

^ = Me,t)ar, (4.16) 

where, according to the conventional parton model, the photon-pomeron cross sec- 
tion is given in terms of the quark distribution gjp(/9, Q^) of the pomeron, ^ = 
(vre^/Q^) 2/5gjp(/5, Q^). Here f3 = x/^ is the fraction of the pomeron momentum car- 
ried by the struck quark, and the factor 2 is introduced to account for the antiquark 
contribution. The resulting diffractive structure function reads 

F2''^'\x,Q\^,t) = fjpi^,t)2Pqjpi(3,Q'), (4.17) 

and a non-trivial dependence is naturally expected on the basis of the Altarelli-Parisi 
evolution of the quark distribution. 

The above normalization of pomeron flux and pomeron-particle cross section is con- 
sistent with [Q. However, other normalizations are also frequently used (compare the 



conventions of [p|-[TT|, |90|, pT|] ) 



4.3 Diffractive parton distributions 

Loosely speaking, the analysis of diffraction in terms of diffractive parton distributions is 
equivalent to the analysis in terms of the partonic pomeron of Ingelman and Schlein, 'mi- 
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nus its Regge content' [^]. Diffractive parton distributions provide a direct, perturbative- 
QCD-based approach to the hard part of the process, without introducing any less well 
established non-perturbative concepts. 

The basic theoretical ideas are due to Trentadue and Veneziano, who proposed to 
parametrize semi-inclusive hard processes in terms of 'fracture functions' [n], and to 



Berera and Soper, who defined similar quantities for the case of hard diffraction [13| and 
coined the term 'diffractive parton distributions'. The following discussion is limited to 
the latter, more specialized framework. 

In short, diffractive parton distributions are conditional probabilities. A diffractive 
parton distribution dfP{y, ^, t) / dt describes the probability of finding, in a fast moving 
proton, a parton i with momentum fraction y, under the additional requirement that the 
proton remains intact while being scattered with invariant momentum transfer t and 
losing a small fraction ^ of its longitudinal momentum. Thus, the corresponding 7*p 
cross section can be written as 1031 



*<-«^'«''*""" - E r C^W^) . (4,1* 



c?^ dt i "'^ \ '^^ '^'^ 

where (3"(x, Q^, y)"'*'^ is the total cross section for the scattering of a virtual photon charac- 
terized by X and and a parton of type i carrying a fraction y of the proton momentum. 
The above factorization formula holds in the limit Q"^ ^ oo with x, and t fixed. 

At leading order and in the case of transverse photon polarization, only the quark 
distribution contributes. For one quark flavour with one unit of electric charge, the well- 
known partonic cross section reads 



vre^ 



aT{x,Q',yr*'^ = —5{l-y/x), (4.19) 



giving rise to the diffractive cross section 

da{x,Q\i,tr'P^P''' _ 27re2 xdf!^{x,i,t) 



(4.20) 



d^dt g2 ^^^t 

where the factor 2 is introduced to account for the antiquark contribution. 

The main difference to fracture functions is the requirement of a specific t transferred 
to the final state hadron. Since fracture functions are defined to include the t integration, 
a non-negligible contribution to the relevant cross section arises from processes where 
the outgoing hadron, in this case the proton, belongs to the current fragmentation re- 
gion. This contribution complicates the dependence of fracture functions, but not of 
diffractive parton distributions. Even if the t integration is performed, the above problem 
does not arise in diffractive kinematics, where ^ is small and the production of a very 
energetic forward proton in the fragmentation of the current is strongly suppressed. 



It is now obvious that Eqs. ( [4.16|) and (|4.17| ) of the last section form a special case of 



the present more general framework. There, Regge phenomenology and the assumption 
of a partonic pomeron lead to the specific form 

a:^^^^ = Me,t)/3gip(/5) (4.21) 
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for the diffractive quark distribution. However, the factorizing dependence on 13 = x/C, 
and the expression for the pomeron flux factor /jp(^, t) given in the last section have not 
been derived in the framework of QCD. 

The operator definition of diffractive parton distributions can be obtained as fol- 
lows [Q. Consider first a complex scalar field which can be written as 

(j){x) = [ dk \d{k)e-'^'' + 6^(A;)e^H , dk = , , , (4.22) 
J L J (27r)'^2fi;o 

where a and h are the annihilation operators of particles and antiparticles respectively. In 
this field theory, the number of particles with momentum in the interval d^k is measured 
by the operator 

dN = dk d\k)d{k) . (4.23) 

The conventional inclusive distribution f{y) of partons of the field is given by the 
number of particles in a momentum interval dk^ = P+dy, normalized to the size of the 
interval dy. Note that, in this section, the co-ordinate system is such that the proton and 
photon momenta have large plus and minus components respectively. This is customary 
for the discussion of DIS in the Breit frame and contrasts with the remainder of this 
review, where the proton rest frame is emphasized and, correspondingly, the photon 
momentum is defined to have a large plus component. 

The particle number is measured in the state of a hadron with momentum P, and 
integration over all transverse momenta k± is assumed. Thus, the explicit formula reads 

i2nf2PoS'iP'-P)fiy)dy = J d'k^ {P'\-^\P) , (4.24) 

where the prefactor {2TTf 2Po 5^{P' - P) is a result of the normalization of the hadronic 
state \P), 

{P'\P) = {27i)^2Po5'\P' - P) . (4.25) 

It can then be shown that, in terms of the fundamental field 0, the distribution of scalar 
partons reads 

f{y) = ^J rfx_e-^^^+^-/2(P|0t(o,x_, 0^)0(0, 0,0^)|P). (4.26) 
In the case of a gauge theory, this definition has to be supplemented with a link operator 

U,_,o = Pexp i^-'- £' A+(0,y_,Ox)rfi/_) , (4.27) 

connecting the two scalar field operators. 

Before generalizing this definition of the parton distribution to the diffractive case, 
it is convenient to rewrite it in the form 

f{y) = ^ / dx^e-'y^'+^-Z^ E(^l<^^(0' 0±)f/._,oo|X) (X|f/oo,o 0(0, 0, 0^)\P) . 

(4.28) 
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In the difFractive case, the operators describing the creation and annihilation of the 
parton are the same. However, the proton is required to appear in the final state carrying 
momentum P' . Thus, the above definition is changed to 

<(l/,e,t) yP^ 



, / rfx_e-*^^+^-/2 (4.29) 
X Y.^P\<P\Q, x_, Ox)?7,._,oo|P', X) (P', X|[/oo,o 0(0, 0, Ox) |P) , 

X 

where (P', X| denotes the outgoing state, the only restriction on which is the presence 
of a scattered proton with momentum P'. 

The above formulae for inclusive and diffractive distributions of scalar partons are 
generalized to the case of spinor quarks by the substitution 

0t(O, x_, 0^)0(0, 0, Ox) ;rV^(0' 0^)7+^(0, 0, 0^) , (4.30) 

and to the case of gluons by the substitution 

0t(O, a;_, Ox)0(O, 0, 0^) -^F^{0, x., Ox)+^P(0, 0, 0^)/ . (4.31) 

[yP+r 

The operator expressions appearing in the definitions of both the inclusive and 
diffractive parton distributions are ultraviolet divergent. They are conveniently renor- 
malized with the MS prescription, which introduces the scale as a further argument. 
The distribution functions then read /(x, /i^) and df^{x,^,t,fi'^)/d^dt. 

Accordingly, Eq. ( [4.18| ) has to be read in the MS scheme, with a fi dependence 



appearing both in the parton distributions and in the partonic cross sections. The claim 
that Eq. ( [4.18| ) holds to all orders implies that these /i dependences cancel, as is well 



known in the case of conventional parton distributions. Since the partonic cross sections 
are the same in both cases, the diffractive distributions obey the usual Altarelli-Parisi 
evolution equations, 

m-AhKl _ ^ f . (4.32) 



d{yn.jj?) d^dt j ^ y d^dt 

with the ordinary splitting functions Pij{x/y) 



Clearly, this is equivalent to the assertion that, in the operator definition of Eq. ( [4.29|) 



the ultraviolet divergences are independent of the final state proton P'. If this is the 
case, the Altarelli-Parisi evolution of the distribution functions follows from the operator 
definitions exactly as in the inclusive case of Eq. ( [4.28[ ). 



Thus, for the analysis of diffractive DIS, it is essential to gain confidence in the 
validity of the factorization formula Eq. ( [4.18| ). Berera and Soper first pointed out 



that such a factorization proof could be designed along the lines of related results for 
other QCD processes (see ||9^ for a review). 
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Using Mueller's method of cut vertices ||9^, Grazzini, Trentadue and Veneziano |97] 
proved, in the framework of a simple scalar model, that the above factorization property 
holds for 'extended fracture functions'. These objects differ from fracture functions in 
that they are differential in the momentum transfer t. They are thus equivalent to the 
diffractive parton distributions discussed here. 

Collins showed that factorization holds in full QCD by generalizing the essential 
step of dealing with soft gluon interactions, which lies at the heart of many previous 
factorization proofs |^, to the case of diffractive DIS. A brief summary of the essential 
arguments is presented below. 

To begin with, recall the main ideas of the factorization proof for inclusive DIS. 
It is based on the dominance of contributions from so-called leading regions, shown in 
Fig. [4.7| a. Here the hard subgraph is denoted by H, the subgraphs with momenta collinear 
to P and to the produced jets are denoted by A and Ji ■ • • J„ respectively, and the soft 
subgraph is denoted by S. The analysis is performed in the Breit frame, where lines in H 
have typical virtuality Q^, lines in A and Ji ■ ■ • have small virtualities but may have 
large longitudinal momentum components of order Q, and all components of momenta 
in S are small compared to Q. 




a) 




Figure 4.7 : Leading regions in inclusive (a) and diffractive (b) DIS (cf. ||98|| ). 



The rationale behind the discussion in terms of leading regions is the realization 
that the momentum integrations are dominated by pinch singularities, i.e., poles of the 
propagators that can not be avoided by deforming the integration contours. In the limit of 
large Q^, such singularities are associated with the leading regions shown in Fig. \i.7\ . By 
power counting, only one hard line may connect the A and H subgraphs, while arbitrarily 
many soft (dashed) lines may connect the soft subgraph with other parts of the diagram. 

An essential part of the factorization proof is the demonstration that the soft sub- 
graph can be factored out. More specifically, it has to be shown that the soft lines are not 
important for the subgraphs H and Ji ■ ■ ■ J„, so that a perturbative hard cross section 
with free partons in the final state can be used. This is achieved using Ward identi- 
ties, which, however, can only be applied in the region where all components of the soft 
momenta are small and of comparable size. Difficulties arising in the region where one 
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component of a soft momentum is much smaller than the other components can be solved 
by appropriately deforming the integration contour. 

To see this in more detail, consider the particularly simple diagram of Fig. ^.8| , 
where there is only one current jet, and a single soft gluon connects the corresponding 
jet subgraph with subgraph A. Let the soft gluon with momentum / couple to outgoing 
particles with momenta kj and in the subgraphs J and A respectively. The particle 
propagators 

and -r^ ^— (4.33) 



(/cj — ly — m? + ie (kA + 0^ " + 

attached to the gluon vertices produce poles in the complex /+ and /_ planes that lie above 
and below the real axis respectively. This is also true for further / dependent propagators 
in J and A since, due to completeness, final state interactions can be disregarded in both 
subgraphs. Thus, regions where either Z+ or /„ are too small can always be avoided by 
deforming the integration contour. This takes us back to the genuinely soft region where 
all components of / are small and of comparable size, and where Ward identities can be 
used to factor out the soft subgraph. 




Figure |4.8| : Leading order process in DIS with additional exchange of a soft gluon 



between the target and current jet subgraphs (cf. @]). 

It is precisely this part of the factorization proof that is affected in the case of 
diffraction. The requirement of a final state proton leads to the presence of final state 
interactions in subgraph A of Fig. [4.8| (as shown in Fig. [4.7| b) that can not be neglected. 
Since both initial and final state interactions appear in A, poles on both sides of the /_ 
plane exist, and the contour can no longer be deformed to avoid the dangerous region 



of too small /_. However, Collins was able to demonstrate |^ that deformations of the 
/+ contour can instead be used to show that this dangerous region does not produce a 
non-factorizing contribution. The argument is based on a detailed analysis of the pole 
structure in subgraphs Ji ■ ■ ■ Jn (or the single jet subgraph J in the simple example 
above). It leads to the conclusion that, exactly as in the inclusive case, the soft subgraph 
can be factored out. The result is a convolution of the calculable hard part with diffractive 
parton distributions. 

Note that diffractive factorization does not hold in the case of hadron-hadron col- 
lisions. A general argument for this was given in and the effect was also found in 



the model calculation of [13]. The reason is that, in contrast to the DIS case, the use of 
completeness in the final state can not be avoided in the factorization proof for hadron- 
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hadron collisions. This completeness is lost if a final state proton with a given momentum 
is required, and a breakdown of the factorization theorem results. 



4.4 Target rest frame point of view 

In this section, the connection between the target rest frame point of view, used in the 
semiclassical approach, and the Breit frame point of view, relevant for the two previous 
sections, is established. In particular, the consistency of the semiclassical approach with 
the concept of diffractive parton distributions and with the factorization formulae of the 
last section are demonstrated. 



Even before the all-orders factorization proofs of |^ and [^], it was suggested that 
diffractive factorization could be understood in the semiclassical picture in the proton 
rest frame ||100|| . This was then explicitly shown in the leading order analysis of |^ 



on 



which the present section is largely based. Calculating the cross section with the methods 
of Chapter |^, a result is obtained that can be written as a convolution of a partonic 
cross section and a diffractive parton distribution. Within the semiclassical model, this 
diffractive parton distribution is explicitly given in terms of integrals of non-Abelian 
eikonal factors in the background field. 

To be more specific, it is explicitly shown that the amplitude contains two fundamen- 
tal parts: the usual hard scattering amplitude of a partonic process, and the amplitude 
for soft eikonal interactions with the external colour field. The latter part is determined 
by the scattering of one of the partons from the photon wave function. This parton has 
to have small transverse momentum and has to carry a relatively small fraction of the 
longitudinal photon momentum. In a frame where the proton is fast, this parton can be 
interpreted as a parton from the diffractive structure function. 

The special role played by the soft parton in the photon wave function was also 
discussed in |28, pT, 101 , 102|| in the framework of two gluon exchange. However, the 



present approach has the two following advantages: firstly, by identifying the hard part 
as a standard photon-parton scattering cross section, the necessity for non-covariant 
photon wave function calculations is removed. Secondly, once it is established that the 
main contribution comes from the soft region, non-perturbative effects are expected to 
become important. The eikonal approximation provides a simple, self-consistent model 
for this non-perturbative region. 

The explicit calculation closely follows the calculation of Chapter |^. It is convenient 
to start with the particularly simple case of scalar partons. The kinematic situation is 
shown symbolically in Fig. [4.9| . The process is split into two parts, the hard amplitude 
for the transition of the photon into a virtual partonic state and the scattering of this 
state off the external field. 

To keep the amplitude for the first part (denoted by H) hard, the transverse momenta 
P[j)± ~ l...n) are required to be large, i.e., ~ Q. The momentum k± is small, i.e., ~ A, 
and the corresponding parton carries only a small fraction (~ h? /Q"^) of the longitudinal 
photon momentum in the proton rest frame. While the hardness condition for particles 1 
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proton 



'(1) 



'(n) 



Figure |4.9| : Hard diffractive process in the proton rest frame. The soft parton with 
momentum k is responsible for the leading twist behaviour of the cross section. 



through n is introduced "by hand" , simply to make the process tractable, the softness of 
the last particle follows automatically from the requirement of leading twist diffraction. 

The cross section for the scattering off a soft external field reads 

da = ^ |T|2 27r5(go - go) dX^''^^^ , where q' = k' + Y.p'u) ■ (4.34) 



2g, 







All momenta are given in the proton rest frame, T is the amplitude corresponding to 
Fig. [4.9| , and dX^"^^^^ is the usual phase space element for n + 1 particles. 



According to Eq. ( |3.3| ), each of the particles interacts with the external field via the 
effective vertex 

V{p\p) = 2^5{p',-po)2p, [U{p'^-p^) -{27,f 5\p'^-pJ\ . (4.35) 

The amplitude T can now be built from the hard amplitude (symbolized by H in 
Fig. |4^ ) together with the effective vertices defined by Eq. ( |4.35| ) and the appropriate 
propagators. It is convenient to consider first the amplitude T', which is defined like T 
but without the transverse-space 5-function terms of Eq. ( 4.35| ), 



X 



^ 27r5(A;[, - A;o) 2ko U{k'^ - k^) • (4.36) 



In this equation, colour indices have been suppressed. Notice also that some of the 
produced partons are antiparticles. The corresponding matrices U have to be replaced 
by f/^ To keep the notation simple, this is not shown explicitly. 

The integrations over the light-cone components P(j)+ can be performed using the ap- 
propriate energy 5-functions. After that, the P{j)- integrations are performed by picking 
up the poles of the propagators I/Pq)- The result is 

d'^V( i\ I \ 2kn ~ 



T'= lTn\{[ Uip'u)^ - PU)±) ^) 9^ U{k'^ - k^) . (4.37) 
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Here poles associated with the P{j)- dependence of Tjj have been disregarded since their 
contribution is cancelled by diagrams where part of the hard interaction occurs after the 
scattering off the external field. 

Next, a change of integration variables is performed, 

d'^P{n)L d^k± . (4.38) 

Since the external field is assumed to be soft, it can only transfer transverse momenta 
of order A, i.e., P{j)± — p[j)± for all J- In general, the amplitude Th will be dominated 
by the hard momenta of order Q. Therefore, it can be assumed that Th is constant 
if the momenta P{j)± vary on a scale A. In this approximation, the integrations over 
Pij)± (i = l...n — 1) can be performed in Eq. ( [4.37|) , resulting in (5-functions in impact 



parameter space. These manipulations give the result 

= lB?l^^''Ly^ (nf^(^±))f^(^^+l/±)e-'^^^^-^^-(^--'^-) , (4.39) 

where A is the total momentum transferred from the proton to the diffractive system 
and, in particular, = k'j_ + J2p[j)±- It is intuitively clear that the relative proximity 
of the high-px partons in impact parameter space leads to the corresponding eikonal 
factors being evaluated at the same position x±. 

Now, the colour structure of the amplitude will be considered in more detail. Spelling 
out all the colour indices and introducing explicitly the colour singlet projector P, the 
relevant part of the amplitude reads: 

Zoiour = n'-''"' (llU{x±)y"'^" U{x^ + y^t P,;...,, . (4.40) 

\ j / ai...an 

Using the fact that Th is an invariant tensor in colour space, and introducing the function 
W defined in Eq. (|3.17| ), the following formula is obtained: 

rp rjii rpai...anbr> rpai...a„b Tj^b' p (A AT\ 

-I- colour — J- colour ^ai...anb — H b -< ai...a„6' • l^-^-^J 

Here the first equality states the explicit relation between T' and the true amplitude T, 
where the trivial contribution of zeroth order in A has been subtracted. This subtraction 
corresponds to the unit matrix on the r.h. side of Eq. ( |3.17| ). 

For colour covariance reasons 

T^d-'-''Pa,...a^b' = const. X d^, . (4.42) 
Since the photon is colour neutral, the following equality holds: 

T^H^-""'r^a....a„6 = \Tr-'''''Pa,...a^b? = . (4.43) 

Here the partons are assumed to be in the fundamental representation of the colour 
group SU{Nc). Combining Eq. ( |4.41| ) with Eqs. ( [4.42| ) and (|4.43| ), it becomes clear that 



the colour structure of the hard part decouples from the eikonal factors, 

1 

N, 



Tcolour\' = ^\tT[W]\' \Th\'. {AAA) 
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The hard part will be interpreted in terms of an incoming small-Zc^ parton that collides 
with the virtual photon to produce the outgoing partons 1 through n. Therefore, a factor 
1/Nc for initial state colour averaging is included in the definition of iT/zp. 

In the expression for the cross section, the two functions W appear in the combination 



(tr[^.Jyx)])(tr[W^.^(yl)])*e- 



(4.45) 



with independent integrations over Xj^,x'^,y_[_ and y'^. If the external field is sufficiently 
smooth, the functions W vary only slowly with x_l and x'^. Therefore, after integration 
over x_[_ and x^, the expression in Eq. (|4.45|) produces an approximate 5-function in Aj^. 
Furthermore, it is assumed that the measurement is sufficiently inclusive, i.e., the hard 
momenta are not resolved on a soft scale A. This corresponds to a A_L-integration, 
which gives an approximate 5-function in x±—x'^. Since the expression in Eq. ( |4.45| ) will 
always appear under x^, x'^ and A_l integration, the above considerations justify the 
substitution 

g-i(xx-x'Y)Ax _ (27r)2 6^{x^ - xl) 6\^i_) . (4.46) 

Combining Eqs. ( [4.34| ), ( |4.39| ) and ( [4.44| ), the following formula for the cross section 
results, 



da= — 
2go 



IT, 



H\ 



<Pk^ tr[iy.Jfcl-fcx)] 
(27r)2 



5(go-go)rfX("+i). 

(4.47) 

Note that the soft momentum k'j^ has been neglected in the transverse 5-function. 

To finally establish the parton model interpretation of diffraction, the hard partonic 
cross section based on iT^-p has to be identified in Eq. ( [4.471 ). Consider the process 

7*(g) + q{yP) ^ g(p'(i)) + . . . + g(p'(„)) , (4.48) 

where the photon collides with a parton carrying a fraction y of the proton momentum 
and produces n high-p^ final state partons. The cross section is approximately given by 



da{y) 



1 



2(s + Q2 



(4.49) 



where s = (SPq))^ and the quantities |Tf^P and k are the same as in the previous discus- 
sion. Equation ( [4.49|) is not exact for several reasons. On the one hand, |TffP is defined in 
terms of the unprimed momenta pq), which differ slightly from p'^jy On the other hand, 
the vector k is slightly off shell and has, in general, a non-zero transverse component. 
However, both effects correspond to A/Q corrections, where Q stands generically for the 
hard scales that dominate Th- 



Using Eq. ( 4.49 ), the cross section of Eq. ( 4.47 ) can now be rewritten as 



da 



dk- 



27rgo 



i2kofda{y) 



d^k^ tijW.^k'jrk^)] 
(2vr)2 



d^k' 



{2n)^2k'o 



(4.50) 



The light-cone component k_ is given by — A;_ = yP_ = yrrip. Note that the minus 
sign in this formula comes from the interpretation of the parton with momentum k 
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as an incoming particle in Eq. (|4.49|) . This is, in fact, the crucial point of the whole 
calculation: due to the off-shellness of k, the corresponding parton can be interpreted as 
an incoming particle in both the process of Eq. ( |4.48|) and in the soft scattering process off 
the external field. The latter process, where an almost on-shell parton with momentum 
k scatters softly off the external field changing its momentum to the on-shell vector k', 
is most easily described in the proton rest frame. By contrast, the natural frame for the 
hard part of the diagram is the Breit frame or a similar frame. In such a frame, the 
k_ component is large and negative, so that the above parton can be interpreted as an 
almost on-shell particle with momentum —k, colliding head-on with the virtual photon. 



Substituting the variables y and ^ for k_ and k'^, the cross section, Eq. ( [4.50|) , takes 
the form 

where the diffractive parton distribution for scalars is 



(4.51) 




d^k± tT[W^^{k'j_-k^)] 



(27r)2 k'lP + kl{l 



(4.52) 



This result is in complete agreement with the concepts described in the last section. In 
contrast to the formulae presented there, the above parton distribution is inclusive in 
t. However, as will be seen in an example calculation in Sect. |6.1| , parton distributions 
differential in t can also be obtained in this formalism. 

For an external colour field that is smooth on a soft scale A and confined to a 
region of approximate size 1/A, the function tr[P1/'a;j^(?/x)] is also smooth and vanishes at 
y± = together with its first derivative. From this it can be derived that the k± and fc^ 
integrations in Eq. ( |4.52| ) are dominated by the soft scale. This justifies, a posteriori, the 
softness assumption for one of the partons used in the derivation. 

The qualitative result is that the eikonal scattering of this soft parton off the proton 
field determines the diffractive parton distribution. The hard part of the photon evolution 
can be explicitly separated and expressed in terms of a standard cross section for photon- 
parton collisions. 

Having worked out the kinematics in the simple scalar case, it is straightforward to 
extend the calculation to realistic quarks and gluons. The introduction of spinor or vector 
partons does not affect the calculations leading to the generic expression in Eq. ( [4.51|) . 
However, those parts of the calculation responsible for the specific form of Eq. ( 4.52 ) 
have to be changed if the soft parton is a spinor or vector particle. 

Referring the reader to Appendix y for details, only the final formulae for diffractive 
quark and gluon distributions in the semiclassical approach are presented below. They 
read 



2 



d^k'^{k 



(27r)4Ar, 



d'^k^ k^ii[W^^{k'^k^)\ 



(27r)2 k'll3 + kl{l-l3) 



(4.53) 
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for the case of a realistic spinor quark, and 



di 





{2nr (iV2 - 1) 7. 



(27r)2 k'l(3 + kl{l-(3) 



with 



2k'' 



(4.54) 



(4.55) 



for the case of a gluon. It was checked exphcitly that these distribution functions to- 



gether with the appropriate partonic cross sections reproduce the resuhs of Sects. |0 
and p73| . Henceforth, the simphfied notation dq{y,$,)/d^ and dg{y,^)/d^ will be used for 
the diffractive quark and gluon distribution calculated above. 

The application of these diffractive quark and gluon distributions to diffractive DIS 
is summarized in an intuitive way in Fig. 4.10. On the l.h. side, the two lowest-order 



processes, qq and qqg state production, are shown from the target rest frame point of view. 
On the r.h. side, the same two processes are shown from the Breit frame point of view, 
which is conventionally used to discuss the partonic interpretation of DIS. The essence 
of the calculations presented in the present section is the identification of these two 
viewpoints, as a result of which explicit formulae for the diffractive parton distributions, 
expressed in terms of the target colour field, are obtained. ^ 




a) 




Rapidity 
gap 




Figure |4.1C)| : Diffractive DIS in the proton rest frame (left) and the Breit frame (right); 



asymmetric quark fluctuations correspond to diffractive quark scattering, asymmetric 
gluon fluctuations to diffractive boson-gluon fusion. 



Note that it is incorrect to interpret the r.h. side of Fig. |4.10| as a two-step process, 
where a colour-neutral cluster is first emitted by the proton and then probed by the 
virtual photon. If this was the case, the two-gluon or two-quark cluster relevant in this 
calculation would necessarily lead to parton distributions symmetric in /3 and 1 — j3. A 
counter example to this is provided by the model distributions derived in Sect. |6.2|. 
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< > 



b) 




Figure |4.11| : Inclusive DIS in tlie proton rest frame (left) and the Breit frame (right); 
asymmetric fluctuations correspond to quark scattering (a), symmetric fluctuations to 
boson-gluon fusion (b). 



Thus, the idea of a 'pre-formed' colour neutral cluster, as it is usually associated 
with the pomeron structure function, is not supported by the present calculation. 

At this point, it is appropriate to add a brief discussion of the target rest frame vs 
Breit frame interpretations of inclusive DIS The leading order semiclassical calcu- 
lation of inclusive DIS, which amounts essentially to inclusive qq pair production off an 
external fleld, was given in Sect. p.2| . It was shown that, in contrast to the diffractive case, 
both asymmetric and symmetric qq conflgurations contribute to the leading twist cross 
section. As explained in more detail in Appendix 0, these conflgurations correspond, in 
the parton model, to leading order quark scattering, testing the inclusive quark distribu- 
tion, and boson-gluon fusion, testing the inclusive gluon distribution (cf. Fig [4.11j ). The 



symmetric conflgurations have a small transverse size and test directly the one-gluon 
component of the target colour fleld. This is the reason why, in the semiclassical frame- 
work, the leading order calculation is already sensitive to the gluon distribution of the 
target. 

The inclusive gluon distribution plays a very special role. In contrast to both the 
inclusive quark distribution and the diffractive quark and gluon distributions, it is only 
sensitive to the short distance structure of the proton fleld, and it is enhanced by an 
explicit factor l/a^ (see Appendix 0). As a result, the dominance of the inclusive over 
the diffractive DIS cross section, which is of fundamental importance for the successful 
phenomenological analysis of |^ (cf. Sect. |7. Ij ) , emerges. 
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5 Two Gluon Exchange 



So far, inclusive diffraction, as parametrized, e.g., by the diffractive structure function 
F^, was at the centre of interest of this review. It was argued that, for inclusive processes, 
the underlying colour singlet exchange is soft, and two corresponding approaches, the 
semiclassical framework and the pomeron picture, were described in some detail. 

In perturbative QCD, the simplest possibility of realizing colour singlet exchange is 
via two t channel gluons. In fact, the colour singlet exchange in certain more exclusive 
diffractive processes is, with varying degree of rigour, argued to be governed by a hard 
scale. In such cases, two gluon exchange dominates. Some of these processes are discussed 
in the present chapter. Finally, attempts to approach the whole diffractive cross section 
in two-gluon models are described. 



5.1 Elastic meson production 

Elastic meson electroproduction is the first diffractive process that was claimed to be 
calculable in perturbative QCD llSip!^. It has since been considered by many authors, 
and a fair degree of understanding has been achieved as far as the perturbative calcu- 
lability and the factorization of the relevant non-perturbative parton distributions and 
meson wave functions are concerned. 

To begin with, consider the electroproduction of a heavy qq bound state off a given 
classical colour field. This calculation represents an alternative derivation Ryskin's cele- 
brated result ||15[ for elastic J/ ip production. 

The relevant amplitude is shown in Fig. |5.1| . In the non-relativistic limit, the two 
outgoing quarks are on-shell, and each carries half of the J/i/j momentum. Thus, the two 
quark propagators with momenta p' = k' = g'/2 and the J/ip vertex are replaced with 
the projection operator gj^ jij^' + m). Here 

^ "em 

Fgg is the electronic decay width of the J/ip particle, mj = 2m is its mass, and ej its 
polarization vector ||103| . 



Using the notation and calculational technique of Sect. |3^ , the amplitude of Fig. |0 
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Figure 5J.: Leading order amplitude for the elastic production of a J/\l/ particle off an 
external colour field. 



can now be written as 
i27iS{qQ - qo)Tqg = ie^gj 



(M 



rtr 



'j{^' + m)V,{p',p) 



j/ ~ m — m 



-V,{k, k') 



(5.2) 

where Cc = (2/3)e is the electric charge of the charm quark. The Dirac structure is 
simplified employing the identities 



gj^ji^' + m) 



gj 

2m 



gj_ 

2m 



Y^Vr{k')Vr{k')i,j Y.Us{p')Us{p') (5.3) 



as well as Eqs. (|3.1CI|) and (|3.11| ). The further calculation proceeds along the lines of 
Sect. 3^ using, in particular, the quark scattering vertices Eqs. ( 3.12|) and ( 3.13|) . Adding 
the two contributions where only the quark or only the antiquark is scattered to Eq. ( ^.2|) , 
the full amplitude T = Tqq + Tq + Tq takes the form 



T 



T'dcg.j 

Sn'^m 



q+ / (i^k^ 



a(l — a) 



E [uApyjiq)vr'Ck)] [vrik'yjiq)uAp')] (5.4) 



X 



+ fci ^ 

U{p'±-p±)UHk± - k'^) - (27r)V(pl -pA_)5\k'^ - k^) 



Note that this is very similar to Eq. ( p.l4|) with the U matrix structure replaced according 
to Eq. ( ^.l(j| ). The main difference is that the produced quarks are projected onto the 
J/ip state. 

Since both Q and m are considered to be hard scales while U and f/^ are governed by 
the soft hadronic scale A, the integrand in Eq. ( |5.4|) can be expanded in powers of the soft 
momentum k±. The leading power of the amplitude is given by the first non- vanishing 
term. In the case of forward production, p'^ = /c^ = 0, the dependence on the external 
colour field takes the form 

(fk^klii \u{p'^-pA_)U\k^ - k'^) - {2TxY5^{p'^~p^)5^{k'^ - k^) 



d^kik^, 



tr 



Uix^)U^iy^) - 1 



oik± {y±-'J:±) 



J Xj_ 



y±=o 



(5.5) 



Now, the crucial observation is that precisely the same dependence on the external 
field is present in the amplitude for forward Compton scattering shown in Fig. p.2| . In 
the case of longitudinal photon polarization, the transverse size of the qq pair is always 
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Figure 



The Compton scattering amplitude witliin the semiclassical approach. 



small, and the target field enters only via the second derivative of W that appears in 
Eq. (pD. 

Thus, comparing with the parton model result for longitudinal photon scattering, this 



derivative can be identified in terms of the gluon distribution of the target proton ||4^ 



d'y^ I tiW^^iy^) ^^^^ = 27r'asxg{x) . (5.6) 



Using this relation^] and Eq. (O), the amplitudes for the forward production of trans- 



versely and longitudinally polarized J/ijj mesons by transversely and longitudinally po- 
larized virtual photons are obtained. To go from amplitudes for scattering off an external 
field to usual covariant amphtudes, a factor 2mp has to be introduced. Under the addi- 
tional assumption ^ mj, the covariant amplitudes for longitudinal and transverse 
polarization read 

Tl = -iQ47r^asgje{xg{x))-^ , Tt = ^Tl, (5.8) 

where ^/s is the centre-of-mass energy of the 7*p collision. 

It is not surprising that the gluon distribution of Eq. (p.6|), calculated according 
to Fig. f).2\ shows no scaling violations and only the trivial Bremsstrahlungs energy 
dependence ~ The reason for this is the softness assumptions of the semiclassical 
calculation. Firstly, the eikonal approximation implies that all longitudinal modes of the 
external field are much softer than the photon energy. Secondly, the reduction of the 
field dependence to a transverse derivative is only justified if the scales governing the 
quark loop, i.e., in the case of Fig. 5^ and and m? in the case of Fig. are 



harder than the transverse structure of W . These two approximations, evidently valid 
for a given soft field, are also justified for a dynamical target governed by QCD as long 
as only leading logarithmic accuracy in both 1/x and is required. Thus, a non-trivial 
dependence on and can be reintroduced into Eq. (|578|) via the measured gluon 
distribution, keeping in mind that the result is only valid at double-leading-log accuracy. 

Note that this is precisely what was claimed in the original two-gluon exchange 



calculation of |15[. Note also that, in distinction from this presentation, the calculation 



of [|l5l obtains the gluon distribution by coupling the t channel gluons to the quarks of 



^Note that this is also consistent with the resuh for longitudinal photon scattering given in Eq. ( |3.25 ) 



since 



dlj trW^Av^) tr[dy^W^A^)dy^WlAQ)\ , (5.7) 

i/j_=0 



which follows from unitarity of the U matrices. 



60 



the target and identifying the logarithmic integral over the transverse gluon momentum 
as the logarithm accompanying the usual quark-to-gluon splitting function. 

The first essential extension of the above fundamental result is related to the treat- 



ment of the bound state produced. Brodsky et al. fl^] showed that, at least for longitu- 
dinal photon polarization, a perturbative calculation is still possible in the case of light, 
non-perturbative bound states like the p meson. The calculation is based on the concept 
of the light-cone wave function of this meson. Referring the reader to | |104| j for a detailed 
review, a brief description of the main ideas involved is given below (cf. [ [105|| ). For this 



purpose, consider the generic diagram for the production of a light meson in a hard QCD 



process given in Fig. 5.3 



k 
q-k 

Figure |5.3|: Generic diagram for meson production in a hard process. 



Assume that, as shown in this figure, all diagrams can be cut across two quark-lines, 
the constituent quarks of the meson, in such a way as to separate the hard process H 
from the soft meson formation vertex V, which is defined to include the propagators. 
The amplitude can be written as 

T = j (fkTH{k)V{k) = dzTniz) ^ J dk-d^k±V{k) = dzTH{z)(j){z) , (5.9) 

where z = k^/q'_^_, and the last equality is simply the definition of the light-cone wave 
function of the meson. The two crucial observations leading to the first of these equal- 
ities are the approximate k_ and k±_ independence of Th and the restriction of the z 
integration to the interval from to 1. The first is the result of the hard scale that 
dominates Th, the second follows from the analytic structure of V. In QCD, the k±_ 
integration implicit in ip usually has an UV divergence due to gluon exchange between 
the quarks. Therefore, one should really read = 0(z,yU^), where the cutoff /i^ is of the 
order of the hard scale that governs T^- At higher orders in a^, the hard amplitude Th 
develops a matching IR cutoff dependence. 

A more rigorous definition of the light-cone wave function, required, in particular, for 
the discussion of higher order corrections, can be given in the operator language. Without 
going into further detail, note that the above wave function satisfies the relation 

{meson{q')\ip\y)ifi{-y)\0) = C dz e'(i-2^)(«'j^)0(z) , (5.10) 







where = and (p is the field operator, for simplicity a scalar, corresponding to the 
particle content of the meson. 

The case of the p meson is more complicated because it is a vector particle built from 
spin-(l/2) constituents. The matrix element relevant for the exclusive electroproduction 
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of a longitudinally polarized vector meson reads 



Jo 

(5.11) 

where e is the polarization vector of the meson and fa^b,c are leading twist qq distribution 
functions (see ||106|| for more details on these quantities and related issues). After the 



convolution with the appropriate hard production amplitude, only the combination 

i d 
2d'z 



^(z) = U{z) + - — Uz) (5.12) 



survives. Under the simplifying assumption that the p is built from one flavour of quarks 
with one unit of electric charge (which does not affect the flnal result if normalized by 
the decay width) the longitudinal amplitude reads |T^,|107| 



n = -tAn'ase (xgix))-^ f ^z^f^ . (5.13) 

SQ-^nip Jo Z{1 — z) 

The underlying calculation of the hard scattering process is analogous to the case of the 
J 111) discussed above. Although the shape of the light-cone wave function 0p is not known, 
its normalization can be related to the electronic width F^g of the p via the relations 

^ dz<Pp{z) = V2mpfp , r^, = ■ (5.14) 

If the non-relativist ic approximation used in the J/ip case was to be applied here, a wave 
function 0p proportional to 6{z — 1/2) would result. In this case, Eqs. (|5.13|) and ( |5.14|) 



would give an unambiguous prediction for the p production amplitude in terms of its 
electronic width, precisely as for the J/ ip amplitude. 

Note that, in the transverse case, the hard amplitude generates a z dependence which 
is even more singular at the endpoints than the factor l/z{l — z) of Eq. ( p.l3| ). While the 
expected fall-off of the wave function at the endpoints is sufficient to render Eq. ( |5.13|) 
flnite, this is probably not true in the transverse case. Thus, perturbative calculability 
in the sense of the J/ip and the longitudinal p amplitudes can not be claimed. 

Since the gluon distribution is deflned by the imaginary part of the forward Compton 
scattering amplitude, the above result for longitudinal vector meson production, ~ 
ixg{x), is, strictly speaking, only the imaginary part of the full amplitude. Assuming 
that the energy dependence is given by xg{x) ~ (1/x)", where the intercept a belongs to 
an even signature trajectory, the real part follows from the signature factor of Eq. ( [4.71) . 
If a — 1 is small, the ratio of real and imaginary part is approximately (7r/2)(a — 1), and 
the real part correction can be introduced into Eq. (|5.13| ) via the substitution ||T6|,|lQ 



• / N • / N ^ dixgix)) ,^ ^, 

%xg{x) ixg{x) + - ^^^^^^^^ ■ (5-15) 

The discussion of vector meson production given above was limited to the double- 
leading-log approximation as far as the colour singlet exchange in the t channel is con- 
cerned. To go beyond this approximation, the concept of 'non-forward' or 'off-diagonal' 
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parton distributions, introduced some time ago (see |18,21| and refs. therein) and dis- 
cussed by Ji il^ and Radyushkin [pOl in the present context, has to be used. Recent 



reviews of these quantities and their evolution, which interpolates between the Altarelli- 
Parisi and the Brodsky- Lepage evolution equations, can be found in |p.09| . 



Recall first that the semiclassical viewpoint of Figs. |5.1| and p.2| is equivalent to two 
gluon exchange as long as the transverse size of the energetic qq state is small. So far, 
the recoil of the target in longitudinal direction has been neglected. However, such a 
recoil is evidently required by the kinematics. For what follows, it is convenient to use 
a frame where q_ is the large component of the photon momentum. In Fig. f)A\ the 
exchanged gluons and the incoming and outgoing proton with momenta P and P' are 
labelled by their respective fractions of the plus component of P = (P + P')/2. If A is 
the momentum transferred by the proton, C,P+ = A_|_/2. The variable y is an integration 
variable in the gluon loop. 



J/\\f 



P' 



Figure bA: Elastic J/ip production (three further diagrams with the gluons connected 
in different ways have to be added). The two gluon lines and the incoming and outgoing 
proton are labelled by their respective fractions of the plus component of P = (P + P')/2. 



The lower part of the diagram in Fig. ^.4| is a generalization of the conventional gluon 
distribution (cf. Eqs. ( 4.26 ) and ( 4.31| )). It can be described by the non-forward gluon 
distribution 



1 



dx_e-*^^+^-/2(p'|Ft(0, x_, 0±yF{0, 0, Ox)/|P) . (5.16) 



110 1 that, as in the case of conventional parton distribu- 



Recently, it has been shown 
tions ||111|| , no time ordering of the operators in Eq. ( |5.16|) is required. 

The description of elastic meson production in terms of non-forward parton distri- 
butions is superior to the double- leading- log approach of |jl5|,0 since as corrections to 
the hard amplitude, meson wave function and parton distribution function can, at least 
in principle, be systematically calculated. However, the direct relation to the measured 
conventional gluon distribution is lost. A new non-perturbative quantity, the non-forward 
gluon distribution, is introduced, which has to be measured and the evolution of which 
has to be tested - a very complicated problem given the uncertainties of the experiment 
and of the meson wave functions involved. 

Over the recent years, the theory of non-forward parton distributions has developed 
into an active research field in its own rights, a detailed account of which is beyond the 
scope of this paper (see, however, [p.09|| for recent reviews). Important issues include the 



investigation of different models for non- forward distribution functions | |1 1 2| | , helicity- 
flip distributions [p.l3|| , the further study of non- forward evolution equations ||1 14|| , and 



63 



possibilities of predicting the non- forward from the forward distribution functions | 115 |. 
The latter suggestion relies on the observation that, at sufficiently high Q^, the non-trivial 
^ dependence (cf. Fig. ^]4| and Eq. (|5.16|) ) is largely determined by the evolution. 



Furthermore, following the basic results of |[15|, p!6|, a number of interesting phe- 
nomenological analyses of meson production have appeared. The analyses of ||116|| and 



117 | focus, among other issues, on the effects of the meson wave functions. More details 



of these approaches will be given in Sect. [TBI , when experimental results are discussed. 
The effects of Fermi motion and quark off-shellness have recently also been discussed 
Other interesting topics discussed in the literature include the form of the 

Hi: 



m 



energy dependence 
polarization 112 



shadowing effects [120|, Sudakov suppression [121 , effects of 



the t slope [ 117 , |123| , and the intrinsic transverse momentum of the 
hadron |[L 24|] . Vector meson production at large momentum transfer {\t\ ^ A^) is dis- 
cussed in |T7|,[T25|. 

Note especially the proposal of ||126|| to approach both transversely and longitudinally 
polarized diffractive p production on the basis of open qq production combined with the 
idea of parton-hadron duality. 

Vector meson electroproduction was also investigated in the framework of the model 
of the stochastic vacuum, which emphasizes the non-perturbative aspects of the process 
56|-p8|. The crucial observation is that, at realistic Q^, the asymptotic hard regime has 



not yet set in, and the typical transverse size of the qq fluctuation scattering off the 



target is not small. More details on this approach are found in Sect. |6.3| , as well as in 
Sects. Ol and fO- 



It should be emphasized that the above list of interesting subjects and related papers 
is in no way complete. 



5.2 Factorization 

Having discussed the leading order results for the cases of heavy vector meson production 
and light vector meson production with longitudinal polarization, the next logical step 
is to ask whether the systematic calculation of higher order corrections is feasible. For 
this, it is necessary to understand the factorization properties of the hard amplitude 
and the two non-perturbative objects involved, i.e., the meson wave function and the 
non-forward gluon distribution. In this section, a brief discussion of the general proof by 



Collins, Frankfurt and Strikman [22| is given. In addition, the essential role played by 



gauge invariance is explained in the framework of a simple model ||23|. This illustrates, in 



a particularly intuitive target rest frame approach, the physical mechanism underlying 
factorization properties in the small- a; limit. 



The analysis of is based on the method of leading regions discussed previously 
in Sect. [4.3| . The leading regions relevant for elastic meson production are shown in 
Fig. The main line of reasoning is analogous to the factorization proofs for inclusive 



hard scattering. As before, H is the hard subgraph, A is the subgraph with momenta 
coUinear with the incoming and outgoing proton, and S is the soft subgraph. Subgraph 
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B contains only lines which are coUinear with the produced meson. 




Figure 5.5: Leading regions in elastic meson production (cf. 



Factorization means that, to leading order in 1/Q, the amplitude corresponding to 
the process in Fig. |5.5| can be written as 



T 



dz / dyH{y,x/2)THiQ'^,y/x,z,n'^)(f)viz,iJ,'^ 



(5.17) 



where Th is the hard scattering amplitude, is the light-cone wave function of the vector 
meson produced, and H is the non-forward parton distribution of the proton. It could, 
for example, be the non-forward gluon distribution Hg of the last section. The variable 
X = XBj is the usual DIS Bjorken variable. Equation (|5.17|) implies that the soft subgraph 



in Fig. |5.5| is irrelevant, and only two parton lines, i.e., the minimal number required by 
the exchanged quantum numbers, connect the hard subgraph H with both A and B. 

According to p2[, the essential steps of the prove include the demonstration that 
the leading regions are indeed given by Fig. |5.5|, the derivation of a power counting 



formula showing that H is connected to A and B by the minimal number of lines, and 
the demonstration of factorization of the soft subgraph using Ward identities. 

An important complication compared to conventional factorization proofs results 
from the fact that final state interactions are important for both the outgoing proton 
and the produced meson. Thus, the soft gluons of S have to be factorized from both of 
these subgraphs at the same time. 

Furthermore, the so-called endpoint regions for the two quark lines connecting H and 
B have to be analyzed in detail. Endpoint regions are those kinematic domains where 
one of the two quarks carries almost all of the meson's longitudinal momentum. They 
are thus endpoints of the z integration in Eq. ( |5.17| ). An essential result of is the 
suppression of these endpoint regions to all orders of perturbation theory in the case of 
longitudinal photon polarization. 

In the case of transverse polarization, it was shown that the amplitude is suppressed 
by a power of Q relative to the longitudinal case. However, the endpoints are not sup- 
pressed relative to the intermediate z region, and therefore the factorization formula 
Eq. (^.17| ) can not be established in the transverse case. 

A discussion of helicity and transversity parton distributions, measurable via the 



polarization of the produced meson, is also contained in |22]. This will not be reproduced 



65 



here. More recently, factorization proofs similar to [22| have been given for the process 
of deeply virtual Compton scattering ||127|| , where the situation is simpler since no soft 



meson wave function has to be factorized. Note also that the factorization proofs of |22 



127|| , which use Breit frame kinematics, do not rely on the limit of small x. 



k+t 





Figure |5.6|: The leading amplitude for a point-like meson vertex. 



However, it is instructive to see in a particularly simple situation how factorization 
works specifically at small x, from the point of view of the target rest frame commonly 
used for the description of small-x processes ||2^. For this purpose, consider a very 
energetic scalar photon that scatters off a hadronic target producing a scalar meson 
built from two scalar quarks (see Fig. |5.6|) . The quarks are coupled to the photon and 
the meson by point-like scalar vertices ie and zA, where e and A have dimension of mass. 
The coupling of the gluons to the scalar quarks is given by —ig (r^ + r|^), where r and r' 
are the momenta of the directed quark lines, and g is the strong gauge coupling. 



Under quite general conditions 
f+ < g+, f_ < P_ and f 

effectively has the structure 



31 , the gluon momenta satisfy the relations 
~ Then the lower bubble in Fig. |5]| 



F^^^f., f , P) ~ (5(P_£+) F{t\) P^P" 



{5.U 



which is defined to include both gluon propagators and all colour factors. A similar 



expression was found by Cheng and Wu ||128|| in a tree model for the lower bubble. 

Assume that F restricts the gluon momentum to be soft, <^ Q"^. In the high-energy 
limit, it suffices to calculate 



M 



(2vr: 



_rpfj.v p ^ 



4(27r) 



where 



rpflU _j_ rpflU ^ rpfj,!/ 



is the sum of the upper parts of the diagrams in Fig. 



(5.19) 



(5.20) 



Note that, because of the symmetry of F^^, with respect to the two gluon lines, the 
amplitude T^'^ of Eq. (|5.20| ) is used instead of the properly-symmetrized upper amplitude 



1 



TZ^ii, q) = -[T^^(£, q) + T^^{-£', q)] . 



(5.21) 
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The two exchanged gluons together form a colour singlet and so the symmetrized ampli- 
tude T^^"^ satisfies the same Ward identity as for two photons, 

T,^l.{i,i\q)i,C = 0. (5.22) 

Writing this equation in light-cone components and setting i± = as appropriate for 
forward production, it follows that, for the relevant small values of £'_, i+ and 

Tsym,++ ~ ^± (5.23) 

in the limit i\ — >■ 0. Here the fact that the tensor T^"^, which is built from i', i and q, 
has no large minus components has been used. The i- integration makes this equation 
hold also for the original, unsymmetrized amplitude, 

a2 



~ £1 . (5.24) 

This is the crucial feature of the two-gluon amplitude that will simplify the calculation 
and lead to the factorizing result below. 



Consider first the contribution from diagram a) of Fig. |5.6| to the I- integral of 
which is required in Eq. ( ^.191) , 



Here N"^ = z{l — z)Q'^, z = k+/q+ and the condition = 0, enforced by the 5-function 
in Eq. ( p.l8| ), has been anticipated. 

Now / (i£_Tfe ++ and / (i£_Tc ++ each carry no i± dependence. So, to ensure the 
validity of Eq. (|5.24|) , the sum of the three diagrams must be 



/ 



ii^T,, = 4e/,, / - ^W^,^^ , (5.26) 



where 

r 1 1 1 

Note the behaviour obtained after a cancellation of l/Q"^ contributions from the 

individual diagrams. This cancellation, which is closely related to the well-known effect 



of colour transparency ||129|] , has been discussed in |^ in the framework of vector meson 
electroproduction. 

Introduce the k_[_ dependent light-cone wave function of the meson 

.2N_ ««+ 1,1,. ,,1. '■^ 



4.(.. kl) = -^J dk.ik, ^(*. (5.28) 

The final result following from Eqs. ( p.l9| ) and ( |5.26|) is a convolution of the production 
amplitude of two on-shell quarks and the light-cone wave function: 



^ = ( / 2|§F^^^^^^0 I'^'l '^'^^T^^^n^A^^ kl) . (5.29) 



41 -z) 

(N^ + kiy' 
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Figure |5.7| : Diagram for meson production with the vertex modelled by scalar particle 
exchange. 



This corresponds to the 0(£5_) term in the Taylor expansion of the contribution from 
Fig. pTB| a, given in Eq. ( ^.251 ). 

At leading order, factorization of the meson wave function was trivial since the point- 
like quark-quark-meson vertex V{k'^, {q' — k)"^) = iX was necessarily located to the right 
of the all other interactions. To see how factorization comes about in the simplest non- 
trivial situation, consider the vector meson vertex 



Vik\ {q' - kf) = J 



d^k' 



iXX 



12 



(27r)4 k'-^(q' -k')'^(k-k') 



(5.30) 



which corresponds to the triangle on the r.h. side of Fig. |5.71 . Here, the dashed line 
denotes a colourless scalar coupled to the scalar quarks with coupling strength A'. 

The diagram of Fig. [5.7| by itself gives no consistent description of meson production 
since it lacks gauge invariance. This problem is not cured by just adding the two diagrams 
|5.6| b) and c) with the blob replaced by the vertex V. It is necessary to include all the 
diagrams shown in Fig. 



The same gauge invariance arguments that lead to Eq. ( [5. 24] ) apply to the sum of 

Therefore, the complete result for T++, which is 



all the diagrams in Figs. |5.7| and [ST 
now defined by the sum of the upper parts of all these diagrams, can be obtained by 
extracting the term at leading order in the energy and Q^. Such a term, with a power 
behaviour ~ £\/Q^, is obtained from the diagram in Fig. (replace iX in Eq. ( ^.25|) 
with the vertex V of Eq. ( |5.30| )) by expanding around £_l = 0. It can be demonstrated 
that none of the other diagrams gives rise to such a leading-order l\ contribution (see |2^ 
for more details). 

The complete answer is given by the £^ term from the Taylor expansion of Eq. ( |5.25|) . 
The amplitude M is precisely the one of Eqs. ( [5. 29] ) and (5.28), with iX substituted by 
V of Eq. ( [5.301 ). The correctness of this simple factorizing result has also been checked 
by explicitly calculating all diagrams of Fig. ^.8| . 

The above simple model calculation can be summarized as follows. The complete re- 
sult contains leading contributions from diagrams that cannot be factorized into quark- 
pair production and meson formation. However, the answer to the calculation can be 
anticipated by looking only at one particular factorizing diagram. The reason for this 
simplification is gauge invariance. In the dominant region, where the transverse momen- 
tum £^ of the two t-channel gluons is small, gauge invariance requires the complete quark 
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Figure 5.8 : The remaining diagrams contributing to meson production within the above 
simple model for the meson wave function. 



part of the amplitude to be proportional to The leading £^ dependence comes exclu- 
sively from one diagram. Thus, the complete answer can be obtained from this particular 
diagram, which has the property of factorizing explicitly if the two quark lines are cut. 
The resulting amplitude can be written in a factorized form. 



5.3 Charm and high-^x jets 



In the two previous sections, the exclusive production of vector mesons was described 
as an example of a diffractive process with hard colour singlet exchange. As a differ- 
ent possibility of keeping the colour singlet exchange in diffraction hard, the diffractive 
production of heavy quarks p^-^ and of high-p^ jets [^8H 3T|] was considered by many 
authors. However, as will become clear from the discussion below, both processes can 
be associated with either soft or hard colour singlet exchange, and it is necessary to 
distinguish the two mechanisms carefully The semiclassical approach provides a 

very convenient framework for this analysis. 

One might expect the hard scale, provided by the transverse momentum of the 
jets, to ensure the applicability of perturbation theory. Indeed, the production of fi- 
nal states containing only two high-p^ jets can be described by perturbative two-gluon 



exchange [|13CI| , |131|| . This process has been studied in detail by several groups and higher- 
order corrections have already partially been considered ||28H3l|. 



Below, the two simplest configurations, qq and qqg, are discussed following PH . 
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In both cases, diffractive processes are obtained by projecting onto the colour singlet 
configuration of the final state partons. Although the discussion focusses on high-p^ 
jets, all qualitative results carry over to the case of diffractive charm production p7 



Technically, it is not important whether the hard scale in the diffractive final state is 
p]_ or ml. One has simply to replace the high-p_L qq jets with cc jets, whose transverse 
momentum will automatically be ~ m^.. However, there is a clear phenomenological 
difference. On the one hand, rric is fixed and not very large, while the hard scale can, 
at least in principle, be arbitrarily high. On the other hand, charm production is simpler 
to analyse since it does not require the identification of jets. 



Consider the production of a diffractive qq final state. Using the results of Sect. ^ 
for the transversely polarized photon one easily finds 



dt da dp 



/2 



t=0 



2N,{27tY 



-a 



n'^+pI 



(5.31) 



There are two essential differences compared to Eq. ( p. 201) : firstly, the colour trace is 
taken at the amplitude level to ensure colour neutrality of the qq state; secondly, two 
independent x±_ integrations are applied to the two factors W and W'^ . This is the result 
of Eq. ( |5.31|) being differential in t at t = 0, in contrast to Eq. ( |3.20| ), where the t 
integration has been performed. 

The cross section for large transverse momenta is calculated by expanding the inte- 
grand around p±_ = p'j_, as exercised in Sect. ^]2| for inclusive electropro duct ion in the 
longitudinal case. Note however that, due to the colour singlet condition, the first two 
terms (cf. Eqs. ( |3.22 ) and (p.23 )) do not contribute, so that the leading contribution 
comes from the third term, which is proportional to the second derivative of W at the 



origin (compare the discussion of diffraction in Sect. |3.2| ). Even higher terms of the Tay- 
lor series give rise to contributions suppressed by powers of p'l, thus demonstrating the 
dominance of the short distance behaviour of tTWx^{y±)- The leading order result reads 



dai 



dtdadp' 



12 



3847r2 



d 



p'± 



dp'^J N^+p'l 



(5.32) 



As the derivation illustrates, this cross section describes the interaction of a small qq 
pair with the proton. Hence, it is perturbative or hard. According to Eq. ( |5.6| ), the cross 
section Eq. (|5.32|) is proportional to the square of the gluon distribution. In order to 
obtain the t integrated cross section, one has to multiply Eq. ( p.32[ ) by the constant 



C 




A' 



(5.33) 



where A is a typical hadronic scale. The resulting cross section integrated down to the 
transverse momentum p'^ yields a contribution to the diffractive structure function 



which is suppressed by A'^/p'lcnt- 



As discussed in Sect. |^, a leading twist diffractive cross section for jets with p± ^ Q 
requires at least three partons in the final state, one of which has to have low transverse 
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momentum. It can be written as a convolution of ordinary partonic cross sections with 
diffractive parton distributions. In the case of high-px quark jets there is an additional 
wee gluon. The partonic process is then boson-gluon fusion, and the cross section 



d^dp 



/2 



dp 



/2 



(5.34) 



involves the diffractive gluon distribution of Eq. ( [4.541 ). 

In addition to boson-gluon fusion, the QCD Compton process can also produce high- 
P± jets. In this case either the quark or the antiquark is the wee parton. The corresponding 
cross section 



dai 



d^dp 



n 



■ ' — wi 



(5.35) 



involves the diffractive quark distribution of Eq. ( [4.53| ). An analogous relation holds in 
the antiquark case. 

The cross sections of Eqs. ( |5.34| ), ( p.35| ) for diffractive boson-gluon fusion and diffrac- 
tive Compton scattering can be evaluated along the lines described in |^^. In the leading- 
ln(l/a;) approximation, one obtains for the longitudinal and transverse boson-gluon fu- 
sion cross sections 



da J 



dadp'l 

daT 
dadp'l 

hA 



EgegOemas [a(l - a)]'^Q'^p'l 

(A^2 + p'_2)4 



27r3 



ln(l/x)/i^. 



y + (!■ 



-a] 



{p'l + 



(N^+p'l)^ 



ln(l/x)/i. 



A ■ 



(5.36) 
(5.37) 
(5.38) 



Similarly, one finds for the QCD-Compton cross sections 



daL 
dadp'l 

daT 
dadp'l 



277r3 [a(l - a)]Q^ 



277r3QV^ 



a(l — a) 



a{l — a) 



(5.39) 

(5.40) 
(5.41) 



where the constant /ijf is defined analogously to Eq. ( p.38| ) but with the U matrices in 
the fundamental representation. 

Comparing Eqs. ( p. 36] ), ( |5.37D with Eqs. ( p.39| ), ( |5.40| ), it is apparent that the con- 
figurations with a wee gluon are enhanced by ln(l/x) at small x relative to those with 
a wee quark or antiquark. The origin of this enhancement can be understood as follows. 
Eqs. ( |5.34| ) and ( ^.35| ) provide cross sections differential in ^, p'l and y or, equivalently, 
in ^, p'l and a, since y = x [Q^ + p'l/a{l — a)] /Q'^. The above p_L-spectra are obtained 
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after performing the ^-integration from y to some fixed <^ 1. The integral is domi- 
nated hy ^ ^ y, i.e., beta <C 1, where one has ydg/d^ ~ 1/,^. The ^-integration then 
yields a factor ln(l/x) + constant. In contrast, the diffractive quark distribution behaves 
as ydq/d^ ~ uH'^i and consequently the ^-integration only yields a constant. 

Simple arguments concerning colour, outlined in suggest an additional large 
suppression of the wee fermion contributions due to colour factors {h^ ^ hjr). This is 
in qualitative agreement with model calculations of Sects. ^TI] and |6.2| , which predict the 
diffractive gluon distribution to be much larger than the quark distribution. As a result, 
it can be claimed that configurations with a wee gluon dominate over those with a wee 
fermion in the small-x region relevant to diffraction. The latter shall be ignored from 
now on. 

The differential cross section for the leading order qq fluctuation can be calculated 
as described above (cf. Eq. ( ^.32|) ). The longitudinal and transverse cross sections are 



daL _ 2^,ela,WsA^G{0fC [a{l - a)fQ\a' - p'lf 



dadp'l 3 {m + p'ly 

dar 2J:^ela,^al7r'[^G{OfC {a^ + (1 - a)^)p'lN^ 



(5.42) 



(5.43) 



dadp'l 3 {m+p'D^ 

Identical differential distributions have been found for two-gluon exchange in leading 



order One can easily see that the region of a close to or 1 dominates, and that 
high-px configurations are unlikely. This statement will now be quantified. 

The quantitative differences between the qq and qqg configurations are particularly 
pronounced in the integrated cross section with a lower cut on the transverse momentum 
of the quarks. Since the overall normalization of the contributions is uncertain (it is 
inherently non-perturbative), the shape in p'l of each configuration is compared. Consider 
the quantity 

noo rl flu 

<Pl.cni)= /„ dp'l / ^, (5.44) 

h'l^c^^ JO dp'lda 

which is the fraction of events remaining above a certain minimum p'l. The integrand 
here is obtained by adding the contributions from longitudinal and transverse photons. 



Figure [5.9| shows the dependence of the corresponding event fraction on the lower limit, 
P'l cut- Each curve is normalized to its value at p'^ cut ~ ^ GeV^. One can see that the 
spectrum for the qqg configuration is much harder than that for the qq configuration. 
This is expected since in boson-gluon fusion p± is distributed logarithmically between 
the soft scale and Q thus resulting in a significant high-p^ tail above p'^ 

Let Mj be the invariant mass of the two-jet system in diffractive events containing two 
high-p^ jets in the diffractive final state. The measurement of this observable provides, 
in principle, a clean distinction between qq final states, where Mj = M^, and qqg final 
states, where Mj < M^. In practice, however, this requires the contribution of the wee 
gluon to the diffractive mass, which is responsible for the difference between and Mj, 
to be sufficiently large. To quantify the expectation within the semiclassical approach, 
consider the transverse photon contribution to the differential diffractive cross section 
da/dM'^dMj. 
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P'lcnt (GeV^) 

Figure |5.9| : The fraction of diffractive events witli p'^ above p'^ for of 10 GeV^ 
and 100 GeV^ (lower and upper curve in each pair). 



In the case of a qq final state this cross section can be obtained directly from 
Eq. (ra , 

"^""^ ^.ela^^aymOrC 5{M^ - M^) l^/fl^'^^ • (5.45) 



Here the 5-function setting = M| is only precise up to hadronization effects, which 
are expected to be of the order of the hadronic scale. The dependence on the transverse 
momentum cutoff enters via the variable k = 4p'^ ^^^^/Mj. 

By contrast, the mass distribution for diffractive processes with three-particle final 
states is not peaked at Mj = M^. Concentrating, as before, on the transverse photon 
polarization and on the contribution from the diffractive gluon distribution, the following 
formula can be derived from Eq. ( ^.34|) 

r4 



_ 2vrS e^a a y^^^iy^ Q + ^ r. Arctanhx/T^ - x/T 



K 



(5.46) 

Explicit results can be obtained by using model calculations for the diffractive gluon 
distribution (cf. Chapter P) or utilizing a simple parametrization (cf. the numerical pre- 
dictions of ^,|8l.) 

Note that the In (1/x) -enhancement present in Eq. ( p.37| ) can be recovered if the 
integration is performed in Eq. ( p.46| ). The origin of this enhancement is the integration 
measure (iM^/M^, which appears since y'^ {dg{y,i)/di) ~ l/M^ for sufficiently large. 
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The main contribution to the total cross section comes from the region where is 
significantly larger than M|. Thus it is clear, even without detailed calculations, that the 
distribution of Eq. (|5.46|) differs qualitatively from Eq. (|5.45|) . With sufficient statistics. 



a determination of the relative weight of soft colour singlet exchange, relevant in the qqg 
case, and hard colour singlet exchange, relevant in the qq case, should be feasible. 

The above calculations can be summarized as follows. The diffractive production 
of a qq final state with high or with charm quarks proceeds via hard colour sin- 
glet exchange. The description in the semiclassical framework reproduces the two-gluon 
exchange calculations. By contrast, high-p^ jets or charm in qqg final states are predom- 
inantly produced via boson-gluon fusion. The colour neutralization mechanism is soft, 
and the cross section is proportional to the diffractive gluon distribution. The boson- 
gluon fusion mechanism is distinguished by a much harder p± spectrum and a diffractive 
mass that is, on average, much larger than the invariant mass of the two-jet system. 

The energy dependence of the process could be very different for the above two mech- 
anisms. For hard colour singlet exchange, a steep rise is expected from the known small-x 
behaviour of the gluon distribution. For soft colour singlet exchange, the ^ dependence 
is expected to be less steep. 

In Ref. [^, the importance of higher order corrections to diffractive charm pro- 
duction was estimated in the framework of two-gluon exchange. A sizeable enhancement 
of the cross section was found. Clearly, the ccg final state is part of these corrections. 
However, the above discussion shows that this final state is dominated by the region 
where the gluon is soft. In this case, the ccg contribution is not perturbatively calculable 
and can not be considered an as correction to the hard cc process. Thus, the systematic 
calculability of corrections to hard cc or jet production is an interesting open problem, 
which may be related to the problem of defining 'exclusive' jet production at higher 
orders. 



5.4 Inclusive diffraction 

There are a number of attempts by several authors to approach the bulk of the diffractive 
DIS data, as described by the structure function F2 , from the perspective of two gluon 
exchange in the t channel. The degree to which perturbation theory is taken seriously 
varies significantly in the different investigations discussed below. 

In what can be possibly called the most modest approach, two gluons with an ap- 
propriate form-factor-like coupling to the proton are used as a simple model for colour 
singlet exchange, even if this exchange is believed to be non-perturbative. In the context 
of the diffractive electroproduction cross section at HERA, two gluon exchange calcula- 
tions were performed in [0], where the importance of the soft region was pointed out, and 
the result was parametrized in terms of an effective two gluon form factor of the proton. 

Furthermore, the two gluon calculation was used to describe both diffractive and 
inclusive cross sections in terms of the qq component of the light-cone wave function of 
the virtual photon and of cr{p), the cross section for a qq pair to interact with the hadronic 
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target 0. If the qq component dominates, which is, however, a non-trivial assumption, 
this approach allows one to link the diffractive cross section at t = to the inclusive DIS 
cross section via the optical theorem. 



In Ref. 29 , the diffractive structure function was discussed on the basis of the 



exchange of two gluons with non-perturbative propagators in the sense of the Landshoff- 



Nachtmann model [50|. Further analyses extend the two gluon exchange calculations to 



include the qqg component of the incoming photon pT| , |101| , P^ . 



Going one step further in the direction of perturbation theory, the known relation 



between the qq cross section a{p) and the inclusive gluon distribution ||132 



a(p) = ^«,N(x)K + 0(p^), (5.47) 

may be employed for the calculation of diffractive cross sections. This is similar to what 
was discussed in Sects. ^]T] and |5]^ in the case of meson production and high-p^ jet 
or charm production respectively. However, such a relation to the gluon distribution. 



employed, e.g., in |28, 133|, is problematic since the diffractive structure function is dom- 



inated by large transverse sizes of the qq fluctuation of the photon, where Eq. ( ^.47|) is 



not valid. It was emphasized in ||134|| that small, perturbative values of p become more 
important if the analysis is restricted to small diffractive masses. 

In spite of the evident problems with the applicability of perturbation theory to the 
diffractive structure function, it is still interesting to take the perturbative approach even 
further, calculating the BFKL leading-log corrections [|135|| to the two gluon exchange 
amplitude. A possible formal justification of such a treatment can be obtained by consid- 
ering diffractive DIS off a heavy quark-antiquark state, a so-called onium, which provides 
a perturbative scale in addition to the of the virtual photon. The focus is clearly on the 
energy dependence of the cross section, which so far can not be quantitatively described 
by non-perturbative QCD based methods. 

A detailed discussion of the BFKL technique of summing leading logarithms in the 
high-energy limit of perturbative QCD amplitudes is beyond the scope of the present 
review. The essential qualitative result, relevant for the following brief overview, concerns 
the scattering of two small colour dipoles at very high centre-of-mass energy ^i. It states 
that all corrections of the form as In s to the above process, which proceeds via two gluon 
exchange at leading order, can be summed. This results in a power-like growth of the 
amplitude with s. One may think of the energy logarithms as being associated with 
gluonic ladders, although the ladder topology does not exhaust all relevant diagrams. 

More recently, a colour dipole picture of the BFKL amplitude has been devel- 
oped p8| , |136H"l38|| . In this picture, each of the colliding colour dipoles radiates gluons, 
thus creating new colour dipoles with smaller energy, which are the source of further 
gluon radiation. Eventually, two dipoles, one from each of the two colliding cascades, in- 
teract via simple two gluon exchange. The equivalence with the original BFKL technique 
has been established for the most fundamental, but not for all relevant applications. In 
particular, the role of the large limit, which is used in addition to the leading loga- 



rithmic approximation ||136|,|137||, is not yet fully understood. 
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For illustration, the particularly simple formula for the total cross section of two onia 
with radius R and mass m is given (see, e.g., ||137|| ). It reads 



[s/m 



2\(ajp-l) 



y'(7/2)a,iVeC(3) \n{s/m) 



(5.4J 



where ajp = 1 + (4Ai'ca;s/7r) ln2 is the intercept of the BFKL pomeron. The reader is 
referred to the original papers ||135|| and to the modern introductory text [|17| for further 
details. 

The most straightforward application of the above perturbative techniques to diffrac- 
tion at HERA relies on postulating the exchange of a BFKL pomeron between target 
proton and partonic photon fluctuation. For forward diffraction, the conventional BFKL 
amplitude has to be introduced between, say, the qq fluctuation of the photon and a quark 
of the proton. More generally, the BFKL amplitude at non-zero t is needed. Correspond- 
ing formulae appear as a bye-product in the more general investigations of |^5|,^]. In 
the framework of the colour dipole approach, they are discussed in [ p.39| , |14CI|| . Clearly, 
method provides the desired rise of the cross section with energy or, e 



the method provides the desired rise of the cross section with energy or, equivalently, of 
2 with 1/^. However, with the simple asymptotic BFKL result, the rise appears to be 



far too strong (cf. Sect. |7[T| more details). 



It is a further, even more challenging theoretical problem to consider the double limit 
W"^ ^ ^ Q^, which corresponds to the simultaneous small-.^ and small-/? limit of 

■ The relevant process, illustrated in Fig. p.lO| , is similar to large-mass soft diffraction, 
discussed in Sect. W7^. Several authors proposed to use perturbative BFKL techniques for 



the calculation of this process of large-mass electroproduction. The obvious idea is to use 
the knowledge of the perturbative BFKL amplitude and to interpret the three pomerons 
in Fig. |5.10| as BFKL pomerons. In particular, the cut pomeron in the upper part of the 
diagram corresponds to gluonic radiation being responsible for the large diffractive mass 
created in the process. The dominance of gluons in the large mass region is well-known 
from fixed order perturbative calculations. 




Figure |5.1C1| : Triple pomeron vertex in diffractive electroproduction of large masses. 
In the framework of conventional BFKL summation, the triple pomeron vertex in 



diffractive electroproduction was considered in [^,|36|. It was found that the cut pomeron 
in Fig. is more complicated than the usual gluonic ladder. Four-gluon states were 
found to contribute to the leading amplitude. 
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The triple pomeron vertex was derived in the colour dipole approach and used for 
the calculation of large- mass diffractive electroproduction in p7t|141|| . As expected, an 
enhancement of the cross section both in the limit of large W"^ and large is obtained. 
More details of the phenomeno logical analysis will be given in Sect. |7.1| . 

In spite of the impressive theoretical work discussed above, the perturbative calcu- 
lability of the triple pomeron vertex appearing in large-mass hard diffraction remains 
questionable. On the one hand, it is not clear how far down the gluon ladder the in- 
fluence of the hard scale Q^, introduced by the 7*, extends. On the other hand, it is 
difficult to justify the ad-hoc introduction of a hard scale at the bottom of the diagram, 
as advocated in [^, since the proton is a soft hadronic object. 



Note finally that next-to-leading order results in the BFKL framework have recently 



become available (see ||142|| and refs. therein). For realistic values of as, the obtained 



corrections are larger than the leading order contributions and of opposite sign, thus 
complicating the theoretical status of the perturbative energy dependence even further. 
However, according to the very recent result of ||143|| , the use of BLM scale setting 
improves the situation dramatically. 
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6 Models for the Colour Field of the 
Proton 



In this chapter, three different models for the procedure of averaging over the target 
colour field configurations are described. Using the formulae of Sect. |4.4| , these models 
give rise to predictions for the diffractive quark and gluon distribution of the target 
hadron and thus for diffractive electroproduction cross sections. 



6.1 Small colour dipole 



A particularly simple model based entirely on perturbation theory has recently been 
suggested by Hautmann, Kunszt and Soper ||144|| . The authors study diffraction as quasi- 
elastic scattering off a special target photon that couples to only one flavour of very 
massive (M ^ A) quarks. The large quark mass justifies a completely perturbative 
treatment of the target and the diffractive system. In this situation, the required t channel 
colour singlet exchange is realized by two gluons coupling to the massive quark loop 
of the target. The analysis of ||144|| is based on the operator definitions of diffractive 



parton distributions (Eqs. (|4.29|) - (|4.31|) ), evaluated in an explicit two gluon exchange 
calculation. 



Here, the results of ||144|| will be derived in the semiclassical framework, following Ap- 
pendix B of In the semiclassical approach, the two t channel gluons are understood 
to be radiated by the massive quark loop and are treated as the colour field generating 
trVTtriy^ The semiclassical calculation proceeds as follows. 



Equations ( |4.53|) and (|4.54| ) have the structure 



D 



(6.1) 



where Fj (with i = q,g) is a linear functional depending on J tTWxj^{y±)tTW^^{y'^) 
interpreted as a function of y± and y'j_. To be differential in t, one simply writes 



D 



d^dt 



An 



triy^, trWl, e 



(6.2) 



with = —t. 

The field responsible for trW is created by a small colour dipole which, in turn, is 
created by the special photon that models the target. At leading order in perturbation 
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theory, the colour field of a static quark is analogous to its electrostatic Coulomb field. 
The field of a quark travelling on the light cone in direction at transverse position 0^ 
has therefore the following line integral along the x+ direction, 



2 



y4_ dXA 



(6.3) 



It is exactly this type of line integral that appears in the exponents of the non-Abelian 
phase factors U and f/^ that form W (cf. Eq. (|3.4|) ). A straightforward calculation shows 
that the function tiW produced by a dipole consisting of a quark at p±_ and an antiquark 
at I reads 



x(l 



2 



1 



(1 



-ik' pj_ 



(1 



(6.4) 



where Tp = 1/2 and Ta = Nc have to be used for the fundamental and adjoint represen- 
tation respectively. 

The final formulae for the diffractive parton distributions of the target are obtained 
after integrating over the transverse sizes of the colour dipoles with a weight given by 
the qq wave functions of the incoming and outgoing target photon. They read 



D 



dC, dt 



dzd p± I dz d p_|_ '^Fi 



Jx, Jx' ^-L 



(6.5) 



• \ H [V^7(^' P^^P'^^ 4) i^iiz, p^, 0±, e^)] [iJ*{z' p'^,p'^, e^) i:^{z', p'^, 0±, e±) 



where tYWxj^{y±) is produced by the field of a quark at p± and an antiquark at 0_l, and 
tiWx' {y'j_) is produced by the field of a quark at p'^ and an antiquark at 0_l, as detailed 
in Eq. (0). 

The wave function ip^yi^^ P±j^±j^±) characterizes the amplitude for the fiuctuation 
of the incoming target photon with polarization e and transverse momentum 0^ into a 
qq pair with momentum fractions z and 1 — z and transverse separation p^. Similarly, 
the wave function i/j*{z, p±,p'j_, e'j_) characterizes the amplitude for the recombination of 
this qq pair into a photon with polarization e' and transverse momentum p'j_ = —q±. The 
summation over the helicities of the intermediate quark states, which are conserved by 
the high-energy gluonic interaction, is implicit. 

The required product of photon wave functions can be calculated following the lines 



of [HH] and using the matrix elements of Appendix O. It reads explicitly 



^yiz, p±,p'^, e'^) ijj^{z, p±, 0_L, e_L) 



(6.6) 



2(27r)5 



/ tT^\z,k'^,M,e'^)^{z,k±,M,e±)e 

Jk I ,k', 



ip±(.k'j^-k±+zp'^) 
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where the notation of 11144 



$(z,fcx,M, ex) 



{kl + M2) 



[{1 — z) e± ■ a k_i_ ■ a — z k± ■ a e± ■ a + iM ej_ • a ] , (6.7) 



has been used, M is the quark mass, and ai^2 are the first two Pauh matrices. Note that 
for p'^ = 0, the average of the diagonal elements (e^ = e^) in Eq. ( |6.6| ) reproduces the 
well-known formula for the square of the photon wave function |^. 

Inserting Eq. ( |6.6|) into Eq. ( |6.5|) and introducing explicitly the required functionals 
Fi specified by Eqs. ( 4.53|) and ( 4.54 ), the formulae of | 144 | for diffractive quark and 
gluon distribution are exactly reproduced. For the lengthy final expressions the reader is 
referred to the original paper, where a number of plots, based on the numerical evalua- 
tion of these formulae, is also given. Qualitatively, the behaviour can be summarized as 
follows. The quark distribution (3{df^ / dt) falls off like /3 as /5 — > and approaches a 
constant value, which is small compared to intermediate /? points, as /5 ^ 1. The gluon 
distribution P{df^ / d^ dt) approaches a sizeable constant as /? — >• and a small constant 
as /? — *• 1. For t = 0, the behaviour at (3 —>■ 1 changes - quark and gluon distributions 
vanish approximately as (1 — /?) and (1 — /5)^. The overall normalization of the gluon 
distribution is found to be much larger than that of the quark distribution. 

Even though a real proton is very different from a small colour dipole, it would 
certainly be interesting to perform a phenomenological analysis on the basis of the above 
model. 



6.2 Large hadron 

In this section, the colour field averaging procedure is described for the case of a very 
large hadronic target, where a quantitative treatment becomes possible under minimal 
additional assumptions. The following discussion is based on [^, where the large 
hadron model was developed and applied to a combined analysis of both diffractive and 
inclusive structure functions (see also Sect. [4.4| and Appendix of this review). 



McLerran and Venugopalan observed that the large size of a hadronic target, realized, 
e.g., in an extremely heavy nucleus, introduces a new hard scale into the process of 
DIS |^0|]. From the target rest frame point of view, this means that the typical transverse 



size of partonic fluctuations of the virtual photon remains small [^, thus justifying the 
perturbative treatment of the photon wave function in the semiclassical calculation. 

The basic arguments underlying this important result are best explained in the simple 
case of a longitudinally polarized photon coupled to scalar quarks with one unit of electric 
charge. As far as the Q^-behaviour of the total j*p cross section is concerned, this is 
analogous to the standard partonic process where a transverse photon couples to spinor 
quarks [|145 . 



In analogy to |^, the longitudinal cross section can be written as 

aL= f d^pMp)WL{p), (6. 
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with the square of the wave function of the virtual photon given by 

3 Op, 



Wl{p) 



daN'K^{Np). 



(6.9) 



Here p = \p±\ is the transverse size of the qq pair, a is the longitudinal momentum 
fraction of the photon carried by the quark, A^^ = a{l — a)Q'^, and Kq is a modified 
Bessel function. Note that, in contrast to 0, Wl is defined to include the integration 
over a. 



Within the semiclassical approach, the dipole cross section a{p) is given by 



(fx±tT 1 — U{x±)U''' {x± + p±) 



(6.10) 



but it is convenient to formulate the following arguments in terms of the more general 
quantity cr{p). 

The functional form of a{p) is shown qualitatively in Fig. |6.1| . For conventional 
hadrons of size ~ 1/A (where A ~ Aq^j^), its typical features are the quadratic rise 
at small p {cr{p) ~ with a proportionality constant 0{1)) and the saturation at 
cr(p) ~ 1/A^, which occurs at p ~ 1/A. Consider now the idealized case of a very large 
target of size 77/ A with 77 ^ 1 (77 ~ A^^^ for a nucleus). It is easy to see that at small 
p the functional behaviour is given by a{p) ~ rj^p^ while saturation has to occur at 
cr(p) ~ Iff / for geometrical reasons. It follows that the change from quadratic rise to 
constant behaviour takes place at p ~ 1/y^A, i.e., at smaller p than for conventional 
targets. 




large hadron 
(radius --q/A) 



conventional 
hadron 
(radius ~1/A) 



> P 



'1/A 



Figure |6.1| : Qualitative behaviour of the function cr(p). 



From the above behaviour of (t(p), the dominance of small transverse distances in 
the convolution integral of Eq. ( |6.8D will now be derived. For this purpose, a better 
understanding of the function Wl{p) is necessary. Recalling that Kq{x) ~ ln(l/a;) for 
X ^ 1 while being exponentially suppressed for x 1, it is easy to see that Wl{p) ~ 
Q"^ \t? {1 / p^ Q"^) for p ^ l/Q and Wl{p) ~ 1/p'^Q^ for p > l/Q. Here numerical constants 
and non-leading terms have been suppressed. 

Under the assumption A^ <^ r/A^ <^ Q^, the integral in Eq. ( |6.8| ) can now be esti- 
mated by decomposing it into three regions with qualitatively different behaviour of the 
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functions Wl{p) and cr(p), 

^L = crl + al + af={ I ' + I ' - +1 J vrt/pV(p)W^z.(p) • (6.11) 




Of the three contributions 

~ C'"^'^''^'^ - $HQ'W') (6.12) 



7l/„A2 A2 



the second one dominates, giving the total cross section 



~ ^ In(QV^A^) • (6.13) 



It is crucial that the third integral is dominated by contributions from its lower limit. 
Therefore, the overall result is not sensitive to values of p that are larger than l/^/rjA. 
Phrased differently, for sufficiently large targets the transverse size of the qq component 
of the photon wave function stays perturbative. 

This result can be carried over to the realistic case of a transverse photon and spinor 
quarks, which is obtained simply by substituting K^^pN) with 2[a^ + {1 — ay'\K\{pN) 
in Eq. ( |6.9|) . From the above derivation, one can also expect the result to hold for the 
transverse size of the qqg component of the photon wave function, both in the case of 
inclusive and diffractive scattering. 

Note that this does not imply a complete reduction to perturbation theory since the 
long distance which the partonic fluctuation travels in the target compensates for its 
small transverse size, thus requiring the eikonalization of gluon exchange. 

Within this framework, it is natural to introduce the additional assumption that the 
gluonic flelds encountered by the partonic probe in distant regions of the target are not 
correlated (cf. and the somewhat simplifled discussion in ||61|). Thus, one arrives at 
the situation depicted in Fig. |6.2| , where a colour dipole passes a large number of regions, 
each one of size ~ 1/A, with mutually uncorrelated colour flelds Ai ... An. 

Consider the fundamental quantity Wi.j^(?/_L)jj[A] Wl^{ii'^ki[A\ which, after specify- 
ing the required representation and appropriately contracting the colour indices ijfc/, 
enters the formulae for inclusive and diffractive parton distributions (cf. Sect. |4.4| and 
Appendix 0). According to Eq. (|3.17 ), this quantity is the sum of four terms, the most 



complicated of which involves four U matrices, 

[U^AA] Ul^+yAA]}^^ [u.^+y'^A] Ul^A]}^^ (6.14) 

= {U,, K] [A,] Ul^^y^ [A,] ■ ■ ■ Ul^+y^ [An] } . . 

X {[/.,+,^ [An] ■ ■ ■ U^^+y.^ [A,] Ul [A,] ..-Ul [An] 



82 



^1 ^2 ^3 



Figure 6.2: Colour dipole travelling through a large hadronic target. 



The crucial assumption that the fields in regions 1 ... n are uncorrelated is implemented 
by writing the integral over all field configurations as 

(6.15) 

product of independent integrals. Here the appropriate weighting provided by 
the target wave functional is implicit in the symbol £4. 

Under the integration specified by Eq. (|6.15| ), the V matrices on the r.h. side of 
Eq. ( |6.14D can be rearranged to give the result 



{f/.,[A] Ul^^yAA^^^ {u.,^y',\A\ f/tjA]}^^ (6.16) 

J A\ J An J 



kl 

To see this, observe that the Ai integration acts on the integrand {U^^ [Ai] Ul_^^y^ [Ai]}iij/ 
{Ux^+y'AAi]UlAAi]}j:'ii transforming it into an invariant colour tensor with the in- 
dices i'i'k'V. The neighbouring matrices f/j;^[y42] and f/j!^[y42] can now be commuted 
through this tensor structure in such a way that the expression {U^^ [^2] Ul^j^y^ [A^WiUjii 
{Ux^+y'J[A2\ UlJ[A-2\}k"i" emerges. Subsequently, the A2 integration transforms this ex- 
pression into an invariant tensor with indices i"j"k"l". Repeating this argument, one 
eventually arrives at the structure displayed on the r.h. side of Eq. (|6.16|) . 



To evaluate Eq. ( |6.16D further, observe that it represents a contraction of n identical 
tensors 

Fijki = [ {U.A^m] f/J^+^JAJ},, {f/,^+,^[AJ UlJA^]}ki , (6.17) 

where the index m refers to any one of the regions 1 ... n into which the target is 
subdivided. At this point, the smallness of the transverse separations y± and y'j_, enforced 
by the large size of the target, is used. In fact, for a target of geometrical size ~ n/A 
(where n ^ 1), the relevant transverse distances are bounded by y"^ ~ y'^ ~ 1/raA^. 

Assuming that size and x± dependence of typical field configurations Am are charac- 
terized by the scale A, it follows that the products UxJJl^+y^ and Ux^+y'JJl^ are close 
to unit matrices for all relevant y± and y'^. Therefore, it is justified to write 

UxAAm] Ul^+yA^A = exp{2r'^r(x±,i/^)[AJ} , (6.18) 
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where are the conventional group generators and are functions of x± and y± and 
functional of A^- Equation ( |6.18|) and its analogue are expanded around y± = y'^ = 
(which corresponds to f"'{x±, 0) = 0) and inserted into Eq. (|6.17|) . At leading non-trivial 
order, the result reads 



ijkl 



1 - ^iCniy' + y": 



+ i{y±y'^mT^i 



where Cr is the Casimir number of the relevant representation {Cn 
constant 7 is defined by 



/ nx±,y^)f{x^,y'^) = iS'^'iyW^) + 0{y^") 

J A 



(6.19) 

Cf^a) and the 

(6.20) 



Note that the absence of terms linear in f°- and the simple structure on the r.h. side of 
Eq. ( |6.20|) are enforced by colour covariance and transverse space covariance. The absence 
of an explicit x± dependence is a consequence of the homogeneity that is assumed to 
hold over the large transverse size of the target. Neglecting boundary effects, the x± 
integration is accounted for by multiplying the final result with a parameter Q ~ ri^ / h? 
that characterizes the geometrical cross section of the target. 



Substituting the n tensors Fiju on the r.h. side of Eq. ( |6.16D by the expression given 
in Eq. ( |6.19| ) and contracting the colour indices as appropriate for the inclusive and 
diffractive case respectively, one obtains, in the large- A^^c limit. 



1 - ^lCR{y± - y'^)^ 



(6.21) 



4 



^iCniyl 



(6.22) 



where dpi is the dimension of the representation. 

Since n is assumed to be large and the typical values of y"^ and y'"^ do not exceed 
l/nA^, the formula (1— x/n)" ~ exp[— x] can be applied to the r.h. sides of Eqs. (|6.21| ) and 
( |6.22| ). Furthermore, contributions proportional to {UxJJl^+y^}ij5ku ^ij{Ux^+y'JJl^}ki 
and 5ij5ki have to be added to obtain the complete expression for Wx^{yA_)ij {y'±)ki- 
The corresponding calculations are straightforward and the result reads 



X I J A 



tr 



Qdf 



l-e-^^y -e-^^y' + e-"«(s^^-2^^) ,(6.23) 



X I J A 



-o-rV 



(6.24) 



where an = n'^CR/2 plays the role of a saturation scale. 



The above calculation, performed at large and for the case of a large target 
subdivided into many uncorrelated regions, has no immediate application to realistic 
experiments. However, it provides a set of non-perturbative inclusive and diffractive 
parton distributions which are highly constrained with respect to each other. For the 
purpose of a phenomenological analysis, it is convenient to consider Vt and a = n'jNc/4: as 
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new fundamental parameters, giving rise to the following formulae for the basic hadronic 
quantities required in Sect. and Appendix 0, 



1 

1 



tTW^(y^)trW^Uy'^ 



= nN^ 


1 - 


e-^y' - 


= 


1 - 




= 


"l - 





-ay 



,(6.25) 
(6.26) 
(6.27) 



Here the indices and A stand for the fundamental and adjoint representation. A similar, 
Glauber type y'^ dependence has been recently used in the analyses of [p.46| , |147|| . Note 
that according to Eqs. ( |6.25| )-( |5.27| ) the diffractive structure function is not suppressed 
by a colour factor relative to the inclusive structure function, as originally suggested 
in H. 

Using Eq. ( |6.26| ) in the generic semiclassical formula Eq. ( [4.53| ) for the diffractive 
quark distribution, the explicit result 



27r3^2 

is obtained. Here hq{(3) is an integral over two Feynman-type parameters. 



(6.28) 



X 



+ x) Y 



(x + x')v^+(l-/3) 



(6.29) 



X 



+ 



X 



V^+x V^+x', 



which can be evaluated analytically at /5 = and (3=1 yielding /ig(0) = 1/2 and 
hq{l) = SvrVS - 2. 

Analogously, Eqs. (|6.27|) and ( [4.54| ) give rise to an explicit formula for the diffractive 
gluon distribution. It reads 



d^ 27r^e/3 
where hg{(3) is given by the two-dimensional integral 

( 1-(3 + 3{1 + xYP 

roo I n -J- 



MP) 



(6.30) 



doc ditXj 



Mi + xYii- p + {i + xYP) 



I- 13 + ?,{! + x'f (3 
{l + x')'\l-P+{l + x'fl3f 



{x + x')P+{l- 13) 



X 



rr^ + 1 + x' 



X 



(6.31) 

This integral is easily evaluated for /5 = and (3 = 1 yielding hg{Q) = 41n2 and hg{l) = 
457r2/32 — 17/2. For general (3, hq{(3) and hg{(3) can be evaluated numerically. The results 
can be inferred from the solid curves in Fig. |6.3| , where the distribution dH/dC, = Qdq/d^ 
is displayed to account for the 3 generations of light quarks and antiquarks. The total 
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Figure |6.3| : Diffractive quark and gluon distributions in the large hadron model at the 
initial scale Qo and after evolution. 



normalization, the value of ^, and the evolution given in Fig. |6.3| are not relevant for 
the present section and will be discussed in the context of the phenomenological analysis 
of Sect. |71|. 

The diffractive distributions displayed in Fig. |6.3| are multiplied by (5 and thus reflect 
the distribution of momentum carried by the partons. The quark distribution is peaked 
around ^ 0.65, thus being harder than the distribution — /5) suggested in |T0[. It 
vanishes like /3 for /3 — and like (1 — /3) at large /3; the gluon distribution f3dg/d^, on 
the other hand, approaches a constant for /3 — ^ and falls off like (1 — /?)^ at large (3. This 
asymptotic behaviour in the small- and large-/? region is in agreement with the results 
obtained in the perturbative approach of ||144|| at t = 0. In spite of the (1-/5)^ behaviour, 
gluons remain important even at large (3, simply due to the large total normalization of 
this distribution. 

Very recently, a closely related discussion of diffractive and inclusive structure func- 
tions has been given in [|148|| . The authors focus on the process where the target hadron 
remains intact, discussed in the first part of Sect. |3.4| , and its relation to inclusive DIS. 
Technically, the results obtained are very similar to those of the large hadron model as 
described above It is also emphasized in ||148 | that the results can be generalized to 
the case of conventional hadrons if the assumption of a Gaussian distribution of colour 



sources of [oUI is correct. 



6.3 Stochastic vacuum 



A further fundamentally non-perturbative approach to high-energy hadronic processes, 
which has recently been applied to diffractive structure functions by Ramirez is 
based on the model of the stochastic vacuum. The model was originally developed by 
Dosch and Simonov in Euclidean field theory ||5^. A detailed description of the stochas- 
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tic vacuum approach to high-energy scattering, introduced originally in |53| and closely 



related to the eikonal approach of [|^, can be found in ||149|] (see also ||150|| for a com- 



prehensive review). The method was first applied to diffractive DIS in [Q, where vector 
meson production processes were considered. Here, only a brief description of the main 
underlying ideas will be given. 



The fundamental assumption underlying the model of the stochastic vacuum of |K 
is that of a convergent cumulant expansion for the vacuum expectation value of path 
ordered products of field operators. To calculate the average (...) of the path ordered 
exponential 



P exp { 0{s)ds] = Y] dsi ds2... dsnO{si)...0{sn) , (6.32) 

\Jo / „_n "'0 JO JO 



E 

n=0 

the so-called path ordered cumulants ((••.)), defined by 

(1) = ((1)) 
(1,2) = ((1,2)) + ((1))((2)) 

(1,2,3) = ((1,2,3)) + ((1))((2,3)) + ((1,2))((3)) + ((1))((2))((3)) 

(1,2,3,4) = ((1,2,3,4))+ ■■■ +((1))((2))((3))((4)) 



(6.33) 



are introduced. An expansion in the above cumulants can be applied to the Wegner- 
Wilson loop in a non-Abelian gauge theory. The supposition that the cumulants are 
decreasing sufficiently fast with increasing distance between the operators leads to the 
area law in the purely gluonic case. Neglecting all cumulants higher than quadratic in the 
fields amounts to the assumption of a Gaussian process, where all higher correlators can 
be obtained from the two-point Green's function. All of the above is naturally formulated 
in Euclidean space. The two-point correlator or, more precisely, its analytic continuation 
to Minkowski space, is the fundamental object in applications to high-energy scattering. 

In the treatment of diffractive DIS, the stochastic vacuum approach describes the 
hadron as well as the virtual photon in terms of fast partons moving in opposite direc- 
tions. The colour field facilitating the interaction is soft in the centre-of-mass frame of 
the 7*p collision (cf. the general discussion of soft hadronic processes in |^3[). In this 
situation, the partons from both sides interact with the field in an eikonalized way, and 
the actual model of the stochastic vacuum is used to evaluate the correlation function of 
the resulting oppositely directed light-like Wegner- Wilson lines. 

In the simplest case, where the photon fiuctuates into a qq pair and the hadron is 
modelled as a quark-diquark system, the interaction amphtude of two Wegner- Wilson 
loops has to be calculated (see Fig. p?^ ). Introducing the notation 



W = P exp 



-ig j A^{z)dzf' 

JdS 



(6.34) 



for the path ordered integral around the loop dS, which is extended to light-like infinity 
in both directions, this amplitude can be written as 

J(6,i?i,4) = (^-^tT{W,{b/2,R^) - 1) ^tr(W^2(-V2,i?2) - 1))^ • (6.35) 
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Figure |6.4| : Wegner- Wilson loops formed by the paths of quarks and antiquarks inside 
two dipoles. The impact parameter b connects the centres of the two loops, while Ri and 
R2 point from the quark to the antiquark line of each dipole. All three are vectors in the 
transverse plane of the collision (figure from [p3 ) . 



The notation of this section follows the original papers reviewed (see, e.g., 



which means that the function defined in Eq. (|6.34| ) differs by a unit matrix from the 
closely related function W used previously in the context of the semiclassical approach. 
The brackets {■■■)a denote the vacuum expectation value or, in the functional language, 
the integration over all colour field configurations. 

The line integrals in Eq. ( |6.35| ) are transformed into surface integrals with the help 
of the non-Abelian Stokes theorem 



Pexp 



JdS 



5exp 



S 



(6.36) 



Here F^y{z,uj) are the field strength tensors F^u{z), parallel-transported to a common 
reference point u. The operator P5 denotes an appropriate surface ordering of these 
matrix valued tensors (see ||149|| and refs. therein for further details). The surface S that 
is chosen for each of the two Wegner- Wilson loops is the upper side of the pyramid with 
the loop as base and the origin of the co-ordinate system as apex (see Fig. |6.4D . The 
reference point u is also chosen to be at the origin. 

The interaction amplitude for two Wegner- Wilson loops, in its explicit form of 
Eq. (|6.35|) or expressed through surface integrals according to Eq. ( |6.36| ), can not be 
evaluated by direct application of the cumulant expansion methods discussed above. 
The problem is that those methods are adopted for one path ordered integral, in con- 
trast to the two path (or surface) ordered integrals required for J(6, Ri, R2). In Ref. fS^ 
it was suggested to apply the Gaussian factorization hypothesis directly to products of 
the fields F^^{z,uj), 



(F(l)...F(2n)) 



Y: {F{i,)F{j{))...{F%)F{j^)) 



(6.37) 



all pairings 
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where the single argument of F stands for Lorentz and colour indices and space-time 
co-ordinates. Note, however, that this factorization assumption is not equivalent to the 
original approach of (cf. the discussion in ||15CI| , |151|| ). 



Equation (|6.35|) can now be evaluated if the fundamental correlator 

g'F^^,ix,u;)F^^{y,uj))^ (6.38) 

is known. Under the further assumption that this correlator depends neither on the path 
used for transporting the field to uj nor on the position of the reference point lj itself, 
the most general form reads 

9'^F^uix,uj)F;l^{y,uj))^ = ^^^^^ _ _^^ ^K{g^pg^^ - gp^g^p)D{z^ /a^) (6.39) 

+ (1-k)^ [dp{zpgua-Zagup) + du{z„g pp- Zpg p„)] Di{z^/a^) 

where D and Di are, a priori, two independent functions, and (g'^FF) is the gluon 
condensate. Note that the model of the stochastic vacuum is formulated in Euclidean field 
theory and the analytic continuation to Minkowski space is non-trivial. The intricacies 
of this process and, in particular, the constraints it imposes on the shape of the functions 
D and Di will not be discussed here. For completeness, the correlation functions used, 
e.g., in are given: 

^ ^ ^ V / ; ^ j (27r)4 (p - (37r/8)2)4 ^ ' 

Detailed information about the parameters entering Eq. ( |6.39|) and the functional form of 
D and Di is available from low energy hadronic physics as well as from lattice calculations 



(see, e.g., the recent results of ||152|| ). Given these low-energy parameters, a large number 
of soft hadronic high-energy scattering processes is successfully described by the model 
of the stochastic vacuum. For the present review, it is sufficient to state that, within the 
present model, the behaviour of the correlator in Eq. ( |6.39| ) is quantitatively known. 

Thus, the method for evaluating the fundamental dipole amplitude of Eq. ( |6.35| ) 
is now established. The amplitude vanishes for small dipoles, as expected from colour 
transparency, and grows linearly for large dipoles, as suggested by the geometric picture 
of string-string scattering. 

The calculation of diffractive processes is straightforward as soon as a specific quark- 
diquark wave function of the hadronic target is chosen. The procedure of folding the 
dipole amplitude with the virtual photon and the hadron wave functions introduces a 
certain model dependence on the hadronic side. 

An essential difference between the present framework and the semiclassical approach 
discussed above is the frame of reference where the model for the hadronic target is 
formulated. From the perspective of the virtual photon, the situations are similar: the 
qq fluctuation scatters off a soft colour fleld. However, in the semiclassical approach (cf., 
e.g.. Sect. |6^ ), this field is modelled on the basis of ideas about hadron colour fields in the 
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rest frame of the hadron. By contrast, the stochastic vacuum model uses an intermediate 
frame, where both projectile and target are fast (e.g., the centre-of-mass frame), as the 
natural frame for the colour field. The target hadron is described in terms of quarks 
interacting with this field in an eikonalized way, and the field, as seen by the projectile, 
is a result of the distribution of these partons and of the gluon dynamics in the stochastic 
vacuum model. 

As an interesting consequence of the above discussion, the conditions required for 
the eikonal approximation to work are different in the two approaches. While, in the 
semiclassical model, the partons from the photon wave function have to be fast in the 
target rest frame, the stochastic vacuum approach requires them to be fast in, say, 
the centre-of-mass frame. The latter is a far more stringent condition, which can be 
numerically important in phenomenological applications [p^,|59[ . 



The reader is referred to Sections ^^1] and for a brief discussion of recent applica- 
tions of the presented model to diffractive electroproduction data. 
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7 Recent Experimental Results 



In this chapter, a brief discussion of experimental resuhs in diffractive electroproduction, 
focussing in particular on the diffractive structure function, on specific features of the 
diffractively produced final state, and on exclusive meson production is given. The discus- 
sion is aimed at the demonstration of interesting features of different theoretical models 
and calculational approaches in the light of the available data (see also ||153 | for recent 
reviews). From an experimentalist's perspective, what follows can only be considered a 
very brief and incomplete overview. 



7.1 Diffractive structure function 

As already explained in Chapter ^ the observation that diffractive DIS is of leading 
twist and the ensuing measurement of diffractive structure functions lie at the heart of 
increased recent interest in the field. While the basic experimental facts were stated at the 
beginning of this review, the present section supplies some additional details concerning 
F,p and compares observations with theoretical ideas. 

To begin with, a brief discussion of the most recent and most precise data on F^, 



produced by the HI and ZEUS collaborations |173| , [7^ , is given. As mentioned previously, 
the most precise analyses have to be based on measurements where the scattered proton 
is not tagged. In Ref. this problem is handled by starting from a cross section 



dxdQ^d^dtdM^' ^ ■ ^ 

where, by definition, the two clusters X and Y are separated by the largest rapidity gap 
in the hadronic final state of the event. This definition allows for the possibility that 
the proton breaks up into a final state Y with mass square P'^ = My (cf. Fig. 
A structure function F^^^^ is then defined, as in Sect. |2.3| , after integration over t and 
My- The event selection of is based on the requirement that X is fully contained 
in the main detector while Y passes unobserved into the beam pipe. This implies that 
My $ 1.6 GeV and |t|^lGeV^, thus approximately specifying the relevant integration 
region for My and t. A rapidity cut ensures that X is well separated from Y. 

The most recent ZEUS analysis is based on the so-called Mx method, introduced 



in [|7T|. Similarly to what was discussed above, it is assumed that X is contained in the 
detector while Y escapes down the beam pipe. Thus, Mx is simply defined to be the full 
hadronic mass contained in the main detector components. In contrast to the HI analysis, 
the diffractive cross section, integrated over all t and My, is defined by subtracting from 
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this full cross section the non-diffractive contribution. The subtraction is based on the 
ansatz 



dM 



d In Ml- 



D + cexp(6 1nMj.) (7.2) 



for the event distribution in the region of not too large masses. Here D is the diffractive 
contribution and the second term is the non-diffractive contribution. The exponential 
In Mx dependence is expected as a result of the Poisson distribution of particles emitted 
between current and target jet regions in non-diffractive DIS. This leads to an exponential 
distribution of rapidity gaps between the detector limit and the most forward detected 
particles and therefore to the above Mx dependence. The diffractive contribution is 
not taken from the fit result for D but is determined by subtracting from the observed 
number of events the non-diffractive contribution that corresponds to the fit values of h 
and c. 

For more details concerning the data analysis and the experimental results themselves 
the reader is referred to the original papers |]73|,^] and to the figures in the remainder 
of this section, where theoretical models are compared with measured values of ■ It 
can be seen from the direct comparison of HI and ZEUS data found in [|7^] that the two 
methods give similar results, although the HI values for F2^^^ have a tendency towards 
a faster rise for any given (3. 

The Mx method of ZEUS has the advantage of subtracting events which happen 
to have a rapidity gap in spite of non-singlet colour exchange. It also allows for diffrac- 
tive events with relatively forward particles, which are excluded by the rapidity cut of 



HI ||154|| . However, the subtraction of the non-diffractive background also introduces new 
uncertainties, in particular the dependence on the region in Mx where the fit according 
to Eq. ( |7.2| ) is performed. 

Note also that a measurement of F^^^\ based on the use of the leading proton spec- 



trometer of ZEUS, was reported in |73]. The results for F:[^^^\ obtained by explicit t inte- 



gration, are probably the cleanest ones from both the theoretical and experimental per- 
spectives. They agree with the latest F^^^^ measurements discussed above, but at the mo- 
ment they are, unfortunately, less precise because of limited statistics. The observed t de- 
pendence of F^^^^ can be parametrized as e^*, with b = 7.2 ± l.l(stat.) j:Q;g (syst.) GeV~^. 



Leading twist analysis 

As mentioned previously, the leading twist behaviour of small-x diffraction is one of its 
most striking features. For this reason, and in keeping with the general perspective of this 
review, it is convenient to start the discussion with calculations, such as the semiclassical 
approach, that focus on the leading twist nature of the process. 

Detailed numerical predictions exhibiting this behaviour were first made in based 
on simple aligned jet model calculations with soft two gluon exchange. Further theoret- 
ical and numerical developments by these and other authors [p^ , |131| , p.39|| , focussing, in 



particular, on high-mass diffraction and higher twist effects, are discussed below. 

Data analyses based on the idea of parton distributions of the pomeron p| and their 
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Figure |7.1| : The structure function F'^^^\^, ^jQ"^) computed in the semiclassical ap- 
proacli with HI data from [^. The open data points correspond to < 4 GeV^ and 
are not included in the fit. 



evolution were performed early on by many authors [ 155 , FSS]. The most recent HI 
publication includes such a partonic analysis as well. 



Here, a detailed description of the approach of will be given, which is based on 
diffractive parton distributions in the sense of Berera and Soper calculated in the 
semiclassical approach with a specific model for the averaging over all target colour fields 
(cf. Sect. p.'2\ ). As can be seen from Figs. |7[T] to [7751 , a satisfactory description of F^^^^ is 
achieved in this framework. 
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Figure |7.2| : The structure function F:^^^\^, f^^Q"^) computed in the semiclassical ap- 
proach with ZEUS data from [^. The open data points correspond to < 4 GeV^ 
and are not included in the fit. 



In the presented approach, the diffractive distributions of Sect. and the inclusive 
distributions of Appendix are used as non-perturbative input at some small scale Q^. 
They are evolved to higher using the leading-order DGLAP equations ||157||. The non- 



perturbative parameters of the model as well as the scale Qq are then determined from a 
combined analysis of experimental data on inclusive and diffractive structure functions. 

At first sight, the semiclassical description of parton distribution functions always 
predicts an energy dependence corresponding to a classical bremsstrahlung spectrum: 
q{x),g{x) ~ 1/x. However, one expects, in a more complete treatment, a non-trivial 
energy dependence to be induced since the field averaging procedure encompasses more 
and more modes of the proton field with increasing energy of the probe (cf. the discussion 
at the end of Sect. p.4|) . This energy dependence is parametrized in the form of a soft, 
logarithmic growth of the normalization of diffractive and inclusive parton distributions 
with the collision energy ~ consistent with the unitarity bound. As a result, the 
additional parameter L is introduced into the formulae of both Sect. |6.2| and Appendix 0, 
where the overall normalization factor Q has to be replaced according to 

n^n{L-\nxf . (7.3) 
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Figure |7.3| : Dependence of the diffractive structure function F^*-^^ on [5 and Q^, com- 
pared to data from HI (left) and ZEUS (right) [^,|74|. Open data points correspond to 
< 4 GeV^. The charm content of the structure function is indicated as a dashed hne. 



The following expressions for the diffractive parton distributions are obtained, 



(7.4) 
(7.5) 



where the functions hq^g{f3) are defined in Sect. |6.2| . Corresponding expressions for the 
inclusive input distributions of quarks and gluons are given in Appendix 0. 

Thus, the input distributions depend on a, fi, L, and on the scale Qq at which these 
distributions are used as a boundary condition for the leading-order DGLAP evolution^. 

^Note that the two variables a and Vl can not be combined into one since the inclusive quark distri- 
bution depends on them in a more complicated way. 
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At this order, the measured structure function F2 coincides with the transverse structure 
function. In defining structure functions and parton distributions, all three light quark 
flavours are assumed to yield the same contribution, such that the singlet quark distribu- 
tion is simply six times the quark distribution defined above. Valence quarks are absent 
in the semiclassical approach. Charm quarks are treated entirely as massive quarks in 
the fixed flavour number scheme | 158 |. Thus, = 3 in the DGLAP splitting functions. 



and only gluon and singlet quark distributions are evolved. 

The structure functions F2 and F^^^^ are then given by the singlet quark distribution 
and a massive charm quark contribution due to boson-gluon fusion. Explicit formulae can, 
for example, be found in ||156|| . For the numerical studies, the values ALo,n/=3 = 144 MeV 



{as{Mz) = 0.118), rric = 1.5 GeV, = 4.5 GeV are used, and the massive charm quark 
contribution is evaluated for a renormalization and factorization scale fic = 2mc. 

The resulting structure functions can be compared with HERA data on the in- 
clusive structure function F2{x,Q'^) ||159|| and on the diffractive structure function 



F^^^\C,, PjQ"^) [^,0. These data sets from the HI and ZEUS experiments are used 
to determine the unknown parameters of the model. The following selection criteria are 
applied to the data: x < 0.01 and ^ < 0.01 are needed to justify the semiclassical 
description of the proton; with Qq being a fit parameter, a sufficiently large minimum 
Q'^ = 2 GeV^ is required to avoid that the data selection is infiuenced by the current value 
of Ql; finally > 4 GeV^ is required in the diffractive case to justify the leading-twist 
analysis. 

The optimum set of model parameters is determined from a minimization of the total 
(based on statistical errors only) of the selected data. As a result, 

Ql = 1.23 GeV^ , 



'0 

L = 8.16, 

n = (712 MeV)- 

a = (74.5 MeV) 



(7.6) 

2 



All parameters are given with a precision which allows reproduction of the plots, but 
which is inappropriate with respect to the crudeness of the model. The distributions 
obtained with these fitted parameters yield a good qualitative description of all data 
on inclusive and diffractive DIS at small x, as illustrated in Figs. p.l| , |7.1| and |7l2 . 



The starting scale Qq is in the region where one would expect the transition between 
perturbative and non-perturbative dynamics to take place; the two other dimensionful 
parameters QL^ and a are both of the order of typical hadronic scales. 

The approach fails to reproduce the data on ^2^*^^'' for low M^. This might indicate 
the importance of higher twist contributions in this region (see, e.g., ||3^). It is interesting 
to note that a breakdown of the leading twist description is also observed for inclusive 
structure functions ||160|| , where it occurs for similar invariant hadronic masses, namely 
W^<A GeV\ 

The perturbative evolution of inclusive and diffractive structure functions is driven 
by the gluon distribution, which is considerably larger than the singlet quark distribution 
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in both cases. With the parameters obtained above, it turns out that the inclusive gluon 
distribution is about twice as large as the singlet quark distribution. By contrast, the 
relative magnitude and the (3 dependence of the diffractive distributions are completely 
independent of the model parameters. Moreover, their absolute normalization is, up to 
the slowly varying factor l/as{Ql), closely tied to the normalization of the inclusive 
gluon distribution. 

In spite of the (1—/?)^ behaviour, gluons remain important even at large /?, simply due 
to the large total normalization of this distribution. As a result, the quark distribution 
does not change with increasing for (3 ~ 0.5 and is only slowly decreasing for larger 
values of /?. 

The dependence of the diffractive structure function on /5 and is illustrated in 



Fig. |7.3| , where the predictions are compared with experimental data [|73|,|7^ at fixed 
^ = 0.003 (HI) and C, = 0.0042 (ZEUS). The underlying diffractive parton distributions 
dg/d^ and d'E/dC, = 6dq/d^ at input scale Qq and their evolution are shown in Fig. |6.3| . 

It is an essential feature of the above semiclassical analysis that the rise of F2{x, Q^) 
and of i^2^^^^(C, /3, Q^) small x and ^ has the same, non-perturbative origin in the energy 
dependence of the average over soft field configurations in the proton. With increasing 
Q^, the X dependence is enhanced by perturbative evolution in the case of the inclusive 
structure function while, in the diffractive case, the ^ dependence remains unchanged. 



Importance of higher twist 



To illustrate the importance of higher twist contributions, consider the parametrization 
of Bartels et al. ||3^ (cf. Figs. [7.4| and [7.5|) . In this analysis, four different contributions 
to the diffractive structure function F.F^^^ are considered. These are and - the 

z qq qqg 

leading twist contributions of qq and qqg components of the transverse photon wave 
function - supplemented by AF^^ and AF^ - the higher twist contributions of the qq 
component of longitudinal and transverse photon wave function. 



As has already been explained in Sect. p.2| , there is no leading twist contribution to 
diffraction from the qq component of the longitudinal photon. Nevertheless, this compo- 
nent is important because of its dominance in the region (3^1. 



The approach of [0 is to write down generic expressions for the above four contribu- 
tions to the structure function. These expressions incorporate the qualitative knowledge 
about /5, ^ and dependence of each contribution. This knowledge is based on the two 
gluon exchange calculations of pl| , |101|| and is largely in agreement with what can be 
learned from the semiclassical treatment presented in this review. For example, the lead- 
ing twist contributions are expected to have a softer 1/^ behaviour than the higher twist 
contributions, which are expected to behave as {^giC])"^ (cf. Sect. ^.3|) . Furthermore, the 
longitudinal qq component is the only one known not to vanish at t = and /? — >■ 1. 

The relative weights of the four contributions, which are, in principle, determined by 



photon wave function and structure of the target, are allowed to vary in the fit of p4| 



The results of the two independent fits to ZEUS and HI data (cf. Figs. and [731) are 



97 



Mv= 5 GeV 



M,= 1 1 GeV 




100 



200 



100 



200 



100 



200 



Q = 14GeV' 
150 

100 















1 1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 1 





100 



200 



100 



200 



100 



200 



Q'= 27GeV' 
50 




100 



200 



100 



200 



100 



200 



Q'= 60GeV 

10 r- 




Figure |7.4 

result 
also [g 



Fit to ZEUS data for F2 

dashed line: FT^; dotted line: FT 



by Bartels et al. g. Solid line: total 
ggg, dashed-dotted line: AF^^ (data from ||161|] , see 



98 




99 



in good agreement with physical expectations. Note especially the dominance of 
at /9 — > 1 and of at /5 — > 0. Note also the steep ^ dependence of AF^, the need 



<?<?' 



of which is visible, in particular, in the top right-hand bins of Fig. |7.5| . Thus, the data 
clearly has room for the higher twist contributions introduced in . It would certainly 
be desirable to include such contributions in the leading twist analysis of with its 
consistent treatment of the evolution. An obvious problem to be solved before such 
a development can be realized is the ^ dependence of leading and higher twist terms. 



As an interesting feature of the data, the authors of |34| notice that the HI results 



can be described by two different sets of parameters: one in which is enhanced at 



large (3 (corresponding to Fig. |7.5|) and a second one which is more similar to the ZEUS 
fit. The solution with enhanced qqg contribution can be interpreted as the analogue, 
within the framework of this parametrization, of the singular gluon proposal of HI ||73|| . 



As can be seen from the QCD analysis in |]73[, the HI data exhibit some preference 
for a gluon distribution that is large for high such that for low values of most of the 
pomeron's momentum is carried by a single gluon. This singular gluon proposal, which 
has some similarity with the naive boson-gluon fusion model of is hard to justify 
theoretically. In particular, the validity of the partonic leading twist picture at /3 ~ 1 and 
not too large is questionable. Given the successful analysis of [^] and the alternative, 
less singular parametrizations in p^fTSf 



a singular diffractive gluon distribution is not 



unavoidable with the present data. 



Note that the QCD analysis of [|7^] has the important advantage of explicitly deal- 
ing with the Reggeon contribution, which, although suppressed in the high-energy limit, 
represents an important correction in the upper region of the ^ range. A Reggeon con- 
tribution is also included in the recent, more detailed analysis of HI data ||162|| in the 
framework of the parametrization of p3 . For a discussion of Reggeon exchange in diffrac- 
tive DIS the reader is referred, e.g., to [|163|| . 



A further interesting analysis, based on the model of the stochastic vacuum (cf. 
Sect. |67^ ) and including both leading and higher twist contributions, was reported by 
Ramirez . As an example, results for F^^^'' and for the transverse contribution F^^^^ 
at = 12 GeV^ and x = 0.0075 are shown in Fig. 7.6 



The presented calculation is limited to the qq fluctuation of the photon and does 
not include Altarelli-Parisi evolution. Therefore, the author limits the discussion of data 
to the high-/3 region, where gluon radiation effects are not yet dominant. The model as 
used in previous analyses tends to overshoot the data. Modifications proposed by Rueter 
(see [ ^ , [164] ]) subtract integration regions where one of the quarks is too slow for the 
eikonal approach to work, and where the transverse size of the pair is too small for the 
non-perturbative model to be applicable. The modified model gives a good description 
of the data (cf. Fig. fTBi ). 

It is particularly interesting to see the absolute magnitude of the longitudinal struc- 
ture function at large /3 as a prediction of an explicit calculation. Furthermore, the model 
shows a no n- vanishing contribution from the qq component at /5 — > 1. This is an improve- 
ment compared to the treatment of which is related to the better description of the 
dynamics of the outgoing proton at non-zero t in the calculation of . 



100 
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Figure |7.6| : Diffractive structure functions at = 12 GeV^ and x = 0.0075 in the model 
of the stochastic vacuum, calculated by Ramirez compared to HI data from 



^F^^^\^, f3,Q'^) and ^F^^"^' (3,Q^) in the modified model (long-dashed and solid lines 
respectively); ^F^^^\c,, /SjQ"^) in the model without modifications (short-dashed line); 
figure from 15^ . 
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Note that the above approach allows, in principle, for the inclusion of explicit qqg 
contributions of the photon wave function as well as for the all orders resummation of 
InQ^ terms, in a way similar to ||6^. It would be very interesting to see the results of 
such an extended analysis. 



The analysis of Bertini et al. I^S] is devoted to the influence of higher twist terms 
on the dependence of diffractive DIS. It is found that, similar to what is known 
about the longitudinal qq contribution, the higher twist terms in the transverse cross 
section are short distance dominated and are therefore perturbatively calculable. These 
contributions affect the dependence at /3 — >• 1 and should be subtracted from the 
cross section before a conventional DGLAP analysis of the large (3 region is performed. 

The data analyses of [|165| , p.66|| are based on the dipole picture of the BFKL pomeron 



discussed briefly in Sect. The diffractive structure function is described as a sum of 
two components: the elastic component, where the qq fluctuation of the virtual photon 
scatters off the proton ||140|| , and the inelastic component, where, based on the original 



qq fluctuation, a multi-gluon state develops, giving rise to high-mass diffraction [37 



While the first contribution is modelled by BFKL pomeron exchange in the t channel, 
the second one involves the triple-pomeron coupling and a cut gluon ladder representing 
the final state. 



In Ref . ||165|| , a combined analysis of both inclusive and diffractive structure functions 
was performed. The inclusive structure function was fitted on the basis of a colour dipole 
BFKL calculation, which, however, contains a number of free parameters and a smaller 
pomeron intercept than naive perturbation theory would suggest. Introducing further 
parameters for the normalization of the elastic and inelastic components discussed above, 
quahtative agreement with both the diffractive and inclusive DIS data was achieved. 
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A further phenomenological investigation of diffraction on the basis of colour dipole 
BFKL was performed in [p.66|| , where a good 7 parameter fit to the diffractive structure 
function was presented. It is certainly a challenging problem to gain a better understand- 
ing of the dynamics, thus reducing the number of required parameters. In particular, the 
applicability of perturbative methods needs further investigation. 

Very recently, an analysis of diffraction has appeared [|147|| that is based on the pre- 



vious inclusive DIS fit with a Glauber type model for the dipole cross section a{p) | |146 |. 



Motivated by the idea of saturation (compare the discussion in ||167|| ), an energy or x 
dependence of a{p) is introduced according to the formula 



a{x, p) = (To i 1 - exp , \ , (7.7) 




where the function Rq{x) vanishes as a: — > 0. The authors find a similar energy depen- 
dence of the diffractive and inclusive cross section. Two new parameters required for the 
analysis of diffraction are the diffractive slope and the fixed coupling constant a^. The 
very successful description of F^f is based on the qq and qqg fluctuations of the virtual 
photon. It includes higher twist terms but no Altarelli-Parisi evolution. 

In the different approaches discussed above, separate predictions for transverse and 
longitudinal photon cross sections can be derived. Corresponding measurements, which 
are, however, very difficult, would be immensely important in gaining a better under- 
standing of the underlying colour singlet exchange. An interesting possibility of obtaining 
the required polarization information is the measurement and analysis of the distribution 
in the azimuthal angle, which, at non-zero t, characterizes the relative position of the 



leptonic and hadronic plane of a diffractive event ||168|| 



7.2 Final states 

A further tool for studying the mechanism of diffraction is provided by the details of the 
hadronic final states. In this section, the focus is on events with a diffractive mass of the 
order of or larger. Diffractive production of single particles forms the subject of the 
next section. 

The first topic is the production of jets with high transverse momentum relative 
to the 7*p collision axis within the diffractively produced final state. Such diffractive 
dijet events were reported soon after the beginning of the investigation of diffraction 
at HERA [169[ |. Given the kinematic constraints and the available number of events, a 
direct analysis of jet cross sections in diffractive DIS proved difficult. Therefore, event 
shape observables such as thrust were first considered by both collaborations [|17CI|] . One 
of the main results was the alignment of the diffractive state along the 7*p axis and, 
furthermore, the presence of significant transverse momentum components. If this is to 
be interpreted in terms of diffractive parton distributions, a large gluonic component 
seems to be required. 

A direct analysis of diffractive dijet events has recently been reported by the HI 
collaboration ||171|| . In the following, this analysis is discussed in more detail since, for 
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Figure [7.7| : Differential cross section in transverse jet momentum for events where the 
diffractively produced hadronic state contains two jets. Photon virtualities vary in the 
range 7.5 < < 80 GeV^. The data are compared to predictions of RAPGAP Monte 
Carlo models with leading order pomeron and reggeon parton distributions: quarks only 
(fit 1), 'flat' gluon dominated (fit 2) and 'peaked' gluon dominated (fit 3). The figure is 
from ||171 | 



where further details can be found. 



such measurements, a qualitative comparison with different theoretical models is simpler 
and less affected by hadronization effects. As an illustration, one of the distributions 
pubhshed in [ p.71| | is reproduced in Fig. fLl\ . 

As explained in Sect. ^.3| , jets in the diffractive final state can be the result of either 
boson-gluon fusion, with the gluon coming from the diffractive gluon distribution of the 
proton, or of so-called exclusive dijet production, where the jets are associated with the 
pure qq component of the photon wave function. The former case, viewed in the target 
rest frame, is associated with the qqg component of the photon wave function, where the 
gluon has small p± and the colour singlet exchange in the t channel is soft. 

In Fig. \r7[ the data points are compared to predictions based on the boson-gluon 
fusion scenario. They are well described if the diffractive parton distributions are domi- 
nated by the gluon, a situation reahzed, e.g., in the semiclassical calculation based on the 
colour field model of Sect. |6]^. It would certainly be interesting to substitute the diffrac- 
tive parton distributions underlying the Monte Carlo results of Fig. |7.7| (see, e.g., [|I72 



isee, e.g , 

for a recent review of Monte Carlo generators used for HERA diffraction) with the model 
distributions of Chapter ||. However, given the precision of the data, which is at present 
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consistent with the two very different gluon distributions of fit 2 and fit 3, it appears 
unhkely that any gluon dominated model can be ruled out. 

The most important qualitative conclusion to be drawn is that the data favours 



boson-gluon fusion rather than exclusive dijet production. As discussed in Sect, p^ , the 
exclusive process has a much steeper p± distribution, which would be reflected in the 
quality of the fit in Fig. [7.7| . This conclusion is made even clearer by the distribution 
of the variable z^^^, which characterizes the fraction of the proton's momentum loss 
absorbed by the jets. In exclusive jet production, z-^p^ = 1, which, even taking into 
account the smearing caused by hadronization, is in sharp contrast to the data. Thus, 
the observations suggest that only a small fraction of the observed dijet events are related 
to exclusive jet production with hard colour singlet exchange. It would be important to 
try to quantify this statement. 

A further interesting aspect of the diffractive final state is the presence of open 
charm. So far, no published results by either the HI or ZEUS collaboration are available 
on this subject. However, preliminary data have recently been presented by both col- 
laborations ||173| , |174j] . Cross section measurements for the production of D*^ mesons in 
diffractive events are reported. It is found in ||174|| that both and dependences of 
diffractively produced open charm as well as the fraction of charmed events are consistent 
with the observations made in inclusive DIS. This is in agreement with the semiclassi- 



cal scenario of or, equivalently, with a boson-gluon fusion scenario based on gluon 



dominated diffractive parton distributions. 

It should finally be mentioned that a number of more inclusive diffractive final state 
analyses focussing, in particular, on charged particle distributions ||175||, transverse en- 



ergy flow ||176|| and the multiplicity structure ||177|| were reported by both collaborations. 



Although a detailed discussion of the obtained results can not be given here, it is im- 
portant to note that all observations seem to be consistent with the physical picture 
of a partonic pomeron. Clearly, this is also in accord with the more general concept 
of diffractive parton distributions. As expected, particle distributions and energy flow 
support the by now well-established picture of colour singlet exchange with the target. 
In agreement with the jet oriented analyses, a significant gluonic component is required 
by the data. It was emphasized by HI (see ||175| , |177| ) that a relatively hard diffractive 



gluon distribution is favoured. 

As already mentioned in the Introduction, an approach describing diffractive DIS 
by the assumption of soft colour interactions in the final state was proposed by Edin, 
Ingelman and Rathsman in EO]. The idea of soft colour exchange is similar to the naive 



boson-gluon fusion model of |39|, but the Monte-Carlo based implementation is quite 
different. After the hard scattering, the physical state is given by a number of quarks 
and gluons - those created in the scattering process and those from the proton remnant. 
It is assumed that, at this point, colour exchange can take place between each pair 
of colour charges, the probability of which is described by a certain phenomenological 
parameter R. This changes the colour topology and leads, in certain cases, to colour 
singlet subsystems, giving rise to rapidity gaps in the final state. The number of gap 
events initially increases with R, but saturates or even decreases at larger values of this 
parameter since further colour exchanges can destroy neutral clusters previously created. 
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A more detailed description of the soft colour interaction model can be found in ||178 



where the consistency of the approach with the large transverse energy flow in the proton 
hemisphere, a feature of HERA data that appears to be 'orthogonal' to the rapidity gaps, 
is demonstrated. In a more recent development ||179|| , the colour exchange is described on 
the level of colour strings rather than on the partonic level discussed above. Furthermore, 
e"'"e~ data is included in the analysis. In conclusion it can be said that an impressive 
agreement with many features of hadronic final states in small-x DIS has been achieved 
within the framework of soft colour interactions. It is an interesting question in how far 
the perturbatively known suppression of interactions with small colour-singlet objects, 
i.e., colour transparency, is consistent with this approach. 



7.3 Meson production 

Exclusive meson production in DIS is, both theoretically and experimentally, a very 
interesting and active area that certainly deserves more space than is devoted to it 
in the present review. The main theme of the previous chapters was the semiclassical 
approach, the applicability of which to meson production processes is limited by the 
lack of a genuine prediction for the energy dependence of the amplitude. Although, as 
discussed in Sect. ^A\ , Ryskin's double-leading-log result is reproduced, two essential 
problems in going beyond it, the dynamics of the energy growth and details of the meson 
wave function, have not been addressed. Nevertheless, this review would be incomplete 
without at least a brief discussion of the most important new experimental results and 
their implications for theory. 

Recent analyses of exclusive vector meson electroproduction have been published by 
ZEUS for p and J/ ip mesons ||180|] and by HI for p mesons ||181|| . New preliminary results 



on exclusive J/tjj production were also reported by HI [|182|| , with previous measurements 



published in ||183||. Agreement is observed between the most recent results of the two 



experiments [ |181 |. 



For illustration, the ZEUS results for a particularly interesting quantity, the forward 
longitudinal p meson cross section, are shown in Fig. [778| . As discussed in Sects. and 



572| , this quantity can, at least in principle, be calculated in perturbation theory. If the 
asymptotic form of the meson wave function (j)p[z) ^ z{l — z) is used, the amplitude of 
Eq. (|5.13|) gives rise to the cross section formula [p!6t|107|| 



da 



t=o acmQ^ 



The main qualitative predictions of this formula are in good agreement with experi- 
mental results. Taking into account the anomalous dimension of the gluon distribution, 
the observed dependence of roughly is consistent with Eq. ( [7.8|) . The data indi- 
cate that the energy dependence becomes more pronounced with growing Q^, as expected 
from the square of the gluon distribution. However, taking the above formula at face value 
and using, e.g., the MRS(A') gluon distribution [ |184|| , the absolute normalization of the 
cross section comes out too high. 
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Figure [7.8| : Forward longitudinal p electroproduction cross section da/dt\t=Q as mea- 
sured by the ZEUS collaboration ||180|| . Shaded areas indicate normalization uncertain- 
ties in addition to the error bars shown. The curves are based on calculations by Martin, 
Ryskin and Teubner | |126|| (solid line), by Frankfurt, Koepf and Strikman ||116|| (dashed- 
dotted line) and on colour dipole model calculations of Nemchik et al. |[L17|| (dashed line). 
The figure is from 



180 , where further details can be found. 



Several improvements on the result of [|1^] were discussed in ||116|| . In particular, it 
was found that the Fermi motion of the quarks in the vector meson leads to a suppression 
factor in the cross section. Furthermore, in comparing the recent ZEUS data with the 
perturbative calculation, the gluon distribution is evaluated at a scale Ql^ < Q^. The 
new scale is determined from the average transverse size of the qq pair, which turns out 
to be larger in the p production process than, say, in longitudinal electroproduction. 
Such a rescaling is certainly legitimate in view of the leading logarithmic nature of the 
approach. It leads to a further reduction of the cross section. While both suppression 
effects together produce results that are too small , either one of them gives a reasonable 
fit (see Fig. and ||180|| ). Note also the claim ||185|| that, at least in the J/ip case, Fermi 
motion effects might be considerably smaller than expected. 
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The approach of ||126|| , briefly mentioned in Sect. also gives a reasonable fit to the 
p production data. Both the results of FKS ||116|| and MRT ||126|| shown in Fig. |7.^ are 
based on the ZEUS 94 NLO gluon distribution, which, in the relevant parameter region, 
is not too different from, e.g., MRS(A') ||184|| . Note that all calculations which derive 
their total normalization from the identification of the two gluon exchange amplitude 
with the gluon distribution are affected by the large experimental uncertainties of the 
latter. 



The results of ||117|| , also displayed in Fig. [778| , are based on calculations that em- 



phasize the non-perturbative aspects of the process. The interaction of the qq pair with 
the target is described in terms of the colour dipole cross section a^p,^), where u is 
the photon energy. The amplitude for the production of a qq pair of transverse size p is 
convoluted with a meson wave function with non-trivial p dependence. In the relevant re- 
gion, where is not yet very large, these non-perturbative effects and the contributions 
from a(p, u) at relatively large p are important. It is proposed that the energy and 
dependence of the measured cross section will allow for the 'scanning' of the dipole cross 
section cr(p, u) as a function of p (cf. | ]186|| ) and provide evidence for BFKL dynamics. 

The violation of s channel helicity conservation is found to be small. However, evi- 
dence for a helicity flip amplitude on the level of 8 ± 3% of the non-flip amplitude was 
reported in ||181|| . A particularly interesting feature of the data is the relatively small 
ratio R = ai/cTT of longitudinal and transverse cross sections, which reaches the value 
20GeV2 |8T1. This is in sharp contrast to a naive extrapolation of 
leading one to expect R ~ Q"^ /m^. Even though 



i? ~ 3 at 

the perturbative result of Eq. 
the transverse cross section is not perturbatively calculable in the case of light vector 
mesons, one would still expect a linear growth of R with Q^, which is not favoured by 
the data 
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In the stochastic vacuum approach of |]56|-[58[], the supposition that, at presently 
accessible values of Q^, the perturbative regime is not yet reached is taken even more 
seriously. The analysis is based on a fundamentally nonperturbative model of the in- 
teraction with the target hadron (cf. Sect. |0| ) and is thus ideally suited for studying 
the transition region to photoproduction, as well as the transverse cross section, where 
QCD factorization theorems fail. Clearly, the approach depends on the model for the 
p meson wave function, but certain observables, such as the ratio of longitudinal and 
transverse cross sections or elastic slopes, represent relatively robust predictions. In par- 
ticular, the ratio of longitudinal and transverse p meson production cross sections is well 
described 
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In its original form, the model of the stochastic vacuum predicts constant high-energy 
cross sections. This limits the applicability of the model in the HERA regime, where the 
energy growth of meson electroproduction cross sections can not be neglected. Recently, 
the model has been extended by introducing a phenomenological energy dependence 
based on two different pomerons coupled to small and large dipoles, very much in the 



spirit of | [187|| . Furthermore, for small dipole configurations, a perturbative two gluon 
exchange contribution was added. This allowed for a good description of both photo- 
and electroproduction of vector mesons and of F2{x, Q"^) in the HERA regime. 
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The presently available data can also be described in the generalized vector domi- 
nance approach of ||188|| . In contrast to the calculations discussed so far, this approach 
predicts both longitudinal and transverse vector meson production cross sections to fall 
like ~ at asymptotically large virtualities. 



Note also the recent perturbative model calculation of ||118|| , which finds an asymp- 
totic behaviour of R = cxL/arp that is qualitatively different from the expected linear 

growth. The calculations of | ]118| , |126| , |188| reproduce the experimentally observed flat 
behaviour of R at high Q^. 

The most striking feature of the measured J/ip production cross sections is the strong 
energy growth ~ {W"^)^'^. Its approximate independence up to ~ 13 GeV^ suggests 
that, in the present data, the charm quark mass rather than the photon virtuality play 
the dominant role in making the colour singlet exchange hard. The flavour-symmetric 
ratio a{J /iIj) / a{p) = 8/9 is not yet reached at Q'^ = 13GeV^. 

Fitting the t dependence of p and J/ip electroproduction cross sections with the 
function e''*, values of 6 ~ 8 GeV"^ and 6 ~ 5 GeV~^ respectively were obtained 
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The exclusive electroproduction of excitations of both the p and the J/ ip meson have 
also received much theoretical interest. However, experimental data on this subject is 
only beginning to emerge [|189| . 
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8 Conclusions 



In the above review, different methods employed in the treatment of diffractive electro- 
production processes were put into context. 

The semiclassical approach, which was chosen as the starting point of the discussion, 
is well suited to develop a basic understanding of the underlying physical effects. In this 
approach, non-diffractive and diffractive DIS are treated along parallel lines. Diffraction 
occurs if, after scattering off the target, the partonic fluctuation of the virtual photon 
emerges in a colour singlet state. Altarelli-Parisi evolution of both the diffractive and 
inclusive structure function is related to the presence of higher Fock states of the photon. 
These fluctuations have a small transverse size and do not affect the soft colour exchange 
with the target. The application of the semiclassical approach to experimental results 
is particularly simple if the approach is used to derive both diffractive and inclusive 
parton distributions at some small input scale. In this case, the analysis of all higher-Q^ 
data proceeds with standard perturbative methods. Different models for the underlying 
colour fields can be compared to diffractive and inclusive structure function data in a 
very direct way. 

If the scattering amplitude of the partonic fluctuation of the photon and the target 
grows with energy, then, in the semiclassical framework, this growth enters both the 
diffractive and inclusive input distributions in precisely the same way. Such an energy 
growth is expected to be generated by the process of averaging over all relevant colour 
field configurations. However, no explicit non-perturbative calculation exists so that, 
at present, the energy dependence remains one of the most challenging aspects of the 
method. 

The idea of the pomeron structure function combines methods of Regge theory with 
the partonic description of hard scattering processes in QCD. In the treatment of diffrac- 
tive electroproduction, both the semiclassical approach and the pomeron structure func- 
tion are equivalent as far as the hard part of the process, i.e., the evolution, is 
concerned. In this respect, they represent different realizations of the more general and 
less predictive concept of diffractive parton distributions. However, as far as the colour 
singlet exchange mechanism is concerned, the pomeron structure function idea is very 
different in that it assumes the existence of a 'pre-formed' pomeron state, off which the 
virtual photon scatters. By contrast, in the semiclassical approach the diffractive charac- 
ter of an event is only determined after the photon fluctuation has passed the target. No 
pre-formed colour singlet object is required. Prom this point of view, the semiclassical 
description of diffraction shows similarity to the idea of probabilistic soft colour exchange 
added on top of standard Monte Carlo models. However, in contrast to this method, the 
semiclassical calculation keeps track of the relevant transverse distances in the parton 
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cascade, ascribing soft colour exchange only to large size configurations. 

In an apparently quite different approach, the t channel colour singlet exchange, 
characteristic of diffractive processes, is realized by the exchange of two t channel gluons. 
This is justified for certain specific final states, like longitudinally polarized vector mesons 
or two-jet systems, where the colour singlet exchange is governed by a hard scale. Prom 
the semiclassical perspective, this effect manifests itself in the small transverse size of 
the photon fluctuation at the moment of passing the external colour field. Thus, only 
the small distance structure of the field is tested, which corresponds to the two gluon 
exchange amplitude discussed above. 

The fundamental advantage of two gluon exchange calculations is the understanding 
of the energy dependence which, one may hope, can be achieved in this framework by 
the resummation of large logarithms. However, leading order resummation predicts a 
far too steep energy growth, and the recently obtained next-to-leading correction is very 
large, thus calling the whole method into question. It is also possible that the mechanism 
responsible for the energy dependence in soft processes is fundamentally different from 
the above perturbative ideas. 

Nevertheless, leading logarithmic two gluon calculations were successful in relating 
certain diffractive cross sections to the conventional gluon distribution, the energy growth 
of which is measured. Going beyond leading logarithmic accuracy, the non-forward gluon 
distribution was introduced. This is a promising tool for the further study of diffractive 
processes with hard colour singlet exchange. 

On the experimental side, much attention has been devoted to the diffractive struc- 
ture function, which represents one of the most interesting new observables in small-x 
DIS. Leading twist analyses, based on the concept of diffractive parton distributions, are 
successful in describing the bulk of the data. It has been demonstrated that a simple large 
hadron model for the target colour fields allows for the derivation of a consistent set of 
diffractive and inclusive parton distributions in the semiclassical framework. Other colour 
field models, such as a perturbative dipole field and the model of the stochastic vacuum, 
promise a successful data analysis as well. It can be hoped that future work with the 
data will expose the qualitative differences among the available models, thus advancing 
our understanding of the proton bound state and non-perturbative QCD dynamics. 

The available data shows the importance of higher twist contributions in the small- 
mass region of diffractive structure functions. A consistent theoretical description of 

both leading and higher twist effects has yet to be developed. Furthermore, the energy 
dependence of the diffractive structure function and, even more importantly, its relation 
to the energy dependence of the inclusive structure function are interesting unsolved 
problems. While the naive summation of logarithms appears to be disfavoured, no new 
standard framework has emerged. The proposal of an identical, non-perturbative energy 
dependence of both diffractive and inclusive cross sections at some small virtuality Qq 
is phenomenologically successful. This energy growth is reflected in the ^ dependence of 
the diffractive structure function, which is not altered by the perturbative evolution. 

The detailed investigation of diffractive final states provides a further broad field 
where different approaches can be tested. Here, considerable improvements on present 
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data are expected in the near future. So far, the analyses performed emphasize the 
dominance of the soft colour exchange mechanism and the importance of the diffractive 
gluon distribution. Diffractive meson production is, on a qualitative level, well described 
by perturbative QCD wherever the dominance of the hard scale is established. However, 
large uncertainties remain, and non-perturbative contributions are important even at 
relatively high values of Q^. In particular, the ratio of longitudinal and transverse p 
meson production cross sections is far smaller than suggested by simple perturbative 
estimates. A better understanding of the dynamics of meson formation appears to be 
required. 

In summary, it is certainly fair to say that diffractive electroproduction at small x 
proved a very rich and interesting field. Its investigation over the recent years has lead, 
in many different ways, to the improvement of our understanding of QCD dynamics on 
the interface of perturbative and non-perturbative physics. However, many important 
problems, such as the energy dependence and the systematic treatment of higher twist 
contributions, remain at present unsolved. 
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Appendix 

A Derivation of Eikonal Formulae 



A particularly simple way to derive the eikonal formulae of Sect. |3.1| is based on the direct 
summation of all Feynman diagrams of the type shown in Fig. 3^ in the high-energy 
limit. In the case of scalar quarks, the one gluon exchange contribution to the S'-matrix 
element reads 

Siip',p) = -19{P,+V') I d'xA>^ix)e'^^^P'-P^ . (A.l) 



Working in a covariant gauge and in the target rest frame, all components of A are 
expected to be of the same order of magnitude. If, in the high-energy limit, p+ and p'_^_ 
are the large momentum components, the approximation (p^ +p'^)A^ ~ (p+ + p'_^_)A-/2 
can be made. Assuming that the x dependence of A is soft, one can write 

rfx_e^^-(PV-P+)/2^_(a;+,x_,xx) ^ An 6{p'^ - p+) A4x+,x±) , (A.2) 

where the x_ dependence of A has been suppressed on the r.h. side since it is irrelevant 
for the process. More precisely, one has to think of the incoming and outgoing particles 
as wave packets, localized, say, at x_ ~ 0, which is the x__ position where A has to be 
evaluated. Equation ( |A.1|) takes the form 

S^{p',p) = 2ir6{p'^-p+) 2p+J dVe-^^^(^^-^^) (-| / rfa:+A_(a;+, x^)) , (A.3) 

which, within the approximation k^6{k^) ~ A;o'^(^o); is precisely the first term of the 
expansion of Eq. ( p.3| ) in powers of A. 

The n gluon exchange contribution to the S'-matrix element can be written as 

^ nl-f/<^v-./i|^(p?+pr').-(.«)ViipTir 



1=2 



^ / (rf^ +pf ) A_(x«)e-'^^(^'^'-^'°M , (A.4) 



I 2 

where x^^^ is the space-time variable at vertex i, the momentum of the quark line between 
vertex i and vertex i+1 is denoted by p^^\ and the initial and final momenta are p = p^^^ 
and p' = p^^\ With Eq. ( |A.2[) , all p+ integrations in Eq. ( A.4 ) become trivial, and all 
momentum plus components become identical with p^^ . Integrations over the momentum 
minus components are performed using the identity 

g-ip_{j/+-2+)/2 2'jii 

dp- = -Oiy+ - z+) — . (A.5) 

P+P- - PI + «e P+ 
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The resulting step functions translate into path ordering of the matrix valued fields A 
along the plus direction, 

J dx^l\.. dxf6{xf - x^l^)... e{xf - ■■} = ^J dxl\.. dxf P{- ■ •} . (A.6) 

Now that all dependence on the intermediate transverse momenta pj^ ...p]^ origi- 
nating in the propagators has disappeared, the corresponding integrations are straight- 
forward. They produce (5-functions in transverse co-ordinate space, 

£p^e'P^iv^-'^) = (2nf 5\y^ - za_) . (A.7) 



Making use of Eqs. (|A.5| ) - (|A.7|) , the expression for the nth order contribution to 



the ^'-matrix element finally takes the form 

= 27r5(p'+-p+)2p+|| rfVe-'"^(^^-P^)}^p(-| / dx+A.{x+,xS) • 

(A.8) 

Since this is precisely the nth term of the expansion of Eq. (|3.3|) in powers of A, the 
derivation of the eikonal formula in the scalar case is now complete. 

To obtain the eikonal formula for spinor quarks, Eq. ( |3.5| ), write down the analogue of 
Eq. (|A.4|) using the appropriate expressions for quark propagator and quark-quark-gluon 
vertex. Applying the identity of Eq. ( |3.11| ) to each of the quark propagators, neglecting 
the 7+ term, which is suppressed in the high-energy limit, and making use of the relation 

Us'ik'){-ig4)us{k) ^ -igk+A^ 6ss' , (A. 9) 



valid if fc+ ^ k'_^ are the big components, the exact structure of Eq. ( |A.4| ) is recovered 



The only difference is an additional overall factor Sss', corresponding to the conservation 
of the quark helicity. The further calculation, leading to the analogue of Eq. ( |A.8| ), is 
unchanged, and the eikonal formula follows. 

The further generalization to the gluon case, Eq. ( |3.(j| ), proceeds along the same lines. 
Begin by writing down the analogue of Eq. ( |A.4| ), i.e., the S*- matrix element for Fig. |3]2 



with a gluon line instead of the quark line. The expression contains conventional three- 
gluon vertices and gluon propagators in Feynman gauge. Contributions with four-gluon 
vertices are suppressed in the high-energy limit. The colour structure is best treated by 
introducing matrices 

= A''{T;\f' = -iAT^^ , (A.IO) 
where f°''"^ are the usual structure constants appearing in the three-gluon vertex. Now, 



the product of matrices A in Eq. ( |A.4| ) is, in the gluonic case, simply replaced by an 
identical product of adjoint representation matrices A^. 

Next, the g'^'^ tensor of each gluon propagator —ig^'^/k'^ is decomposed according to 
the identity 

a'^^-(Ye^e^ +^ + ^- ^kA (A 11) 
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where m is a light-like vector with a non-zero minus component, m = (0, 2, 0^), and the 
polarization vectors are defined by ek = em = and = —1. An explicit choice is given 
by 

e« = (0,^^^^,e(,,), (A.12) 

where the transverse basis e(i)^ = (1, 0) and e(2)± = (0, 1) has been used. If, in the 
high-energy limit, g^'^ appears between two three-gluon vertices, the last three terms of 
Eq. ( [A.ll[ ) can be neglected. Note that, for the second and third term, this is a non- 



trivial statement since the vector k in the numerator could, in principle, compensate for 
the suppression by the fc+ in the denominator. However, this is prevented by the gauge 
invariance of the three gluon vertex. 

Thus, the analogue of Eq. ( |A.4|) is written down with three-gluon vertices and prop- 



agators proportional to 

a'-'^Ee^fl.). (A.13) 

i=l 

Now, each of the three-gluon vertices Vfj_ucr{—k', k, k'—k) (where all momenta are incoming 
and colour indices are suppressed) appears between two transverse polarization vectors 
and simplifies according to 

e^^^,^{k')V^,^{-k',k,k'-k)A''^e'(,){k) ~ -igk+AA^S,i, , (A.14) 

where is the external field. This is similar to what was found in the quark case 
in Eq. ( |A.9| ), with the difference that helicity conservation is replaced by polarization 



conservation. Therefore, as before, the structure of Eq. ( |A.4]) is recovered, but with the 



external field in the adjoint representation and with an additional polarization conserving 
(5-function. The further calculation, leading to the analogue of Eq. ([A.8|) , is unchanged, 
and the eikonal formula in the gluonic case follows. 
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B Spinor Matrix Elements 



In this appendix, the spinor matrix elements of the type Us'{k')^{q)us{k), required, e.g., 
for the calculation the transition from virtual photon to qq pair, are listed. 

Using light-cone components for vectors, a = {a^,a-,a±), and orienting the pho- 
ton momentum along the positive 2;-axis, q = {q^,—Q^/q^,0±), the longitudinal and 
transverse polarization vectors can be defined as 

eL = {q+/Q,Q/q+,0±) , e±i = (0, 0, e^(±)) , (B.l) 
where e±{±) = {l,±i)/V2. 



The Dirac representation of 7 matrices and the conventions of ||7g] for Dirac spinors 
are used. Introducing the matrix 



1 
-1 



(B.2) 



and the two-component spinors 



the two independent Dirac spinors, with s = ±1/2, for negative and positive frequency 
solutions can be explicitly written as 

( Xs \ / ak \ 

Us{k) = ^ko + m , Vs{k) = -^Jko + m\ t^o + m ^ (b.4) 

\ ko + m^' I \ exs I 

where a is the vector formed by the three Pauli matrices. 

The result for longitudinal photon polarization can be obtained from the relation 



Us'{k')'yoUs{k) = 6s's'2\/k)K, (B.5) 

which holds at leading order in the high-energy expansion. From this, the matrix ele- 
ment Us'ik')^ L{Q)us{k) is obtained using the gauge invariance of the gg-photon vertex, 
Usi{k')^Us{k) = 0. Since the quark mass and transverse momentum do not enter the lead- 
ing order relation, Eq. (|B.5|) , they may be neglected so that the relation Vs{k) = —u^sik) 
holds. Thus, the complete result takes the form 

Usi{k')tLUs{k) = v^,{k')eLVs{k) =^^s'{k')eLUs{k) =^_s'{k')eLUs{k) = 5^/^2(5 y'a(l-a) , 
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(B.6) 

where a = ko/qo and = ei^q). 

Neglecting terms suppressed in the high-energy hmit, the relevant matrix elements 
for transverse photon polarization read 



Us'{k')(f+i,i{q)u,{k) 



kok'o 



Vs'{k')(^ +1 1 _i{q)vs{k) 



kok'o 



Us'{k')(f+i/-i{q)vs{k) 



+ 1 



Vs'{k')(f+i/-i{q)us{k) 



+ 1 



kok'Q 



' k±+k'(^ 1 k'^_ko 


m{k'o-ko) / 


y / m{ko - k'o) 
' k'j^_^ko 1 k^_k'^ 


^±+^0 / k±-kQ 
/ m{k',-ko) 


y m{ko -k'^) 1 
^^{k^ + k'^) / 


k±+kQ / k'j__kQ 
kj_+kQ / k'^_ko 


y k'^^ko 1 k_L_k'Q 
/ / -m{ko + k'Q) 


/ m{k'Q + ko) 
^j_+^o / k±—kQ 


V k±+k'Q 1 k'j^_ko 


m{k'f^ + ko) 1 



(B.7) 



(B.8) 



(B.9) 



.(B.IO) 



Here the notation k^± = /ci ± ik2 is used, and the four entries (11), (12), (21) and (22) 
of the two- by-two matrices on the r.h. sides of the above equations correspond to the 
combinations (s's) = {^\ -|- |) , (+| — |) , {—\ + |) and (— | — |). For each of these 
four entries, the expression before and after the oblique stroke ' / ' corresponds to positive 
and negative photon polarization respectively. For example. 



n+i(fc')/+i(g)M_i(A;) = -^^^ "^(^0 - ^o) 



(B.ll) 



The calculation of the transition from virtual photon to gg-gluon configuration 
requires, in addition to the gg-photon vertex, the knowledge of the gg-gluon vertex 
(see Sect. p.3|) . The corresponding spinor matrix elements are easily obtained from 
Eqs. ( |B.7| ) — ([B.lOp if the configuration is rotated in such a way that the gluon mo- 
mentum is parallel to the 2;-axis. In the high-energy limit, such a rotation corresponds 
to the substitutions 



/c_L ^ fc± - g_L(fc+/g+) , fc^ ^ fc^ - g±(fc+/g+) 
where q is now interpreted as the gluon momentum. 



(B.12) 



Note that similar matrix elements are commonly used in light-cone perturbation 
theory (see, e.g.. Tables E and 11 of | |190| ]). 



116 



C Derivation of Diffractive Quark 
and Gluon Distribution 



The explicit formulae for diffractive quark and gluon distributions, Eqs. ( |4.53| ) and ( [4.54| ), 
can be derived by appropriately adapting the scalar calculation of Sect. to the case 
of spinor or vector particles. The most economic procedure is to first identify the piece of 
the old squared amplitude |Tp that depends explicitly on the spin of the soft parton. This 
piece, which is essentially just the squared scattering amplitude of the soft parton and 
the external field, is symbolically separated in Fig. |C.1| . The two independent integration 
variables for the intermediate momentum of the soft parton are denoted by k and k in 
the amplitude T and its complex conjugate T* respectively. Note also that k+ = k+ and 
k- = fc- at leading order. 







— - 




k\\\\\\\\\l 




k 






k 











B 



Figure |C.1| : Symbolic representation of the square of the amplitude for a hard diffrac- 
tive process. The box separates the contributions associated with the soft parton and 
responsible for the differences between diffractive distributions for scalars, spinors and 
vector particles. 



It is straightforward to write down the factor that corresponds to the box in 
Fig. |C.1| for the scalar case. Since the off-shell denominators fc^ and k"^ as well as the 
two eikonal factors and energy 5-functions are present for all spins of the soft parton, 
they are not included into the definition of B. All that remains are the explicit factors 
2fco — fc+ from the effective vertex, Eq. ( [4.35| ). Therefore, the result for the scalar case 
reads simply 



B, = ki 



(C.l) 



The next step is the calculation of the corresponding expressions Bq and Bg, given 
by the box in Fig. in the case where the produced soft parton is a quark or a gluon. 
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Introducing factors Bq/Bg and Bg/Bg into Eq. (|4.52|) will give the required diffractive 
parton distributions (df^/d^) and {df^ /d^). 



Diffractive quark distribution 



Observe first, that the analogue of Eq. ( |4.35|) for spinors is simply 



(C.2) 



The Dirac structure of Vg follows from the fact that, in the high-energy limit, only the 
light-cone component y4_ of the gluon field contributes. The normalization is consistent 
with Eqs. and (|l3]). 

Consider now the Dirac propagator with momentum k. The hard part Th requires 
the interpretation of the quark line as an incoming parton, which collides head on with 
the photon. Therefore, it is convenient to switch to the Breit frame and to define a 
corresponding on-shell momentum /, given by /_ = — l± = —k± and /+ = The 
propagator can now be written as 



1 _ T,sUs{l)Us{l) _ J_ 



(C.3) 



This is technically similar to Sect. p.2| , where an on-shell momentum was defined by 
adjusting the minus component of k. In the present treatment, however, the plus compo- 
nent of / = —k is adjusted, so that a partonic interpretation in the Breit frame becomes 
possible. It has now to be shown that the second term on the r.h. side of Eq. ( |C.3|) can 
be neglected, since it is suppressed in the Breit frame by the hard momentum /_ = yP-. 
This can be intuitively understood by observing that this term represents a correction 
for the small off-shellness of the quark, which is neither important for Th, nor for the 
soft high-energy scattering off the external field. 



To see this more explicitly, write the relevant part of the amplitude in the form 



u{k h+jTn = -u{k )7+ + — 



H 



-(C1 + C2) 



(C.4) 



It will now be shown that the 7_ term C2 is suppressed with respect to Ci. 
The spinor matrix element in Ci can be estimated using the relation 

4(pYp±) 



J2usip'h+Usip) 



PLP- 



(C.5) 



valid in the limit where the minus components of p and p' become large (cf. Table II 
of [119011 ). In the soft region, k"^ k\ k'1 A^, and — A;_ k'_ Q are the large 
components in the Breit frame. Thus, [u{k')'-)j^u{l)\ ~ A^/Q, and the first term on the 
r.h. side of Eq. (|C.4|) is estimated to be 



1 r 

Q 



u{1)Th 



(C.6) 
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where factors 0{1) have been suppressed. 

Introducing a vector a = (1, 0, 0_l), so that 7-/2 = ^, the 7_ term can be written as 



C2 = Y1 [u{k')-f+u,{a) — Us{a)TH . (C.7) 

Since [u{k')'y^u{a)] vanishes for k'^ = k'l/k'_ 0, in the Breit frame, where k'_ is large, 
the estimate [u{k')'y^u{a)] ~ \k'^\/ Jk'_ ~ A/ ^/Q can be made. With /_ ~ Q, one obtains 



C, A [uia)TH] A 

c.^VQ m)TH] ^ g ' ^ ■ ^ 

where it has been assumed that no specific cancellation makes u{T)Th small, i.e., 
[u{l)TH\/[u{a)TH\ ~ \/Q- Equation (|C.8|) establishes the required suppression of the 
7_ term in Eq. ( |(J.3| ). 

To proceed with the evaluation of Bs note that, when the soft part is separated in 
Fig. |(J.1| , the spinor Us{l) — V-s{l) from Eq. ( |C3| ) has to be considered a part of the hard 
amplitude Th- Therefore the analogue of Eq. reads 

Bs = Yl uAk') ^ us{l) Us{!) ^ uAk') , (C.9) 

s,s' 

where / is defined analogously to /, but using the momentum k instead of k. The spin 
summation decouples from the hard part if the measurement is sufficiently inclusive. 

The above expression can be evaluated further to give 

= lk^Y.^Mh+Usil) = (C.IO) 



where the last equality again uses Eq. ( p.5|) . Simple kinematics leads to the result 



Comparing this with Eq. (|C.1D , the diffractive quark distribution, Eq. ( [4.53|) , is straight- 
forwardly obtained from the scalar case, Eq. ( |4.52| ). 

Note that the virtual fermion line corresponds to a right-moving quark with momen- 
tum k in the proton rest frame and to a left-moving antiquark with momentum / in the 
Breit frame. Therefore, the above result has, in fact, to be interpreted as a diffractive 
antiquark distribution. The diffractive quark distribution is identical. 



Diffractive gluon distribution 

To obtain the diffractive gluon distribution, the procedure of the last section has to be 
repeated for the case of an outgoing soft gluon with momentum k' in Fig. |C.1|. Calculating 
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the contribution separated by the box will give the required quantity Bg, in analogy to 
Eqs. dai]) and 



It will prove convenient to introduce two light-like vectors m and n, such that the 
only non-zero component of m is m_ = 2 in the proton rest frame, and the only non-zero 
component of is = 2 in the Breit frame. Since Breit frame and proton rest frame 
are connected by a boost along the z-axis with boost factor 7 = Q/ [rripx) , the product 
of these vectors is {mn) = 27. 

Furthermore, two sets of physical polarization vectors, e(j) and e(j) (with i = 1,2) 
are defined by the conditions ek = ek = 0, = = —1, and em = en = 0. An explicit 
choice, written in light-cone co-ordinates, is 

e(i) = (0, , e(i)-L j and e(i) = ( ^ , 0, e(i)_L j , (C.12) 

where the transverse basis e(i)^ = (1, 0) and e(2)± = (0, 1) has been used. Note that the 
above equations hold in the proton rest frame, in the Breit frame, and in any other frame 
derived by a boost along the 2;-axis. 

These definitions give rise to the two following representations for the metric tensor: 
- (^A;) + (mk) {mky" ) ^^'^^^ 

The amplitude for the process in Fig. |C1| , with the lowest parton being a gluon in 
Feynman gauge, is proportional to 

A = e^{k') Vg{k', k),,T^H = Vy{k' , k)^,g^Pgp„T^H , (CIS) 

where V^'^ is the effective vertex for the scattering of the gluon off the external field. 
Next, the first and second metric tensor appearing in this expression for A are rewritten 
according to Eq. ( |C.13| ) and Eq. ( p.l4| ) respectively. In this situation, only the first terms 
from Eqs. (|C.13| ) and (|C.14|) contribute at leading order in x and A/Q. The intuitive 
reason for this is the relatively small virtuality of k, which ensures that for both the hard 
amplitude Th and the soft scattering vertex V only the appropriately defined transverse 
polarizations are important. 

To see this explicitly, consider the expression 

nk kn nn 



A = eVg 



, , mk km mm ,,9, 

[mk] [mk) [mky 



El 1/ IV IV I 1/ I L I L / 1 0\ 



H 



(C.16) 



where the appropriate contractions of vector indices are understood. Several estimates 
involving products of Vg and with specific polarization vectors will be required. 
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Note that both e and n are 0{1) in the Breit frame. For appropriate polarization, 
the amphtude (eT^/) involves no particular cancellation, i.e., it has its leading (formal) 
power behaviour in the dominant scale Q. Therefore, {nTn) is not enhanced with respect 
to (eTff), 

M^l^Od). (C.17) 

By analogy, it can be argued that {eVgm) is not enhanced with respect to (eV^e): since 
both e and m are 0{1) in the proton rest frame and (eV^e) has the leading power 
behaviour for appropriate polarizations e and e, the following estimate holds, 

(C.18) 



Gauge invariance requires {kTu) to vanish if k'^ = 0. Since the amplitude Th is 
dominated by hard momenta 0{Q), and /c^ ~ <^ Q^, this leads to the estimate 

(kT„) e- 

Analogously, from [eVgk) = at /c^ = 0, the suppression of this quantity at small 
virtualities k"^ can be derived. 

For this estimate it is also important that none of the soft scales involved in Vg, like k'^ 
or the gauge field ^4^, can appear in the denominator to compensate for the dimension 
of P. 

All the vector products nk, mk, ne, me, mn, and ee can be calculated explicitly. 
Using the relations in Eqs. ( |U.17 ) - ( |C20|) , it is now straightforward to show that 



[eVgc) (ee) {eTff) is indeed the leading term in Eq. ( P.16|) . The other terms are suppressed 
by powers of Q or k^. 

Thus, the leading contribution to \A\'^, with appropriate polarization summation 
understood, reads 

1^1' =E [(e{o(^')^se(i))(e(i)e(,))(e(,)TH)] [(e(i)(A;')V^^e(i,)) (e(i,)e(,v)) (e(,v)Tj^)] *, 

(C.21) 

where the arguments k and k of the polarization vectors in the first and second square 
bracket respectively have been suppressed. 

In the high-energy limit, the scattering of a transverse gluon off an external field is 
completely analogous to the scattering of a scalar or a spinor, 

ei'ip') Vg{p',p) ei{p) = 2n6{p'o - po) 2po 5n U^{p'i_ - P±) , (C.22) 

the only difference being the non-Abelian eikonal factor, which is now in the adjoint 
representation. 
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In analogy to the spinor case, the polarization sum decouples from the hard part for 
sufficiently inclusive measurements, so that the squared amplitude is proportional to 



Ul 



(C.23) 

Note that there is no summation over the index /. Recall the definition of B, the soft 
part of the amplitude square, given at the beginning of the last section and illustrated 
in Fig. |C.1| . The corresponding expression in the case of a soft gluon can now be read off 
from Eqs. (10221) and ([0231) : 



Bg = (e(i)(A;)e(i) (/;;)) (e(i)(A;)e(j)(Ai)) . 



(024) 



This is further evaluated using the explicit formulae in Eq. ( |012| ) and the identity 



6' 



ab 



(a, 6 G {1,2}) 



(025) 



Comparing the resulting expression. 



B, = kl 



6'^ + 



k? 



5'^ + 



2k\k{ f l-P ' 



k'l 



(026) 



to Eq. ( |01| ), the diffractive gluon distribution of Eq. (|4.54| ) is obtained. 

Note that the factor appearing in the denominator of Eq. ( [4.52| ) has been replaced 
by the dimension of the adjoint representation, iV^ — 1. 
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D Inclusive Parton Distributions 



This appendix is devoted to the calculation of inclusive parton distributions and inclu- 
sive structure functions in the semiclassical framework. Inclusive DIS was discussed in 
Sect. ^.2| , where the cross section for qq production off a colour field was derived, at the 
where the parton model interpretation of this cross section was outlined, 



end of Sect. 

and in Sect. [7.1| , where it was part of the semiclassical analysis of diffractive and inclusive 
structure function data. Here, the above scattered information is collected, and a more 
coherent account, including further technical details, is given (cf. |]62[| ). 

Recall the leading twist cross section ax for qq pair production by a transversely 
polarized photon obtained in Sect. ^]2| (cf. Eq.( p.27|) ). The corresponding transverse 
structure function reads 




il[dy,W,,{Q)dy,Wl^{Q) 



m + k"^. 



(D.l) 



To map this calculation onto the conventional parton model framework, identify the 
result as Ft^x^Q"^) = 2xq{x,Q'^). The corresponding quark distribution reads 



xq{x, Q' 




1 



J X I ^ 



Here the modified Bessel function Ki has been introduced so that, in both terms, the 
functions W appear in co-ordinate space. This makes it particularly clear that the first 
term is only sensitive to the short distance behaviour, while the second term depends 
on the non-perturbative long-distance structure of the colour field. Note that the sum of 
both terms is independent of /i^. 

The corresponding gluon distribution at small x is most easily calculated as 



xg{x,Q^ 



3n dFT{x,Q^ 
a. 



1 



(91nQ2 



27r2a, 



tr {dy^W,^iO)dy^W^^iO)) . (D.3) 



Equations ( p.2| ) and ( p.3| ) can serve as the starting point for a conventional partonic 
analysis of inclusive DIS. 
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To gain more physical insight into the correspondence of the semiclassical and the 
parton model approach, return to the starting point, Eq. (p.l| ). It is instructive to view 



Ft as a sum of two terms: F^^^^, the contribution of asymmetric configurations where 
quark or antiquark are slow, a < fi'^/Q'^ or 1—a < f^^/Q^ (Fig. [4.11| a), and F^^^, 
the contribution of symmetric configurations where both quark and antiquark are fast, 
a, 1 — a > /i^/Q^ (Fig. 4.11[ b). In a frame where the proton is fast, say, the Breit frame. 



the asymmetric and symmetric contribution to Ft correspond to photon-quark scattering 
and photon-gluon fusion respectively. 

The symmetric part is dominated by small qq pairs, i.e., by the short distance con- 
tribution to the Wilson-loop trace, 

/ tr {W,Ay±WUy±)) = Iv' [ tr {dy,W,AO)dy,WliO)) + 0{y') . (D.4) 

The corresponding contribution to the structure function is related to the first term on 
the r.h. side of Eq. ( p.2|) , which generates the gluon distribution, Eq. ( p.3|) . It can also 
be written as 

F^^^-(0, Q^) = ^,\l dzP,,{z) (in ^ - J^^ tr {dy,W^^^{0)dy,W^^^{0)) , (D.5) 

where Pqg{z) is the conventional gluon-quark splitting function. 

The other splitting functions appear if as corrections to Ft and Ft, associated with 
higher Fock states of the virtual photon, are considered in the semiclassical approach. 
For example, the qqg parton configuration involves, in the case where one of the quarks 
carries a small fraction of the photon momentum, a InQ^ term associated with Pqq{z). 

The splitting function Pgg{z) is most easily derived by considering an incoming virtual 
scalar which couples directly to the gluonic action term F^^F^'^ . Such a current was used 
previously in [ |191| ] to study small-x saturation effects. The relevant lowest order Fock 
state consists of two gluons. As reported in [^, it can be checked explicitly that the 



semiclassical calculation of the corresponding high-energy scattering process yields the 
usual gluon-gluon splitting function. 

Since the semiclassical approach exactly reproduces the well-known DGLAP split- 
ting functions, the large logarithms In (Q^/ /i^) can be resummed in the conventional way, 
by means of the renormalization group. To this end, the parton distributions q{x, Q"^) 
and g{x, Q^) are evaluated using DGLAP evolution equations, with the input distribu- 
tions q{x,Ql) and g{x,Ql) given by Eqs. ( p.2| ) and (p.3| ). Here Qo some small scale 
where logarithmic corrections are not yet important. The parton model description of 
the structure function at leading order includes only photon-quark scattering. The lead- 
ing logarithmic term from the photon-gluon fusion process appears now as part of the 
resummed quark distribution. 



The large hadron model of Sect, provides an expression for the basic function 
fx ^ tr (Wx^{y±)W^_^{y±)y Inserting the explicit formula of Eq. ( |6.25| ) into Eqs. ( p. 21 ) 
and ( |D.3D , the following compact expressions for the inclusive parton distributions at a 



low scale Qq are obtained 
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Figure |D.1| : The inclusive structure function F2{x,Q'^) at small x computed in the 
semiclassical approach with data from [|159|| . The data with = 1.5 GeV^ are not 
included in the fit. 



xq{x,Ql) 
xg{x,Ql) 



aQNc 
37r3 ^ 



In 



0.6424^ 



(D.6) 
(D.7) 



As discussed in Sect. [TJ] , the analysis of |^ introduces a soft energy dependence into 
these input distributions by ascribing a logarithmic growth to the total normalization, 



n ^n{L-\nx) 



(D.8) 



Note also the observation of ||192|| that the small-x structure function is well described by 
a simple ln(l/a;) with an additional InQ^ enhancement, similar to the effects of Altarelli- 
Parisi evolution. 

The DGLAP evolution of the above inclusive parton distributions and corresponding 
diffractive distributions given in Sect. |7.1| provides predictions for both the inclusive and 
diffractive structure functions. In the numerical analysis, the inclusive distributions are 
multiplied with (1 — x) to ensure vanishing of the distributions in the limit x ^ 1, which 
is required for the numerical stability of the DGLAP evolution. 
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A combined fit to small- a; HERA data gives a good description of the experimental 
results. While plots of ^2^^^"* are given Figs. |7.1| and |7.2| , here, corresponding results for 
the inclusive structure function are presented (see Fig. p.l| ). Since the underlying model 
is only valid in the small-x region, data points above x = 0.01 are not considered. To 
appreciate the quality of the fits, recall that, within the large hadron model, the diffractive 
(Eqs. (|7^) and ( |7.5|) ) and inclusive (Eqs. (p.6|) and (p.7|) ) parton distributions are highly 
constrained with respect to each other. 
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